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delay  time  ranges  from  about  2  to  100  ^s.  The  results  of  a  computer-modeling  calcu¬ 
lation  of  the  chemical  kinetics  suggest  that  the  photolytically  produced  hot  H 
atoms  and  vibrationally  excited  NH2  play  important  roles  in  promoting  the  ignition 
chenjistry. 

The  theory,  by  a  combination  of  numerical  and  approximate  analytic  methods, 
examines  the  interaction  of  premixture  stoichiometry,  preferential  diffusion  (i.e., 
differing  diffusivities  for  heat  and  species),  and  flame  radius,  with  respect  to 
the  outward  propagation  of  a  laminar  flame  from  a  single  spheroidal  kernel  (if  a 
flame  develops).  Previously  reported  observations  of  the  variation  of  the  flame- 
front  speed  with  flame  radius,  sometimes  interpreted  in  terms  of  preferential  - 
diffusion  effects  but  probably  better  interpreted  in  terms  of  radiative  heat  loss, 
have  been  examined;  the  need  exists  for  further  spherical -flame-propagation 
experiments,  with  emphasis  on  measurements  taken  when  the  flame  radius  is  more 
nearly  comparable  to  the  flame  thickness.  A 

The  theory  also  considers  the  configuration  of  the  (possibly  multiply  con¬ 
nected)  unburned-gas/burned-gas  interface  »tfvat  evolves  when  a  train  of  rapidly  re¬ 
peated,  brief  pulses  initiate  a  sequence  of  reacting  kernels  in  a  flowing  mixture. 
Such  nonintrusive  deposition  seems  a  means  of  stabilizing  a  combustion  wave  in  a 
fast-flowing  reactive  mixture  with  minimal  energy  and  power  requirements.  Many 
coplanar  sites  of  energy  deposition  would  be  required  to  span  (within  a  short  down¬ 
wind  distance)  a  supersonic  premixed  stream  with  fTame,  if  the  lateral  rate  of 
flame  propagation  from  each  irradiated  "blob"  is  typified  by  the  modest  deflagra¬ 
tion-wave  speed.  However,  if  the  nonintrusive  energy  deposition  is  rapid  enough 
and  substantial  enough  for  the  direct  initiation  of  a  detonation,  then  the  lateral 
rate  of  effective  flamespread  is  comparable  to  the  streamwise  supersonic  speed  of 
convection.  The  entire  breadth  of  the  premixed  stream  would  be  converted  to  pro¬ 
duct  within  a  relatively  very  short  distance  downwind  of  a  few  coplanar  sites  of 
nonintrusive  energy  deposition.  Thus,  the  direct  photochemical  initiation  of  deto¬ 
nation  seems  relevant  to  supersonic  combustion  with  minimized  entropy  increase. 

The  selfsimilar  flow  immediately  downwind  of  a  standing  oblique  (conical)  detona¬ 
tion  wave,  and  the  subsequent  expanding-nozzle  configuration,  have  been  discussed 
for  a  novel  supersonic  combustor. 
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1.0  RESEARCH  OBJECTIVES 


In  high-speed  air-breathing  combustion  systems,  a  premium  is  placed  on  the 
complete  release  of  chemical  energy  during  the  relatively  brief  residence  time  of 
reactants  within  combustors  of  practical  length:  mixing,  ignition,  and  chemical 
reaction  must  be  achieved  in  a  relatively  short  time.  We  address,  from  a  funda¬ 
mental  chemical -dynamics  and  fluid-transport  point  of  view,  the  use  of  both  al¬ 
ternative,  photochemical  ignition  (by  laser  irradiation)  and  ignition-promoting 
additives,  in  order  to  alter  chemical -reaction  pathways  to  achieve  rapid  ignition 
and  enhanced  combustion  rate.  We  address  the  influence  of  mixture  inhomogeneity, 
and  of  departure  from  stoichiometric  proportion  in  fuel/air  mixtures,  on  the  pro¬ 
cesses  of  ignition,  flame  development,  and  flame  propagation.  Chemical  systems 
of  particular  interest  include  hydrogen/air  and  methane/air,  often  with  trace 
amounts  of  carefully  selected  sensitizers.  We  identify,  by  irradiating  premix¬ 
tures  flowing  faster  than  the  adiabatic  flame  speed,  and  by  carrying  out  support¬ 
ing  approximate  analyses  of  these  experiments,  optimal  circumstances  [minimal 
input  energy,  minimal  amount  of  sensitizer(s),  etc.]  for  achieving  ignition  and 
burnup  with  currently  available  optical  sources. 

A  far  more  detailed  exposition  of  the  goals,  background,  and  relevant  litera¬ 
ture  may  be  found  in  the  two  earlier  annual  reports  issued  on  this  project 
(Fendell,  Chou,  Zukowski,  and  Carrier  1988,  1989).  In  the  interest  of  brevity, 
for  the  discussion  which  follows,  the  reader  is  taken  to  be  aware  of  the  contents 
of  those  earlier  reports. 
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2.0  ANALYSIS 


The  aerodynamics  associated  with  flame  stabilization  in  high-speed  reactive 
flows,  for  propulsion  objectives,  has  been  the  subject  of  the  analytical  efforts 
undertaken  in  this  project.  The  novel  flame-stabilization  mechanism  under  study 
is  based  on  the  efficient  absorption  of  energy  nonintrusivel v  deposited  into  com¬ 
bustible  gaseous  mixtures  by  use  of  modern  optical  sources  (lasers).  The  experi¬ 
ments  carried  out  with  lasers  in  this  project  demonstrate  the  existing  capability 
to  initiate  deflagration  waves;  the  complete  international  literature  on  the  use 
of  lasers  operating  at  many  different  frequencies  to  achieve  ignition  in  sensi¬ 
tized  or  unsensitized  mixtures  of  various  reactants  over  a  wide  range  of  stoich¬ 
iometries  is  now  extremely  voluminous.  Further,  the  direct  photochemical  initia¬ 
tion  of  detonations  in  highly  subatmospheric-pressure  gaseous  mixtures  composed 
of  various  constituents  was  first  reported  a  dozen  years  ago  by  use  of  xenon 
flashlamps  (Lee  et  al .  1978).  Here  we  attempt  to  elucidate  by  analysis  a  few  of 
the  new  flame-stabilization  possibilities  (for  high-speed  air-breathing  propul¬ 
sion)  that  may  become  feasible  as  laser  technology  evolves.  In  particular,  we 
anticipate  that  very  rapidly  repeatedly  pulsed,  energetic  lasers  that  are  both 
compact  and  reliable  may  be  achieved  within  a  couple  of  decades. 

Attention  was  initially  concentrated  on  the  continuous  nonintrusive  deposi¬ 
tion  of  laser  radiation,  efficiently  absorbed  by  a  uniform  reactive  mixture. 

The  mixture  is  taken  to  be  flowing  very  subsonically,  but  still  supercritically, 
i.e.,  still  faster  than  the  deflagration-wave  speed,  so  that  a  special  provision 
such  as  energy  deposition  is  required  to  preclude  flame  blowoff.  The  situation 
is  the  counterpart  of  the  more  familiar  use  (e.g.,  in  a  flat-flame  burner)  of 
energy  extraction  to  preclude  flashback  in  a  subcritically  flowing  mixture.  With 
energy  deposition,  the  burned-gas  temperature  is  in  general  in  excess  of  the  adia¬ 
batic  flame  temperature,  whereas,  with  energy  extraction,  the  burned-gas  tempera¬ 
ture  is  in  general  less  than  the  adiabatic  flame  temperature.  The  early  studies 
in  this  project  concentrated  on  planar,  line,  and  point  radiant-energy  sources, 
such  that,  respectively,  one-dimensional,  two-dimensional,  and  axisymmetric 
flames  were  stabilized  (Carrier  et  al .  1990;  Fendell  et  al.  1990).  Only  in  the 
planar  case  is  the  solution  analytically  tractable  for  the  large-heat-release 
conditions  characteristic  of  combustion;  for  the  line  and  point  sources,  analytic 
solution  is  tractable  only  for  the  limit  of  vanishingly  small  heat  release, 
because  large  chemical  exothermicity  introduces  aspects  to  the  dynamics  that 
necessitate  numerical  treatment.  Such  investment  in  numerical  solution  may  not 
be  warranted  because  stabilization  by  a  continuous  nonintrusive  source  is  not 
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likely  to  be  practically  exploited,  though  such  a  demonstration  in  the  laboratory 
may  soon  be  feasible,  and  would  mark  a  noteworthy,  instructive  milestone.  The 
practical  difficulty  with  a  continuous  deposition  for  flame  stabilization  is  that 
the  deposition  is  into  burned  mixture  (since  the  flame  is  situated  upwind  of  the 
deposition  site),  and  needs  to  be  relatively  large  (since  every  element  of  near¬ 
source  mass  throughput  must  be  irradiated  to  preclude  blowoff).  For  a  planar 
flame,  this  means  the  energy  deposition  must  be  achieved  uniformly  across  the 
entire  deposition  plane;  for  a  line  or  point  source,  only  that  mass  throughput 
passing  near  the  source  need  be  augmented  in  enthalpy,  but  the  energy-deposition 
requirement  is  still  demanding  relative  to  that  for  a  rapidly  repeatedly  pulsed 
source. 

For  the  intermittent  source,  only  enough  radiant  energy  need  be  deposited  to 
initiate  discrete  flame  kernels;  as  each  flame  kernel  is  convected  downwind  by 
the  flow,  the  flame  evolving  from  the  kernel  propagates  toward  the  flame  evolving 
from  the  prior  kernel  and  the  flame  evolving  from  the  subsequent  kernel,  until 
the  intervening  reactive  mixture  is  converted  to  product  at  some  downwind  dis¬ 
tance.  The  flame  of  ea»_h  kernel,  then,  envelopes  already-burned  mixture,  and 
propagates  into  unburned  mixture  (see  Appendices  1  and  2).1  Since  the  same 
element  of  mixture  cannot  be  burned  twice,  the  initially  multiply-connected 


^The  interplay  of  stoichiometry,  differing  diffusivities,  and  flame  curvature 
in  the  early  radially-outward  propagation  of  a  single  spherical  flame  kernel 
in  an  unbounded  expanse  of  combustible  mixture  has  been  the  subject  of  experi¬ 
mental  and  theoretical  inquiry  over  the  last  decade,  but  misinterpretations 
(Frankel  and  Sivashinsky  1983;  Strehlow  1984)  may  have  arisen.  It  is  true 
that,  if  a  reactant  is  both  stoichiometrically  deficient  and  relatively  fast- 
diffusing,  then  theoretically  the  speed  of  a  spherical  flamefront  may  be  great¬ 
er  at  a  smaller  flamefront  radius,  at  which  curvature  effects  are  significant, 
than  at  a  larger  flamefront  radius,  at  which  the  propagation  is  effectively 
planar.  What  is  less  clear  is  whether  such  an  effect  occurs  for  parameter 
values  arising  in  practice;  what  is  even  more  doubtful  is  whether  certain  ex¬ 
periments  (Palm-Leis  1966;  Palm-Leis  and  Strehlow  1969)  are  evidence  of  such  an 
effect,  since  the  data  are  taken  at  flame  radii  much  larger  than  the  flame 
thickness;  furthermore,  the  effect  is  observed  only  at  sufficiently  fuel -rich 
stoichiometries  in  propane-air  mixtures  (such  that  radiative  heat  loss  seems  an 
alternative  explanation  of  the  decrease  of  flamefront  speed  with  increasing 
flamefront  radius),  and  is  not  observed  at  any  stoichiometry  in  methane-air  mix¬ 
tures  (which  have  a  smaller  propensity  to  produce  soot).  At  any  stoichiometry 
tested,  the  observed  flamefront  speeds  in  the  cited  spherical -fl ame  experiments 
either  fail  to  approach  an  asymptote  even  at  8-9-cm  radius,  or  approach  an 
asymptote  inconsistent  with  the  density-ratio-adjusted  laminar-planar-flame- 
speed  value.  There  appears  to  remain  a  need  for  definitive  spherical -fl ame- 
propagation  experiments,  with  emphasis  on  data  taken  at  small  flamefront  radii, 
in  order  to  elucidate  unresolved  issues  concerning  flame-kernel  development. 
These  suggested  combustion  experiments  are  not  concerned  with  the  exotic  photo¬ 
chemistry  of  energy  deposition;  in  fact,  ideally,  no  residuals  of  the  ignition 
process  would  remain  even  at  very  early  times  in  the  spherical -flame-kernel 
evolution. 
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flames  merge  to  form  a  "scalloped"  reacting  front  circumscribed  by  the  flame 
locus  engendered  by  a  continuous  line  or  point  source  depositing  energy  at  the 
same  site  as  the  intermittent  (say,  periodic)  point  or  line  source.  Even  if  some 
of  the  kernels  do  not  spawn  flames,  eventually  all  the  Intervening  mixture  is 
converted  to  product  --  but  the  conversion  occurs  further  downwind.  However,  if 
the  convective  speed  of  the  flowing  reactants  far  exceeds  the  propagation  speed 
of  the  premixture  (and,  at  one  atmosphere,  the  laminar  propagation  speed  of  even 
a  stoichiometric  hydrogen/air  mixture  is  but  about  2  m/s),  then  significant 
lateral  spread  to  span  a  stream  with  flame  via  deflagration  occurs  only  far,  far 
downwind.  One  might  seek  to  deposit  pulses  at  several  lateral  sites  to  shorten 
the  streamwise  distance  for  reactant  conversion  to  product;  this  procedure  not 
only  Incurs  more  complication  with  the  source  technology,  but  also  may  incur 
large  pressure  loss  and  high  entropy  rise  owing  Is  thl  Possible  Initiation  of  a 
strong  (normal )  detonation! 

Thus,  the  most  viable  procedure  to  span  a  supersonic  stream  of  reactant  with 
flame  within  a  short  streamwise  distance  of  a  site  of  nonintrusive  energy  deposi¬ 
tion  would  appear  to  be  the  direct  initiation  of  a  detonation.  As  with  the  stabi 
lization  of  a  deflagration  wave  in  a  stream  of  mixture  flowing  at  a  modestly 
supercritical  speed,  so  with  the  stabilization  of  a  detonation  in  a  stream  of 
mixture  flowing  supersonically:  the  least  energy  and  power  requirement  is  associ¬ 
ated  with  a  very  rapidly  repeated  source,  rather  than  with  a  continuous  source 
(which  might  lead  to  a  locally  normal  detonation).  Also,  analogously,  the  spheri 
cal  detonations  initiated  by  a  long  train  of  pulses  interact  to  form  a  scalloped 
detonation  wave.  As  the  interval  between  laser  pulses  vanishes,  the  pressure 
losses  owing  to  the  reflected  shocks  that  arise  from  the  interaction  of  neighbor¬ 
ing  spherical  detonations  also  vanish,  and  a  nearly  conical  detonation  wave  is 
formed.  The  word  nearly  is  introduced  in  that  a  cellular  structure  arises  in 
connection  with  spherical  detonations  (Lee  1977,  1984);  however,  the  finite  struc 
ture  of  this  conical  detonation  seems  a  detail  of  minor  practical  consequence. 

The  flow  downwind  of  this  (weak,  oblique)  conical  detonation  may  be  modeled  as 
steady  supersonic  inviscid  and  axlsymmetrlc;  in  fact,  it  is  also  selfsimilar  and 
readily  described  In  quantitative  detail  (see  Appendix  3). 

The  concept  of  basing  a  supersonic  combustor  (i.e.,  a  combustor  in  which  the 
flow  is  everywhere  supersonic)  on  a  standing  (i.e.,  stabilized)  oblique  detona¬ 
tion  wave  was  first  proposed  decades  ago  (e.g..  Gross  and  Chlnitz  1960).  The 
concept  seems  a  viable  alternative  to  the  mixing-controlled  supersonic-diffusion- 
flame  designs  that  have  received  the  preponderance  of  attention  from  supersonic- 
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combustor  designers.  However,  previously,  emphasis  has  been  placed  on  the  intru¬ 
sion  into  the  reacting  stream  of  a  symmetrically  positioned  right  circular  cone 
of  finite  length.  The  potential  difficulties  with  such  a  design  (or  its  two- 
dimensional  analogue,  the  intrusion  of  a  symmetrically  positioned  wedge  of  finite 
length)  are  many:  additional  frictional  drag  is  incurred;  all  such  bodies  in 
practice  will  be  somewhat  blunted  and  subject  to  jitter,  such  that  an  at  least 
locally  normal  detonation  occurs,  with  associated  pressure  losses;  and,  for  the 
high-altitude-flight  conditions  of  interest,  the  pressure  achieved  even  in  the 
weakly  shocked  flow  may  be  low  enough  that,  in  the  absence  of  a  positive  ignition 
device,  induction  times  for  reaction  may  result  in  impracticable  long  combus¬ 
tors.  Even  if  a  conical  detonation  wave  is  supported  by  a  pointed  cone,  the  wave 
is  likely  to  incur  a  larger  entropy  increase  than  the  increase  associated  with  a 
pulsed-laser-induced  wave.  Whenever  a  laser  with  the  energy,  power,  and  repeti¬ 
tion  rate  is  available,  it  seems  clearly  a  superior  alternative  to  the  introduc¬ 
tion  of  an  intrusive  cone. 

Downwind  of  the  axial  plane  at  which  the  conical  detonation  wave  intersects 
the  surrounding  container,  the  selfsimilar  nature  of  the  flow  in  the  shocked  and 
reacted  gas  no  longer  holds;  that  is,  information  is  passed  through  the  post¬ 
detonation  flow  of  a  finite  physical  length  (the  radius  of  the  combustor  wall). 

In  fact,  the  container  downwind  of  the  detonation-wave/combustor-wall  interaction 
must  flare  appropriately  to  accommodate  an  expansion:  too  great  a  flaring  results 
in  flow  separation  and  effectively  no  thrust,  and  too  little  a  flaring  results 
in  compressional  shocks,  regions  of  locally  subsonic  flow,  and  large  pressure 
losses.  Detailed  refined  design  of  a  nozzle  (to  expand  efficiently  the  shocked 
and  reacted  gas  from  the  high-pressure  levels  of  the  selfsimilar  flow  to  a  flow 
that  is  laterally  uniform  at  the  ambient  pressure)  seems  not  properly  the  func¬ 
tion  of  a  basic  inquiry  into  supersonic  combustors.  Nevertheless,  it  seems  cru¬ 
cial  to  establish  that  one  can  design  a  short  nozzle  that  can  convert  high-pres¬ 
sure  conditions  to  high-flow-speed  conditions,  such  that  the  associated  thrust 
far  exceeds  the  pressure  drag  incurred  by  the  ambient  streaming  past  the  combus¬ 
tor.  [At  the  hypersonic  conditions  of  interest,  the  pressure  drag  does  far 
exceed  the  frictional  drag,  and  the  pressure  drag  can  be  estimated  simply  but 
adequately  by  the  Newtonian  approximation  (no  rebounding  of  ambient-flow  mole¬ 
cules  that  make  contact  with  the  enveloping  combustor  container)].  No  proof  of 
principle  of  the  feasibility  of  a  laser-initiated  oblique-detonation-wave  ergine 
seems  complete  without  evidence  that  such  a  properly  flared  nozzle  can  be  identi¬ 
fied.  Thus,  this  inverse  nozzle-flow  exercise  (given  the  "entrance"  and  "exit" 
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conditions,  design  a  nozzle,  so  the  flow  can  undergo  the  transition  between  the 
conditions  within  a  short  streamwise  distance)  is  pursued  by  trail -and-error 
numerical  calculations  with  candidate  container  configurations. 

The  laser-initiated  oblique-detonation-wave  engine  presumes  the  prior  forma¬ 
tion  of  a  fuel -air  mixture.  However,  while  there  must  be  sufficiently  large  frac 
tion  of  the  flowing  gaseous  mixture  that  is  detonable,  there  is  no  requirement 
that  the  entire  mixture  be  so  constituted.  In  fact,  we  have  speculated  that 
there  may  be  performance  advantages  if  one  could  so  tailor  the  flowing  mixture 
that  a  fraction  would  consist  of  mill imeter-sized  blobs  of  unmixed  fuel  and  air 
in  stoichiometric  proportion.  Such  blobs  would  be  compresseu  but  unburned  during 
passage  through  the  conical  detonation,  in  which  the  well-mixed  portion  of  the 
mixture  is  converted  to  hot  product  gas.  Downwind  of  the  detonation  the  blobs 
would  be  converted  in  the  hot  surrounding  gas,  by  a  diffusionally  controlled 
reaction. 

On  the  basis  of  the  above  line  of  inquiry,  we  believe  th:t  one  remaining 
task  of  high  priority  concerns  the  efficient  mixing,  without  reaction  or  strong 
shocks,  of  fuels  of  interest  (such  as  hydrogen)  into  a  supersonic  airstream;  such 
a  task  is  probably  most  productively  pursued  experimentally.  Another  task  of 
high  priority  concerns  systematic  laboratory  experimentation  to  establish  what 
amount  of  energy,  what  time  of  deposition,  and  what  focal  volume  permits  the 
direct  photochemical  initiation  of  detonation  in  a  (stagnant)  mixture  of  a  compo¬ 
sition,  and  at  a  thermodynamic  state,  of  interest.  Also,  the  interaction  of  two 
such  spherical  detonations,  set  off  in  closr  proximity  in  space  and  time  (but  not 
precisely  coincidentally  or  simultaneously)  in  a  stagnant  mixture,  should  be  ex¬ 
perimentally  examined  to  be  certain  that  no  overlooked,  large-entropy-producing 
phenomena  would  occur  to  frustrate  the  achievement  of  small -pressure-loss, 
low-entropy-rise  supersonic  combustion  via  the  above-sketched  laser-initiated 
oblique-detonation-wave  engine.  Further,  for  completeness,  the  interaction  of  a 
spherical  detonation  wave  with  a  container  wall  might  be  worth  laboratory  examina 
tion  to  be  certain  that  no  here-unanticipated  complications  arise  in  connection 
with  the  accompanying  noise,  vibration,  or  rejected  shock  --  though  we  propose 
to  flare  the  container  wall  to  avoid  detonation-wave/container-wall  interaction. 
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3.0  EXPERIMENT 


The  results  achieved  under  the  experimental  task  are  reported  in  Appendix  IV. 
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APPENDIX  I. 


Evaluating  a  Simple  Model  for  Laminar-Flame-Propagation  Rates. 

I.  Planar  Geometry 
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A  comparison  with  experimental  flame-speed  data  is  made  for  results  from  a 
simple  formulation,  conveniently  reducible  to  quadrature,  of  laminar  isobaric 
flame  propagation  in  an  initially  homogeneous  gaseous  mixture.  A  direct  one-step 
irreversible  bimolecular  second-order  chemical  reaction  with  large  Arrhenius 
activation  energy  Is  adopted,  but  account  is  taken  of  the  modified  exothermicity 
owing  to  partial  dissociation  of  the  product  species  and  to  other  causes  of  incom¬ 
plete  oxidation.  The  effects  arising  from  differing  diffusivitles  for  heat  trans¬ 
fer  and  reactant-species  transport  are  developed  to  within  the  limitations  of  the 
model.  A  tractable  general  expression  is  obtained  for  the  steady-laminar-flame- 
propagation  speed,  by  exploitation  of  the  two-zone  (convective-diffusive,  diffu¬ 
sive-reactive)  deflagration-wave  structure.  For  simple-fuel/air  mixtures,  the 
predicted  variation  of  flame  speed  with  equivalence  ratio  4  agrees  fairly  well 
with  experimental  data,  with  the  provision  that  (1)  the  cold-mixture  transport 
properties  are  evaluated  for  the  pertinent  value  of  4,  and  (2)  the  overall 
activation  energy,  taken  to  be  Invariant  with  4,  is  ascribed  values  In  the  10-15 
kcal/mole  range.  In  particular,  the  experimentally  observed  off-stoichiometric 
(often  fuel-rich)  condition  for  peak  flame  speed  in  simple-fuel/air  premixtures 
is  recovered.  However,  for  simple-fuel/oxygen  mixtures,  the  predictions  deviate 
from  the  data,  although  the  sense  of  the  stoichiometry  for  which  the  peak  speed 
occurs  is  recovered.  Finally,  observations  concerning  the  transient  planar 
deflagration  wave  are  made,  as  background  useful  for  a  study  of  spherical  flame 
propagation,  to  be  reported  separately. 


1.  INTRODUCTION 

1.1  Objectives 

This  work  stems  from  an  interest  in  the  multi parameter  dependence  of  the 
flame  speed  of  a  spherical  flame  on  its  radius,  on  the  transport  properties  of 
the  reactants,  and  on  the  stoichiometry  of  the  combustible  mixture  (Strehlow 
1984).  For  the  ease  of  interpretation  of  results,  In  view  of  this  multiparameter 
dependence,  it  Is  advantageous  first  to  present  an  analysis  of  (planar)  one-dimen¬ 
sional  flames;  spherically  propagating  flames  are  then  treated  in  a  companion 
study  (Carrier  et  al.  1990). 

Ve  adopt  a  basic,  well-established  but  approximate  formulation.  It  does  not 
serve  our  purpose  to  dwell  at  length  on  Identifying  the  most  general  formulation 
tractable  for  a  reduction  to  quadratures,  the  form  to  which  we  seek  to  express 
the  solution;  thus,  the  model  could  be  generalized  at  several  turns.  Ve  note 
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that,  for  the  convenience  of  nondimensionalization,  we  refer  temperature,  den¬ 
sity,  and  several  parameters  to  their  ambient  (unburned)  values,  because  w£  are 
interested  ijj  a  comparison  flf  ihg  predicted  flame-propagation  speed  with  families 
of  experiments  wherein  the  individual  members  of  each  family  involve  a  (mostly) 
common  initial  (unburnt)  thermodynamic  state. 

We  assign  parameters  observed  or  plausible  values  and  then  evaluate  the  pre¬ 
dictive  accuracy  for  simple-fuel/oxidizer  mixtures.  For  example,  we  do  not 
assign  the  orders  of  reaction  to  have  fractional,  even  negative  values,  In  order 
to  try  to  match  flammability-limit  data  (an  adiabatic  model  inherently  cannot  pre¬ 
dict  such  limits  —  cf.  Westbrook  and  Dryer  1981).  Neither  do  we  assign  the  acti¬ 
vation  energy  a  variation  with  equivalence  ratio  to  curve  fit  the  results  of  an 
equal -diffusion-coefficient  model  to  flame-propagation-speed  data  (cf.  Coffee  et 
al.  1983,1984). 

1.2  Formulation 

We  adopt  a  direct  one-step  bimolecular  second-order  irreversible  chemical 
reaction.  The  entirely  gaseous  product  species  are  formed  by  burning  a  pure 
gaseous  fuel  F  (e.g.,  hydrogen  or  a  simple  hydrocarbon)  either  in  oxygen  or  in 
air  (while  nitrogen  could  be  replaced  by  some  other  inert  gas  for  purposes  of 
inquiry,  we  develop  In  terms  of  nitrogen) 

"F  F  +  *'o  0  "  product(s)  +  Q,  (1.1) 

where  Q  is  the  heat  evolved  from  burning  with  oxygen  0,  with  denoting  the 
stoichiometric  coefficient  of  species  i  for  complete  oxidation  of  the  fuel  to 
C02  and  H20.  In  fact,  the  products  vary  appreciably  as  conditions  change 
(e.g.,  from  significantly  fuel-rich  to  significantly  fuel-lean),  so  the  value  of 
Q  is  permitted  to  be  stolchiometrlcally  dependent.  The  heat  Q  Is  assigned  to 
reflect  the  exothermiclty  consistent  with  an  adequate  accounting  of  the  species 
actually  present  In  the  burned  mixture  at  equilibrium  (Coffee  et  al.  1983).  If 
■  •  "to  •'O  +  "F  *T*  where  mi  *s  the  mo1ecular  weight  of  species  1,  and  uref 
Is  a  reference  speed,  to  be  Identified  later, 

-»*  * 

X  “  T  *  1  *  F  »  *here  j2  "  *u  *o  ’  lo  "  "u/uref  ;  0-*) 

o 

_  * 

,3  -  -H—  ,  where  p  e  -  1  ;  (1.3) 

u  'u  ref 
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Q  YF  m  ?c  Q  Y0  m  70 

-  an  r  T  »  Yr  -  = — —  ;  X  -  — -V  .  Vn  -  — . 

HI  C„  Tlt  r  HI*-  v r  ID  C  T..  0  nu 

pu  r  r  p  U  0  U 

Here,  a  superasterisk  Is  used  to  designate  a  dimensional  quantity  if  there  other¬ 
wise  would  arise  an  ambiguity  with  its  dimensionless  counterpart.  The  dimen¬ 
sional  spatial  coordinates  are  designated  x  ,  although  below  we  shall  restrict 
attention  to  one-velocity-component  flows  spatially  dependent  on  only  one  spatial 
coordinate  (designated  x*  in  a  planar  flow,  r*  in  a  spherically  symmetric 
flow);  the  reference  length  scale  i  is  defined  In  terms  of  the  reference  thermal 
dlffusivity  «u  and  a  reference  speed  ureP,  to  be  assigned  below.  The  thermal  dif- 
fuslvity  «u  Is  defined  (as  usual)  In  terms  of  the  thermal  conductivity  ku,  the 
density  pu,  and  the  specific  heat  capacity  at  constant  pressure  cp,  where  sub¬ 
script  u  refers  to  the  (uniform)  state  of  the  cold  (unburned)  premixture  and  cp 
Is  approximated  as  a  constant  throughout  the  flow.  If  we  anticipate  that,  for  a 
given  "pair"  of  reactants  (such  as  a  pure  hydrocarbon  and  air,  the  quantity  ureP 
is  assigned  a  value  which  does  not  vary  with  the  mixture  composition,  then  we 
note  that,  since  *u  does  vary  with  the  mixture  composition,  the  length  scale  and 
the  time  scale  adopted  In  (1.2)  vary  with  mixture  composition.  Also,  T  signifies 
temperature  and  p,  density,  and  both  are  normalized  against  their  values  in  the 
cold  premixture.  The  molecular  weights  of  the  dominant  species  (In  terms  of  mass 
fraction)  are  taken  to  be  comparable  throughout  the  flow,  so  that  the  equation  of 

state  for  the  highly  subsonic,  effectively  isobaric  flows  of  interest  here  Is 

taken  in  (1.3)  to  be  p  e  -  1.  If  Y^  denotes  the  mass  fraction  of  species  i,  (1.4) 
reflects  the  standard  practice  of  Introducing  the  stoichiometrlcally  adjusted 
mass  fractions  Yp,  equal  to  (m  Yp)/(mp  i/p),  and  YQ,  equal  to  (m  Y0)/(m0  j/q); 

for  a  fuel -lean  condition,  the  excess  oxygen  In  the  burned  gas  YQb  ■  YQu  -  Ypu  - 

Y0u(l  -  4),  where  the  fuel-to-oxldlzer  equivalence  ratio  4>  •  YFu/You*  YFu  " 

<«  fFu)/(»F  *F)»  Y0u  -  (m  Y0u)/(m0  *0).  Equivalently,  *  -  (YFu/Y0u)/r,  r  - 
(mp  •'pl/fag  •'o)’  *qu*valently»  If  yn  Is  the  mass  fraction  of  nitrogen,  4  » 

EYFi/(YOu  +  ?Nu,l/(r  N)t  s1nc*  [Y0u/{Y0u  +  YNu}]  *  N»  a  consM«  0.23  for  air) 

—  so  4  also  can  be  defined  In  terms  of  the  fuel-to-alr  ratio.  For  simplicity 

here,  attention  is  limited  to  cold  premixtures  consisting  only  of  fuel,  oxygen, 
and  (sometimes)  nitrogen,  so  YOu  +  YNu  +  YFu  “  1  ;  further,  attention  Is  limited  to 
the  cases  H  •  0.23  (pure-fuel/alr  mixture)  and  N  -  1  (pure-fuel/oxygen  mixture). 

In  (1.4)  we  also  choose  to  Introduce  the  energetically  adjusted  mass  frac¬ 
tion  for  fuel  Y  and  for  oxygen  X;  these  are  so  defined  that  the  Increase  of  one 
unit  of  t  Is  effected  by  the  decrease  of  one  unit  of  (say)  Y,  with  the  understand¬ 
ing  that  the  decrease  is  owing  to  chemical  conversion  In  accord  with  the  reaction 
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(1.1).  In  particular,  If  Yu  -  (Q  Ypu)/(m  cp  Tu)  ■  a,  then  (for  a  fuel -lean 
situation,  1  >  4  >  0),  in  the  fully  burned  ("hot")  gas,  In  which  Y  ■  0,  e  would 
have  increased  from  a  value  of  unity  In  the  unburned  gas  to  a  value  of  (1  +  a)  in 
the  burned  gas,  In  the  absence  of  heat  losses.  The  quantity  (1  +  a)  Is  termed 
the  adiabatic  flame  temperature  if,  and  a,  sometimes  termed  the  second  Damkohler 
number  and  given  the  symbol  Og,  Is  the  ratio  of  chemical  exothermicity  per  mass 
of  fuel  (Q  Ypy/m)  to  the  cold  premixture  enthalpy  cp  Tu.  In  the  absence  of  a 
significant  occurrence  of  soot  (so  that,  with  reservations  discussed  in  a  spherical 
flame  context,  the  radiative  transfer  of  heat  Is  negligible  for  most  purposes,  in 
accounting  for  the  conservation  of  energy),  this  formulation  holds  for  the  case  of 
a  fuel -rich  mixture  by  noting  that  in  the  burned  gas,  in  which  now  X  -  0  and  Y  ■ 

(a  -  b)  >  0,  If  -  (1  ♦  b),  where  Xu  »  (Q  Y0u)/(m  cp  Tu)  -  b.  The  stoichio¬ 
metric  case  (4  »  1)  Is  a  »  b  -  a+,  where  henceforth  the  superscript  +  denotes 
evaluation  at  stoichiometric  conditions.  Note  that  a*  may  take  on  a  range  of 
values  for  the  same  pair  of  reactants  with  Invariant  cold-premlxture  pressure  and 
temperature,  because  of  variable  dilution  with  chemically  inert  species. 

For  the  transport  coefficients,  if  D*  is  the  (binary)  diffusion  coefficient 
of  species  i,  and  k*  is  the  thermal  conductivity,  then  It  is  useful  to  introduce 
dimensionless  counterparts  (<^  is  the  Lewis-Semenov  number  of  species  1): 


*  1 
D1 

M 

ut; 

iuj 

1  u 

Di  *i»  1  *  °’F> 


(1.5) 


* 

where  D^u  is  the  (dimensional)  value  of  in  the  cold  premixture,  and  0.  is 
dimensionless; 


k_ 
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(1.6) 


where  (again)  ku  Is  the  (dimensional)  value  of  k*  in  the  cold  premixture,  and  k 
is  dimensionless.  The  transport  coefficients  and  k  are  taken  to  be  functions  of 
#  only. 

For  the  chemical  kinetics,  a  direct  one-step  irreversible  blmolecular 
second-order  reaction,  subject  to  an  Arrhenius-type  specific  rate  constant 
(adjusted  to  obviate  the  cold-boundary  difficulty).  Specifically,  the  specific 
rate  constant  is  taken  in  the  form 


T°  exp 


T. 

m: 


Bf  Tj  #°  exp[- 


(1.7) 
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where  Ta  is  activation  temperature;  the  effective  frequency  factor  Bf  has  units 
of  n3/(s  nole3  K°);  and  the  dimensionless  activation  energy  E  •  (Ta/Tu),  a 
large  quantity  for  premixtures  of  interest  here. 


2.  THE  ONE-DIMENSIONAL  ANALYSIS 


He  consider  the  one-dimensional  structure  of  a  steady  flame  that  sits  at 
x  •  0  In  a  coordinate  system  In  which  the  flow  is  in  the  positive  x  direction.  We 
consider  only  isobarlc  phenomena  so  that  the  conservation  of  momentum  is  not  at 
Issue;  mass  conservation  requires  merely  that  p  u  -  u(-«)  -  const,  and,  as  is  well 
known,  p  u  will  be  an  eigenvalue  of  the  mathematical  problem. 

The  equations  describing  the  conservation  of  energy  and  of  chemical  species 
are,  in  accord  with  the  discussion  of  Section  1, 


'  “  'x  •  tk(,)  'x5x '  -T-  •* '2  ¥  *  exp[-  rh] 

uref 


YJx  *  -#a'2  RR  » 

l*'1  D0(#)  Xjx  -  -f°*2  RR  . 


(2.1) 

(2.2) 

(2.3) 


Henceforth,  we  adopt  a  common  temperature  dependence  of  the  reactant-species 
diffusion  coefficients,  conditions  consistent  in  accuracy  with  previous  approxi¬ 
mations: 

k (I)  •  r1  Dp(l)  -  r1  D0(O  .  (2.4a) 

The  Lewls-Semenov  numbers  cQ,cf  are  known  input  quantities.  The  numerator  of 
the  first  Damkohler  number  (J  /u|;  p)  is  defined  as  follows: 

J*  "  "u  Bf  'u  Tu+1  v</(Q/V  *  (*-4b) 


The  equations  (2.1)- (2.3)  Involve  dimensionless  quantities  only,  except  for  the 
factors  J*  and  u2gp. 

The  foregoing  equations  can  be  reduced  without  further  approximation  to 
their  constant-diffuslvlty  forms  by  use  of  the  definition 

(  •  Jo  WO]"1  <ix.  (2.4c) 

With  that  substitution,  (2.1)  -  (2.3)  become 


1-5 


p  u  9^  -  »  o(i)  RR  ,  0(9)  ■  9a  ^  k($)  ; 

(2.5) 

p  u  Y?  -  oF  Y^  -  -o(9)  RR  , 

(2.6) 

p  u  X^  -  *o  *££  *  -»(*)  RR  • 

(2.7) 

As  noted  earlier,  we  confine  our  attention  to  mixtures  for  which  E  is  large 
enough  to  imply  that  the  thickness  of  the  zone  in  which  reactions  cannot  be 
Ignored,  (*>uE)‘l,  is  small  compared  to  the  thickness,  (/>u)'*,  of  the  diffusive- 
convective  (preheating)  zone  (Bush  and  Fendell  1970;  Fendell  1972).  It  is  im¬ 
plied  in  all  large-activation-energy  treatments  of  gaseous  deflagrations  that  the 
temperature  to  be  assigned  to  the  hot  end  of  the  convective-diffusive  zone  is  the 
same  as  that  assigned  to  the  hot  end  of  the  reaction  zone.  With  this  constraint 
we  note  that,  both  in  the  reaction  zone  and  in  the  downstream  zone  where  9  •  gf 
»  const.,  (2.5) - (2-7)  imply 


(a  +  cf  Y)^  -  0  ,  (2.8) 

{»  +  CQ  X)u  -  0  .  (2.9) 

When  the  fuel  content  of  the  cold  premixture  Yu(-  a,  for  brevity)  is  smal¬ 
ler  than  the  oxygen  content  of  the  cold  premixture  Xu(»  b),  it  is  recalled 
that  the  burnt  state  of  the  gas  is  one  for  which  Y(«)  »  0,  X(«)  «  (b  -  a), 

#(»)  »  6f  «  (1  +  a).  Accordingly,  (2.8)  and  (2.9)  imply 

Y  •  (1  +  a  -  9)/of  ,  (2.10) 

X-b-a+(l  +  a-  f)/aQ  .  (2.11) 


We  shall  continue  the  analysis  as  though  asb,  but  recall  that  when  b  s  a, 
the  description  displayed  will  be  valid  provided  that,  in  each  recipe,  one  Inter¬ 
changes  a  and  b,  and  also  interchanges  the  subscripts  F  and  0.  It  is  reiterated 
that  the  value  of  Q/cp  (which  is  a  factor  that  enters  in  b,  a,  and  other  quan¬ 
tities)  is  permitted  to  vary  with  stoichiometry  and  dilution  because  in  practice 
the  burnt-gas  composition  changes;  standard  computerized  equilibrium  calcula¬ 
tions  (e.g.,  Gordon  and  McBride  1976)  permit  the  product  species,  and  hence 
Q/cp,  to  be  Identified. 

Also  In  the  reaction  zone  and  downstream  thereof,  we  have,  from  (2.5), 
(2.10),  and  (2.11), 


e(t)  -7 — 
uref 


(2.12) 
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It  follows  that,  upon  multiplication  by  9 ^  and  integration, 


where  the  Integral  is  designated  H(f)  and  the  temperature  distribution  is 
Implicit  in 

df  -  d 6/[z  *(1  +  a)  -i-  H(*)]1/2. 

uref 


']  dT’ 
(2.13) 


(2.14) 


The  details  of  #(()  are  of  no  real  Interest,  especially  in  the  absence  of  a 
chemically  more  accurate  model  for  the  reaction  rate  RR,  but,  In  particular,  the 
downstream  behavior  Is  given  by  *(«)  -  (1  +  a)  and  u(«)  ■  (1  +  a)  u(-»). 

We  do  need,  from  (2.13),  the  value  of  9 ^  at  the  upstream  end  of  the 
reaction  zone.  When  a  *  b,  ^(0-)  is  given  by 


*  * 
92(0‘)  -  2  a{l  +  a)  -i-  H(l)  -  2  0(1  +  a)  -i— 

“ref  “ref 


G(a,b  •  a , E , , ov ) 


In  these  equations  the  use  of  the  lower  limit,  T  ■  1,  in  the  integral  H  can  be 
justified  mathematically  using  conventional  singular-perturbation  techniques. 
Physically,  the  justification  stems  from  the  fact  that,  in  (2.13),  changes  in  the 
value  of  H(f)  associated  with  changes  of  9  In  the  range  1  <  9  <  [1  +  a  -  (A/E)], 
where  A  -  0(1),  are  smaller  than  the  Inaccuracies  inherent  In  the  asymptotic  pro¬ 
cedures  Invoked  In  these  flame -propagation  problems. 

Even  at  this  stage  of  the  analysis,  one  can  see,  by  comparing  (2.15)  and 
(2.16),  that  If  (for  example)  a  deflagration  wave  advances  through  a  fuel-lean 
mixture  of  hydrogen  ( a ^  *  3.5)  and  air  (*7q  •  0.85)  with  a  stoichiometry  given  by 
aj  ■  0.8  bj,  the  value  of  #^(0-)  will  differ  significantly  from  that  of  a 
second  fuel-rich  experiment  In  which  a2  -  bj,  and  b2  •  aj;  here  the  subscripts  on 
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i  and  b  merely  distinguish  the  values  of  a  and  b  in  the  first  experiment  from 
those  characterizing  the  second  experiment.  This  asymmetry  is  clearly  apparent 
in  the  laminar-flame-speed  data,  given  as  a  function  of  stoichiometric  ratio  <t>  by 
Glassman  (1987,  p.  137,  Fig.  14,  and  pp.  461-464),  who  presents  results  not  only 
for  hydrogen/air  mixtures,  but  also  for  a  variety  of  simple  mixtures  of  other 
gaseous  fuel  species  with  air.  The  data  are  mainly  from  Gibbs  and  Calcote 
(1959).  Zabetakis  (1965)  presents  a  limited  amount  of  data  for  fuel/oxygen 
mixtures.  He  shall  compare  our  theoretical  results  with  this  data  in  the  next 
section. 

In  the  region  where  $  Is  so  small  that  reactions  can  be  ignored  [most  of  the 
region  in  which  i  <  (1  +  a)],  (2.5)  implies  that 

#«  '  p  U  m  0  *  (2’17) 
Accordingly,  in  that  region,  if  A  is  a  constant  of  integration, 

l  -  1  +  A  exp(p  u  0  .  (2.18) 

Whenever  E  is  sufficiently  large  that  e as  defined  by  (2.13),  gets  rea¬ 
sonably  close  to  its  asymptotic  value  (that  at  $  ■  1)  while  e  is  still  sensibly 
close  to  the  value  (1  +  a),  one  can  match  e ^  of  the  "Inner  approximation", 
given  by  (2.13),  to  the  e ^  of  the  "outer  approximation",  obtainable  from  (2.18). 
Mathematically,  this  step  again  is  justified  by  an  appeal  to  the  asymptotic  for¬ 
malities  referred  to  earlier  in  this  section;  the  fact  is  that  the  foregoing 
approximation  gives  correct  qualitative  descriptions,  and,  over  a  broad  range  of 
parameter  values  of  practical  interest,  it  gives  quantitative  descriptions  of 
very  acceptable  and  Informative  accuracy.  Therefore,  using  this  matching  speci¬ 
fication,  (2.18)  becomes 


#  -  1  ♦  a  exp(p  u  ()  , 

1 

(2.19) 

#€(0  )  -  p  u  a  -  a  u(-«)  , 

(2.20) 

where  f^(O-)  Is  given  by  (2.15)  or  (2.16),  In  accord  with  the  magnitudes  of  a  and 
b.  In  (2.19)  and  (2.20),  the  factor  a  Is  replaced  by  b  when  a  >  b. 

This  Is  a  particularly  convenient  point  at  which  to  choose  the  dimensional 
parameter  uref  to  be  the  actual  flame  speed  for  equivalence  ratio  unity.  More 
formally,  uref  Is  the  velocity  which  [for  any  pair  of  reactants  In  stoichiometric 
proportion,  with  specified  dilution  by  a  chemical  species  nominally  neutral]  gives 
u(-«)  -  1.  Hence,  u*fif  for  a  given  pair  of  reactants  with  specified  dilution  Is 
given  by 
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(2.21) 


"ref 


2  J*(l,Nref)  *(1+ a+) 

- - G(a  ,0,E,a0,e»f 


) 


where  it  is  recalled  that  a+  is  the  value  of  a,  and  also  the  value  of  b,  that 
defines  the  stoichiometric  mixture  of  specified  dilution.  We  have  already 
anticipated  below  (2.4b)  that,  through  the  factors  *u  and  (Q/cp),  the  quantity 
J*,  defined  in  (2.4b),  is  a  function  of  the  equivalence  ratio  4>  and  the  dilution 
ratio  N,  both  previously  defined  in  Section  1  and  taken  as  specified  input;  thus, 
in  general,  we  henceforth  write  J*(*»N),  and,  in  particular,  for  the  stoichio¬ 
metric  conditions  upon  which  (2.21)  is  predicated,  we  write  J*(l,NreF).  While  in 
general  N  *  Nref  because  it  is  possible  to  conceive  of  various  comparisons  In 
numerical  calculations  presented  below,  we  shall  concentrate  on  cases  in  which  N 
•  Nref,  l.e.,  both  the  degree  of  oxygen  dilution  by  nitrogen  (and  the  examined 
pair  of  reactants)  are  fixed  while  the  equivalence  ratio  4>  is  altered.  It  fol¬ 
lows  that  u2(-«)  for  other  mixtures  of  the  same  reactants  (whether  diluted  or 
not)  is  given  by,  from  (2.15),  (2.16)  (2.20),  and  (2.21), 


«2<~) 


or  by 
u2(-) 


g(l+a)  J*U,N)  f/12  G(a,b  ~  a»E»gQ»gF> 
*(l+a+)  J V.Nref)  J  G(a+,0, E,a0,gF) 

g(l+b)  J*U,N)  fa+l2  G(b,a  -  b,E,gp,g0) 
*(l+a+)  0  U,Nref)  V)  G(a+,0,E,gF,g0) 


b  >  a; 


a  >  b. 


(2.22a) 


(2.22b) 


The  quantity  u(-«)  is  the  ratio  of  the  flame  speed  of  the  mixture  with  a,b 
stoichiometry  to  the  flame  speed  of  the  stoichiometric  (reference)  mixture. 
In  (2.2a)  and  (2.22b),  if  Cl  -  Q/cp, 


jVn) 

J  (l.Nref) 


«U(*,N)  0(l,Nref) 
"u^^ref^  QU,N) 


0(l,Nref) 

QU,N) 


(2.22c) 


If  one  Ignores  the  variation  of  the  thermal  diffusivity  with  stoichiometry. 

Also,  in  (2.22a)  and  (2.22b),  the  factors  a,b  involve  Q(*,N),  whereas  the  factors 
a+,b+  involve  Q(l,Rref). 

Equations  (2.22)  give  the  square  of  the  ratio  (of  the  flame  speed  at  gen¬ 
eral  equivalence  ratio  to  the  flame  speed  at  stoichiometry)  In  the  thermal -con¬ 
ductivity-distorted  (  coordinate,  defined  by  (2.4c).  To  express  the  same  square 
of  the  ratio  in  the  (nondimensionalized)  coordinate  x,  we  multiply  the  right-hand 
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side  of  (2.22a)  by  the  square  of  the  ratio  k(l  +  a)/k(l  +  a+),  and  the  right- 
hand  side  of  (2.22b)  by  the  square  of  the  ratio  k(l  +  b)/k(l  +  a+),  since  dx/dt  < 
(dVdf)(d?/dt).  Thus,  if  k($)  -  8s ,  6  being  a  given  constant,  then,  in  (2.22a), 
for  the  square  of  the  velocity  ratio  pertinent  to  the  x  coordinate,  the  ratio 

r  +  iQ*2+3i 

e(  1  +  a)/c(l  +  a+)  is  replaced  by  the  factor  [(1  +  a)/(l  +  a  )J 


Now  that  we  have  obtained  an  expression  for  u(-®)»  we  recall  that 

1(0  -  l/p(0  »  uU)/u(-), 

and  that:  downstream  of  the  reaction  2one, 

atf\  .  f  1  +  a,  a  <  b 
1  1  +  b,  a  >  b  ; 

in  the  reaction  zone, 

#( 0  Is  Implicit  in  (2.14)  ; 

and  upwind  of  the  reaction  zone, 


-  /  1  +  a  exP[*  » 
t  1  +  b  exp[£  u(— )]  , 


a  <  b 
a  >  b 


(2.23) 


(2.24a) 


(2.24b) 


(2.24c) 


3.  CALCULATIONS  OF  THE  FLAME  SPEED  AS  A  FUNCTION  OF  THE  EQUIVALENCE  RATIO 
3.1  Preliminary  Remarks 

By  use  of  the  simple  theory  presented  In  Section  2,  especially  In  (2.15), 

(2.16),  and  (2.22a-c),  results  were  calculated  for  the  dozen  fuel/air  mixtures 

{n  •  [?ou/[?ou  +  ^Nu)]  *  °’23}  and  fu*Voxygen  mixtures  (N  -  1)  listed  in  Table  1. 

In  Table  1,  the  Lewis-Semenov  numbers  for  the  cold  premixture  have  been  presented 
as  If  they  were  constants  Invariant  with  stoichiometry,  as  Is  done  In  some  simple 
developments  (e.g.,  Frankel  and  Slvashinsky  1983).  This  Is  a  very  rough  approxi¬ 
mation,  since  the  thermal  dlffuslvlty  of  the  cold  premixture  varies  with  the 
equivalence  ratio,  even  if  the  species  binary-diffusion  coefficients  are  some¬ 
times  reasonably  approximated  to  be  fixed  as  the  equivalence  ratio  changes  (Reid 
•t  al.  1977,  pp.  547-548;  Westbrook  and  Dryer  1981,  p.  32).  The  density,  ht?t 
capacity,  and  thermal  conductivity,  and  thermal  conductivity  of  the  cold  mixture 
vary  with  the  mole  fractions  of  the  constituents  (Penner  1957,  p.  250).  if  we 
limit  attention  to  0  s  4  s  2  In  a  slmple-hydrocarbon/alr  mixture,  then  the  mole 
fraction  of  fuel  In  the  cold  premixture  does  not  much  exceed  0.12  in  the  case  of 
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(for  example)  ethane,  though  the  mole  fraction  for  hydrogen  does  reach  0.44  over 
the  range  of  l  just  cited.  Since  the  molecular  weight,  heat  capacity,  and  ther¬ 
mal  conductivity  of  oxygen  and  nitrogen  are  similar  in  value,  if  the  cold  mixture 
has  modest  fuel  content  (as  for  fuel -lean  mixtures),  taking  the  thermal  diffusiv- 
ity  as  constant  at  the  value  appropriate  for  N2  (as  in  Table  1)  might  be  accept 
able  for  pure-fuel/alr  mixtures  --  such  a  procedure  was  discussed  in  (2.22c). 

But  the  approximation  fails  for  pure-fuel/oxygen  mixtures  (since,  as  4  -*  2,  the 
ethane  mole  fraction  approaches  0.36  for  ethane/oxygen  mixtures,  and  0.8  for  a 
hydrogen/oxygen  mixture);  erroneous  values  for  the  thermal  diffusivity  result  in 
erroneous  values  for  the  Lewis-Semenov  numbers,  and  hence  for  the  flame  speed. 

In  the  calculations  to  be  reported,  the  pre-exponential  temperature  depen¬ 
dence  of  the  reaction  rate,  a,  is  taken  to  be  2.0,  and  the  power-law  temperature 
dependence  of  the  thermal  conductivity  s  is  taken  to  be  0.5,  unless  otherwise 
stated.  Calculation  shows  that  for  fixed  N,  increasing  the  algebraic  value  of 
the  exponent  of  the  function  e(^)  ■  (0f)Q"^fi  [see  (2.5)  and  just  above  (2.23)] 
gives  a  slightly  greater  peak  value  for  the  normalized  flame  speed  u(-«),  and  the 
peak  value  occurs  for  slightly  richer  stoichiometry.  In  fact,  the  ratio  u(-»)  is 
augmented  in  value  for  all  4.  However,  since  the  appropriate  value  of  a  is  0(1) 
and  not  well  established,  and  6  •  0.5,  the  factor  does  not  contribute  signifi¬ 
cantly  and  Is  not  discussed  further. 

It  has  already  been  noted  below  (2.11)  that  a  standard  computer  code  (e.g., 
6ordon  and  McBride  1976)  yields  the  equilibrium  (burned-gas)  temperature,  for  a 
given  combustible  mixture  which  is  initially  at  standard  conditions  and  which 
undergoes  an  Isobarlc  process.  By  this  means,  the  dimensionless  burned-gas 
temperature  Tf  is  ascribed  a  physically  suitable  value,  despite  the  adaption  of 
the  one-step  pseudo-mechanism  (1.1).  Since  Tf  -  (1  +  a)  for  4,  s  1  and  Tf  -  (1 
+  b)  for  4  2  1,  self-consistency  of  the  simplistic  theory  with  a  detailed  treat¬ 
ment  (Gordon  and  McBride  1976)  of  equilibrium  in  the  burned  gas  Is  achieved  at 
any  stoichiometry  by  assigning  the  appropriate  value  for  the  ratio  (Q/cp).  The 
Implications  of  differing  dlffuslvitles  on  laminar-flame-propagation  speed  as  a 
function  of  stoichiometry  were  noted  by  Sen  and  Ludford  (1979).  However,  even 
their  ultimate  model  (Sen  and  Ludford  1982)  involves  only  a  single  product 
species,  which  Is  allowed  to  dissociate  (but  weakly)  to  the  original  fuel  and 
oxygen  only.  Such  a  model  seems  quantitatively  inapplicable  over  the  range  of 
Interesting  stoichiometries  to  any  real  chemical  system,  certainly  to  any  system 
used  in  practical  applications,  and  Indeed  no  comparison  with  flame-speed  data 
was  undertaken  in  that  work. 
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For  a  given  chemical  system  of  a  pure  fuel  and  oxidizer  (either  air  or  pure 
oxygen)  at  fixed  Initial  temperature  and  pressure,  results  are  presented  in  terms 
of  the  ratio  u(-«),  which  is  the  flame  speed  at  a  general  value  of  the  equiva¬ 
lence  ratio  4  divided  by  the  flame  speed  for  4  ■  1.  Since  the  frequency  factor 
8f  [see  (2.46)]  can  be  assigned  so  that  the  calculated  flame  speed  recovers  the 
value  experimentally  observed  for  stoichiometric  proportions  of  fuel  and  oxygen, 
for  any  ambient  pressure  and  temperature,  no  compromise  Is  Involved  In  dealing 
with  the  normalized  flame  speed  for  purposes  of  comparison  with  experimental 
data.  Effectively  equivalent  normalization  was  adopted  by  Westbrook  and  Dryer 
(1981). 

Generally,  for  a  given  system,  although  according  to  the  equilibrium  code 
the  peak  burned-gas  temperature  lies  near  stoichiometric  conditions  (meticu¬ 
lously,  slightly  on  the  fuel -rich  side),  the  Inferred  value  of  (Q/cp)  decreases 
monotonically  with  Increasing  4  for  fuel-air  systems.  For  fuel -oxygen  systems 
(N  ■  1),  the  decrease  of  (Q/cp)  with  Increasing  4  Is  appreciably  more  modest  on 
the  fuel -rich  side,  and  for  methanol  (Q/cp)  increases  slightly  with  4  for  4,  >  1; 
the  behavior  for  N  -  1  is  a  consequence  of  the  much  "broader"  maximum  of  the 
burned-gas  temperature  Tf  as  function  of  4,  in  comparison  with  the  more  sharply 
defined  peak  for  fuel -air  (N  -  0.23}  systems. 

From  (2.15),  (2.16),  and  (2.22a-c),  the  asymmetric  behavior  about  stoichiom¬ 
etry  of  the  flame  speed  as  a  function  of  equivalence  ratio  4  for  a  pure-fuel/air 
mixture  Is  owing  to  (1)  the  difference  In  molecular  weights  and  stoichiometric 
coefficients  of  the  fuel  and  oxygen  in  the  stoichiometric  balance  (for  which 
CO2  and  H20  are  the  products);  (2)  the  variation  in  the  equilibrium  product 
species  (l.e.,  the  effective  exothermicity  of  the  reaction)  with  equivalence 
ratio;  (3)  the  nitrogen  content  of  air,  which  distinguishes  It  from  a  pure  fuel 
(except  for  cases  In  which  pure  oxygen  Is  the  oxidizer);  and  (4)  the  difference 
in  values  of  the  (stolchlometrlcally  varying)  Lewis-Semenov  numbers  for  fuel  and 
oxygen.  Only  In  purely  academic  cases  In  which  no  distinction  between  fuel  and 
oxidizer  arises  on  the  basis  of  either  of  the  first  two  items,  and  In  which  pure 
fuel  and  oxygen  are  the  reactants,  Is  the  source  of  the  asymmetric  response  of 
the  flame  speed  as  a  function  of  equivalence  ratio  owing  entirely  to  a  difference 
between  the  values  of  the  Lewis-Semenov  number  for  fuel  and  the  Lewis-Semenov 
number  for  oxygen.  In  particular.  In  such  academic  cases,  If  either  the  Lewis- 
Semenov  number  significantly  exceeds  unity,  then  the  flame  speed  has  a  pronounced 
maximum  at  an  equivalence  ratio  on  the  fuel  (or  oxygen)  side  of  stoichiometry;  In 
fact,  if  both  Lewis-Semenov  numbers  appreciably  exceed  unity,  there  Is  an  evident 
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local  fl ame- speed  maximum  at  a  fuel -lean  equivalence  ratio  and  at  a  fuel -rich 
equivalence  ratio,  and  the  absolute  maximum  occurs  at  a  value  of  4  corresponding 
to  the  larger  of  the  two  Lewis-Semenov  numbers.  In  academic  cases  in  which  no 
asymmetries  about  stoichiometry  arise  on  the  basis  of  any  of  the  four  just-listed 
considerations,  then  increasing  the  dimensionless  activation  temperature  E 
results  in  a  more  rapid  decay  of  u(-«)  as  4  increases  or  decreases  from  unity. 

Of  course,  in  practice,  the  first  two  of  the  four  above-listed  considerations  do 
lead  to  asymmetry,  and  most  combustion  processes  occur  In  air,  rather  than  pure 
oxygen,  to  add  further  asymmetry. 

It  seems  useful  to  anticipate  that  the  magnitude  of  the  normalized  laminar 
flame  speed,  and  the  stoichiometry  of  the  mixture  for  which  that  maximum  occurs 
(often  at  slightly  fuel -rich  conditions),  are  fairly  well  recovered  for  fuel/air 
mixtures,  under  the  provision  that  one  adopts  somewhat  lower  values  for  the  dimen¬ 
sionless  activation  temperature  E  than  are  sometimes  cited.  However,  for  simple 
overall  mechanisms,  the  conventional  values  of  E  often  are  established  on  the 
basis  of  ignition-delay  measurements  (Mullins  and  Penner  1959,  pp.  197-202), 
whereas  appreciably  lower  activation  temperatures  have  been  adopted  In  the  past  by 
several  theoretical  modelers  to  recover  measured  flame-propagation-speed  data 
(Dugger  1951;  Fenn  and  Calcote  1953;  Gaydon  and  Wolfhard  1979,  pp.  118-121;  Coffee 
et  al.  1983).  Also,  typically  about  a  15%  reduction  is  appropriate  for  the  value 
of  E  used  here,  relative  to  the  more  conventionally  assigned  value,  owing  to  the 
approximation  adopted  in  (1.7)  to  resolve  the  cold-boundary  difficulty.  Explic¬ 
itly,  the  dimensionless  activation  temperature  E  here  is  ascribed  values  of  15, 

25,  and  66,  corresponding  to  activation  temperatures  of  about  4500  K,  7500,  and 
20,000  K,  respectively;  the  value  E  -  15  gives  results  closer  to  the  experimental 
data  (e.g.,  Gibbs  and  Calcote  1959).  This  value  Is  marginal  but  acceptable  for 
the  application  of  asymptotic  analysis. 

In  the  absence  of  heat-loss  mechanisms,  the  theory  continues  to  give  finite 
values  for  u(-«)  for  values  of  the  equivalence  ratio  4  well  outside  the  typically 
observed  flammability  limits.  Others  have  postulated  that  prediction  of  a  propa¬ 
gation  rate  below  some  small  value  Is  tantamount  to  a  prediction  of  extinction  In 
practice  (Penner  and  Mullins  1959,  p.  143),  but  attempting  such  rough  rules  of 
thumb  goes  beyond  our  objectives. 

3.2  Computational  Results 

For  the  propane-air  system,  Figure  1  presents  both  the  equllibrium-code-cal- 
culated  values  of  Tf  and  the  Inferred  values  of  (Q/cp)  vs.  4;  Figure  2  presents, 
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for  several  assignments  for  the  parameter  E,  the  resulting  values  of  u(-«)  vs.  <t>, 
together  with  the  observed  values.  Corresponding  graphs  for  the  hydrogen-air 
system  appear  as  Figures  3  and  4;  and  for  the  ethane-air  system,  as  Figures  5  and 
6.  In  all  these  results  the  Lewis-Semenov  numbers  have  been  held  constant  (at 
values  given  in  Table  1)  over  the  equivalence  ratios  examined.  On  Figure  2  are 
also  presented  numerical  results  obtained  by  Warnatz  (1985)  from  a  93-step  mecha¬ 
nism;  this  particular  result  Is  Included  because  it  is  obtained  from  a  scheme 
that  has  been  extended  to  alkanes  from  propane  through  octane,  so  It  seems 
Indicative  of  the  performance  of  a  more  widely  applicable  detailed-chemistry 
scheme  (not  a  scheme  limited  either  conceptually  or  practically  to  just  exceed¬ 
ingly  simple  fuels).  Even  with  93  reactions,  the  scheme  is  deemed  by  Warnatz  to 
be  inapplicable  to  mixtures  richer  than  those  included  on  Figure  2. 

Flame  speeds  calculated  with  the  Lewis-Semenov  numbers  permitted  to  vary 
with  cold-mixture  stoichiometry  are  more  consistent  with  the  data  for  propane-air 
(Figure  7)  and  hydrogen-air  (Figure  8).  For  the  propane-oxygen  system,  Figures  9 
and  10  are  the  counterparts  of  Figures  1  and  2,  which  pertain  to  the  propane-air 
system;  for  the  hydrogen -oxygen  system,  Figures  11  and  12  are  the  counterparts  of 
Figures  3  and  4,  which  pertain  to  the  hydrogen-air  system;  In  Figures  10  and  12, 
the  cold-mixture  Lewis-Semenov  numbers  are  permitted  to  vary  with  the  stoichiom¬ 
etry. 

The  results  make  a  point  apparently  not  previously  noted:  whereas  one-step 
models  can  approximately  predict  flame  speeds  for  simple-fuel/alr  flames  (If  the 
variability  of  the  Lewis-Semenov  numbers  with  stoichiometry  Is  retained,  and  if 
the  appropriate  exothermlcity  Is  assigned  by  use  of  an  equilibrium  calculation 
for  the  burned  gas),  the  adequacy  of  one-step  models  deteriorates  badly  for 
simple-fuel/oxygen  flames.  About  all  the  models  can  do  for  the  pure-oxygen  case 
is  to  predict  on  which  side  of  stoichiometry  the  peak  flame  speed  occurs.  For 
the  hydrogen/oxygen  case,  results  for  the  flame  speed  seem  improved  notably  if 
one  adopts  a  cold-mixture  diffuslvlty  for  the  fuel  species  (and  hence  for  the 
oxidizer  species)  which  Is  decremented  appreciably.  This  observation  suggests 
that  regarding  the  effective  fuel  species  to  be  an  intermediate,  appreciably 
heavier  species  containing  both  hydrogen  and  oxygen  Is  appropriate  for  estimat¬ 
ing  the  flame  speed;  l.e.,  use  of  molecular  hydrogen  as  the  fuel  species  Is  in¬ 
appropriate  for  the  relatively  hot  hydrogen/oxygen  flame,  even  if  it  crudely 
suffices  to  retain  molecular  hydrogen  as  the  fuel  species  for  not-too-rich  hydro¬ 
gen/air  flames.  Decreasing  the  diffusion  coefficients  for  the  cold-mixture 
reactants  in  a  propane/oxygen  flame  also  helps  bring  the  results  more  in  line 
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with  observations,  but  in  this  case  it  is  difficult  to  envision  how  species  with 
diffusion  coefficients  larger  than  those  of  the  original  reactants  may  arise. 
While  the  onset  of  soot  formation  may  be  deferred  to  slightly  higher  equivalence 
ratio  [because  the  soot-precursor-oxidation  rates  rise  more  quickly  with  tempera¬ 
ture  than  the  pyrolytic,  soot-precursor-formation  rates  (Glassman  1988)],  once 
past  the  onset-of-sooting  equivalence  ratio,  the  soot  yield  may  be  enhanced  in 
higher-temperature,  hydrocarbon/oxygen  flames  over  the  soot  yield  in  the  same- 
hydrocarbon/air  flame,  with  all  else  held  fixed.  Thus  the  overestimate  of  ob¬ 
served  hydrocarbon/oxidizer-flame  speed  at  quite  rich  stoichiometries  by  a  theory 
that  omits  radiative  heat  loss  is  to  be  anticipated  (Gaydon  and  Wolfhard  1979,  p. 
205;  Haynes  and  Wagner  1981,  pp.  232-234).  The  overestimation  is  likely  to  be 
pronounced  if  the  oxidizer  is  oxygen  rather  than  air.  Despite  these  caveats,  and 
despite  the  fact  that  empirically-guided  altering  of  the  orders  of  reaction  in  a 
one-step-kinetics  theory  might  be  of  some  avail,  probably  the  one-step  model 
itself  needs  generalization  for  flame-speed  estimation  for  fuel/oxygen  mixtures. 

For  completeness,  we  note  that  elaborate  numerical  calculations  with  multi¬ 
component-diffusion  and  detailed-chemical -kinetic  models  (Warnatz  1981)  can  very 
accurately  yield  the  flame  speed  as  a  function  of  the  equivalence  ratio  for  var¬ 
ious  pure-hydrocarbon/air  mixtures  ("up  to  butane")  and  for  hydrogen  mixtures 
with  various  ratios  of  nitrogen  to  oxygen.  What  is  less  clear  to  the  present 
authors  is  the  degree  to  which  the  success  is  attributable  to  "tuning"  some  of 
the  many  available  chemical -kinetic  parameters  purposely  to  fit  the  particular 
finite  set  of  experimental  data.  Further,  Egolfopoulos  et  al.  (1989)  document, 
for  methane-air  mixtures  initially  at  room  conditions,  the  highly  variable  suc¬ 
cess  of  different  formulations  that  adopt  (still-large)  subsets  of  the  full-chem- 
ical-kinetic  models.  For  example,  one  seriously  proposed,  plausible  subset 
Involving  fourteen  species  and  eighteen  reactions  performed  no  better  than  the 
one-step  scheme  proposed  here,  yet  the  solution  for  a  single  stoichiometry  for 
this  type  of  scheme  Is  reported  to  require  about  130-145  minutes  CPU  time  on  a 
VAX  11/785  computer. 


4.  THE  TRANSIENT  ONE-DIMENSIONAL  PROBLEM 

It  is  useful.  In  anticipation  of  a  companion  study  of  spherical -flame  propa¬ 
gation  (Carrier  et  al.  1990),  to  conclude  with  remarks  on  the  one-dimensional 
counterpart:  the  flame  propagation  that  would  ensue  when,  in  a  gas  mixture  at 
rest  in  -•  <  y  <  •,  a  flame  is  initiated  in  a  small  region  -y0  <  y  <  y0,  so  that 
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It  subsequently  evolves  as  two  flames  propagating  symmetrically  away  from  y  »  0. 
One  such  configuration  (again  with  a  <  b)  can  be  constructed  as  follows. 

Define 

x  »  y  +  (1  +  a)  u(-«)  t  (4.1a) 

and  write  for  y  <  0,  In  view  of  (2.23), 

“transient  *  -<>  +  *>  +  "<x> 

-  («(*)  -  (1  +  a)]  u(— )  .  (4.1b) 

The  coordinates  (  and  x  are  related  through  (2.4c),  and  *({)  Is  displayed  in 
(2.24). 

In  y  >  0  the  configuration  is  the  mirror  Image  of  that  In  y  <  0;  more  pre¬ 
cisely,  u  Is  odd  In  y,  whereas  p  and  8  are  even  in  y. 

The  matchup  at  y  ■  0  of  the  state  variables,  Including  velocity,  Is  fully 
consistent  with  first  principles  and  each  flow  is  merely  that  derived  in  Section 
2,  except  for  the  fact  that,  here,  each  Is  translating  relative  to  the  laboratory 
coordinates  at  a  constant  speed  of  magnitude  (1  +  a)  u(-»).  The  only  feature  of 
this  composite  flow  that  Is  Inconsistent  with  initiation  at  t  ■  0  Is  the  nonzero 
velocity  at  large  y.  Mass  conservation  requires,  of  course,  that  there  be  a  non 
zero  flow  velocity  on  the  cold  side  of  each  reaction  zone,  but  It  should  not  ex- 
tend  outside  of  |y|  ■  Y(t),  where  Y(t)  Is  the  speed  of  the  very  weak  shock  (an 
acoustic  wave,  for  all  practical  purposes)  that  is  consistent  with  the  flow  speed 

behind  the  wave  front.  Thus,  we  expect  the  flow  described  above  to  be  valid  in 

*  * 

|y|  <  Y(t)  <  aQ  t,  where  aQ  -  a0/uref’  and  ao  (dimen$ional)  speed  of  sound 

in  the  ambient  premixture  and  the  reference  speed  uref  has  been  chosen  to  be 
the  stoichiometric-mixture  adiabatic  flame  speed,  typically  less  than  4  m/s. 
Outside  of  these  "precursors",  i.e.,  in  |y|  >  a0  t,  the  motionless  state  of  the 
gas  differs  from  that  of  the  cold  gas  in  the  foregoing  analysis.  The  distinction 
is  completely  negligible,  even  for  vigorous  flames.  For  example,  the  pressure 
change  across  the  acoustic-wave  front  is  given  by 

f  .  T(t  -  l)/»0  ■  0(10'*)  .  (4.2) 

No  Independent  treatment  of  the  time-dependent  equations  is  needed  to  justify 
the  foregoing  construction,  but,  for  comparisons  with  spherical -flame  propagation, 
we  note  that  the  appropriate  conservation  laws  are  [with  k(*)  -  8~l  DF(0  - 

1  1  ®q(0  •  1,  a  •  2] 
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p  *t  +  p  “transient  *y  *  *yy  "  RR’  RR  "  RR 
p  Yt  +  p  “transient  Yy  '  *0  Yyy  "  'RR  ’ 


X,Y,$,E,- 


ref 


(4.3a) 

(4.3b) 


p  xt  +  p  “transient  Xy  '  °F  Xyy  *  *RR  *  (4'3c> 

Further,  with  s  »  y  ♦  r (t)  and  8( y,t)  -*  f(s,f),  Y(y,t)  -  Y(s,f)»  and 
X(y,t)  -*  X(s,r),  the  equations  become 

f  #  #f  ♦  #  f  #,  ♦  #  “transient  *s  *  5ss  “  RR  •  <4'4a> 

?  p  Yr  +  p  Ms  +  p  “transient  Ys  *  Yss  “  ‘RR  •  <4-4b> 

f  p  Xf  +  p  f  X$  +  p  “transient  Xs  ’  Xss  "  ‘RR  •  <4‘4c> 


We  omit  the  supporting  manipulations,  but  we  note  that,  for  the  transient 
phenomenon  described  in  the  foregoing,  6^  ■  Yj.  ■  Xf  «  0.  It  is  also  useful  to 
note  that  this  transient  wave  field  is  one  for  which  no  enthalpy  except  that  of 
the  ambient  gas  was  involved  in  the  ignition  process.  We  explicitly  point  out 
that  “transient  *  0  *n  the  *>“rned  gas. 
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Table  1.  Data  Relating  to  the  Flame  Speeds  of  Fuel -Air  Mixtures 


Fuel 

Species 

mF 

(g/mole) 

"F 

(moles) 

y0 

(moles) 

mm 

aF 

wm 

■CM! 

hydrogen,  Hg 

2 

2 

1 

0.777 

3.506 

1.80 

325 

methane,  CH4 

16 

1 

2 

0.23 

1.004 

1.08 

44.8 

acetylene,  CgHg 

26 

2 

5 

0.17 

0.767 

1.25 

155 

ethylene,  CgH4 

28 

1 

3 

0.16 

0.722 

1.13 

73.5 

ethane,  CgHg 

30 

2 

7 

0.156 

0.703 

1.14 

47.6 

methanol,  CH^OH 

32 

2 

3 

0.133 

0.600 

1.08 

50.4 

butane,  C4H1Q 

58 

2 

13 

0.096 

0.433 

1.03 

44.9 

cis-2-butene,  C4Hg 

56 

1 

6 

0.095 

0.429 

no  data 

no  data 

propane,  C3Hg 

44 

1 

5 

0.088 

0.397 

1.06 

46.4 

acetone,  C3Hfi0 

58 

1 

4 

0.082 

0.373 

.  0.93 

44.4 

benzene,  CgHg 

78 

2 

15 

0.075 

0.338 

1.00 

47.6 

heptane,  C7Hjg 

100 

1 

11 

0.075 

0.338 

1.05 

42.8 

Notes.  The  maximum  adiabatic  flame  speed  um(-«),  and  the  equivalence  ratio  4>  at 
which  it  occurs,  are  for  a  gaseous  mixture  initially  at  298  K  and  one  atmosphere 
(Gibbs  and  Calcote  1959).  Acetone,  benzene,  and  methanol  have  normal  boiling  points 
Just  slightly  above  room  temperature.  For  the  evaluation  of  Lewls-Semenov  number  for 
fuel,  the  ambient  thermal  diffuslvity  *u  »  0.222  c m*/s;  the  Lewis-Semenov  number 
for  oxygen,  oq,  Is  taken  to  be  0.83.  The  stoichiometric  coefficients  for  fuel  and 
oxygen  (*p  and  *0,  respectively)  are  for  complete  oxidation  to  COg  and  HgO.  The  fuel 
molecular  weight  is  m^.  The  ambient  diffusion  coefficients  for  the  fuel  species  are 
mainly  from  Eckert  and  Drake  (1972)  and  Hirschfelder  et  al.  (1954).  The  diffusion 
coefficient  for  propane  was  estimated  and,  In  retrospect,  the  value  may  be  about  13% 
too  small. 
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Figure  1.  For  a  propane-air  mixture  Initially  at  room  conditions,  the  equilibrium 
temperature  Tf  achieved  after  Isobarlc  burning,  and  the  consistent 
effective  chemical  exothermicity  (Q/Cp) ,  vs.  the  equivalence  ratio  <f>. 
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T,  <K) 


Figure  2.  The  laminar  adiabatic  flame  speed  u(-®),  normalized  so  Its  value  Is 
unity  for  a  stoichiometric  mixture  (*  ■  1)>  vs.  the  equivalence  ratio 
4,  for  a  propane-air  mixture  initially  at  room  conditions.  The  activa¬ 
tion  temperature,  nondimenslonalized  against  the  ambient  temperature 
Tu,  is  denoted  E.  Results  for  E  ■  15  are  closest  to  the  experimental 
observations  (Gibbs  and  Calcote  1959;  Yu  et  al.  1986;  Glassman  1987). 
Without  Inclusion  of  heat  losses,  the  simple  model  predicts  no  flamma¬ 
bility  limits.  Some  of  the  parametric  assignments  for  the  theoretical 
model  are  given  In  Table  1.  Also  plotted  are  results  calculated  with 
an  "abbreviated"  detailed-chemistry  model  (Warnatz  1985). 


1-22 


Figure  4.  Same  as  Figure  2,  but  for  a  hydrogen-air  mixture.  The  experimental  ob 
servatlons  are  from  Gibbs  and  Calcote  (1959)  and  from  Lewis  and  von 
Elbe  (1951). 


OtepilO  IgmK) 
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Figure  5.  Same  as  Figure  1,  but  for  an  ethane-air  mixture. 
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.  Same  as  Figure  2,  but  for  an  ethane-air  mixture.  The  experimental 
observations  are  reported  by  Gibbs  and  Calcote  (1959). 


1-26 


Figure  7.  Same  as  Figure  2,  except  that  the  thermal  diffusivity,  and  hence  the 
Lewis-Semenov  numbers,  are  permitted  to  vary  with  the  stoichiometry  of 
the  cold  propane/air  mixture. 
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Figure  8.  Same  as  Figure  4,  except  that  the  thermal  diffuslvlty,  and  hence  the 
Lewls-Semenov  numbers,  are  permitted  to  vary  with  the  stoichiometry  of 
the  cold  hydrogen/air  mixture. 
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Figure  9.  Same  as  Figure  1,  except  for  a  propane-oxygen  mixture. 
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Figure  10.  The  laminar  adiabatic  flame  speed  u(-«),  normalized  so  its  value 
is  unity  for  a  stoichiometric  mixture  *  1),  vs.  the  equivalence 
ratio  4,  for  a  propane-oxygen  mixture  initially  at  room  conditions. 
Results  for  E  -  15  are  closest  to  experimental  observations  (Zabetakis 
1965}.  The  thermal  diffusivity,  and  hence  the  Lewis-Semenov  numbers, 
are  permitted  to  vary  with  the  stoichiometry  of  the  cold  mixture. 
Whereas  DFjj  i  0.088  cm’/s  and  DQu  *  0.18  cma/s,  it  is  noted  that 
arbitrarily  assigning  the  smaller  values  Dfu  a  0.03  cm’/s  and 
D0u  4  0,06  cm  /*  does  better  recover  the  data;  however,  the  need 
for  such  difficult-to-justify  values  suggests  that  the  one-step 
kinetic  mechanism  needs  generalization. 
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Figure  12.  Same  as  Figure  10,  except  for  a  hydrogen -oxygen  mixture.  Again, 

decrementing  the  value  of  the  fuel  diffusivity  (Dpu  *  0.61  cm  /s  is 
the  observed  value)  better  recovers  the  data  (Lewis  and  von  Elbe 
1951).  For  hydrogen-oxygen  mixtures,  such  decrementing  may  suggest 
a  crucial  transport  role  in  flame  spread  for  a  hydrogen-containing 
radical  appreciably  heavier  than  the  hydrogen  molecule. 


1-32 


APPENDIX  II. 


Evaluating  a  Simple  Model  for  Laminar-Flame-Propagation  Rates. 

II.  Spherical  Geometry 

G.  Carrier,  F.  Fendell,  K.  Chen,  and  S.  Cook 

TRW  Space  and  Technology  Group 
Redondo  Beach,  CA  90278 

(Submitted  to  Combustion  Science  and  Technology) 


The  propagation  of  a  spherical  flame  radially  outward  (from  an  unspecified 
ignition  with  purely  local  and  transient  consequences)  is  examined  by  approximate 
analysis  for  an  initially  homogeneous  gaseous  premixture.  A  direct  one-step 
irreversible  bimolecular  second-order  chemical  reaction  with  large  Arrhenius 
activation  energy  is  adopted,  but  account  is  taken  of  the  modified  exothermicity 
owing  to  the  partial  dissociation  of  product  species  and  to  other  causes  of 
Incomplete  oxidation.  The  effects  arising  from  differing  diffusivities  for  heat 
transfer  and  reactant-species  mass  transfer,  and  from  varying  the  equivalence 
ratio  of  the  premixture,  are  considered.  Algebraic  expressions  and  a  simple 
quadrature  are  derived  which  suffice  to  describe  the  evolution  of  the  spherical 
flame  speed  to  the  asymptotic  planar-flame  value,  as  the  magnitude  of  the  flame 
radius  Increases  from  values  modestly  In  excess  of  the  diffusive  scale  to  values 
which  are  large  multiples  of  the  diffusive  scale,  for  a  flame  with  a  two-zone 
(convective-diffusive,  diffusive-reactive)  structure.  Limited  published  data  on 
the  variation  of  flamefront-speed-versus-time  behavior  with  equivalence  ratio  for 
simple-hydrocarbon/air  mixtures  are  examined  in  view  of  the  results. 


1.0  INTRODUCTION 

The  propagation  of  flames  through  combustible  mixtures  depends  in  sig¬ 
nificant  ways  on  several  characterizing  properties  of  the  constituents  and  on 
the  geometry  of  the  configuration.  This  paper  stems  from  an  interest  In  the 
multlparameter  dependence  of  the  flame  speed  of  a  spherical  flame  on  its  radius, 
on  the  transport  properties  of  the  reactants,  and  on  the  stoichiometry  of  the 
combustible  mixture  (Strehlow  1984).  Because  of  this  multi parameter  dependence 
It  Is  very  advantageous  to  so  construct  the  analysis  that  ease  of  Interpretation 
of  the  results  is  optimized.  To  that  end,  we  have  presented  an  analysis 
(Carrier  et  al.  1990)  of  a  planar  (one-dimensional)  flame,  to  which  the  spheri¬ 
cal  flame  evolves  as  It  propagates  to  larger  and  larger  radius  --  we  exclude 
from  our  examination  near-limit  flames  that  do  not  continue  to  propagate.  The 
approximations  adopted  In  these  analyses  are  designed  to  clarify  the  roles  of 
the  several  parameters  of  Interest  and  to  facilitate  comparisons  of  related 
configurations.  In  particular,  however,  we  provide  an  accurate  and  analytical 
treatment  of  these  phenomena  for  flames  of  small  radius  as  is  consistent  with 
the  basic,  well-established  but  approximate  problem  formulation  that  we  adopt. 

As  Is  the  case  with  most  combustion  problems,  It  Is  very  advantageous  to 
conduct  the  analysis  In  terms  of  dimensionless  variables  and  parameters. 
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Our  choices  and  notation  were  introduced  in  the  companion,  planar  study  (Carrier 
et  al.  1990),  and,  for  brevity,  are  not  repeated  here. 


2.0  THE  SPHERICAL  FLAME 
2.1  Formulation  and  Solution 

Me  want  to  find  a  description  of  the  combustion  phenomenon  that  ensues 
when,  at  time  t^,  a  flame  is  ignited  at  r  -  Rm^n  in  an  initially  homogeneous, 
combustible  mixture.  [The  spherical  radial  coordinate  r*  is  nondimensionalized 
against  a  diffusive  scale,  so  r  -  r*/it  l  ■  *u/uref»  "u  *s  thermal  diffusiv- 
ity  of  the  cold  mixture,  and  uref  Is  the  planar-flame  propagation  speed  for  a 
stoichiometric  mixture.]  We  want  the  description  to  be  valid  in  r  >  R^, 
t  >  t.j,  and  we  require  that,  in  r  >  R„^n,  there  be  no  artifacts  of  the  ignition 
process.  That  is,  In  t  >  tj,  the  evolution  of  the  phenomenon  stems  only  from  the 
initially  present,  homogeneously  distributed  supply  of  heat  and  reactants. 

Consistent  with  the  foregoing,  there  is  no  propagation  phenomenon  in  which 
the  flame  radius  R(t)  has  a  constant  speed  [R(t)  -  const.]  and/or  has  an  invariant 
structure  *(r,t)  -  e(x),  where  (in  this  section  only)  x  ■  r  -  R(t).  However,  we 
do  anticipate  that  the  temperature  will  be  most  conveniently  described  by 

*(r,t)  -  e[x,R(t)],  (2.1) 

and  that  its  evaluation  will  be  quasisteady.  That  is,  we  expect  that 

et  i  -  R  eXf  (2.2) 

a  statement  which  implies  that  we  expect 

e^x.R)  «  |ex(x,R)|.  (2.3) 

We  also  expect,  of  course,  that  the  thickness  of  the  reaction  zone,  within 
which  diffusion  and  reaction  are  the  significant  mechanisms,  is  small  compared 
with  either  the  radius  R  or  the  thickness  of  the  preheat  region  within  which 
diffusion  and  convection  are  the  dominant  mechanisms. 

The  conservation  laws  are  the  time-dependent,  three-dimensional,  spherically 
symmetric  counterparts  of  those  used  in  the  planar-flame  analysis  [RR  is  defined 
below  in  (2.11),  and  we  take  k(0)  -  1,  i.e.,  the  thermal  conductivity  is  held 
constant] : 


1 1  -  2 


L(0  -  p  *t  +  P  u  tr  -  *rjr  -  8°'Z  RR,  (2.4) 

P  Yt  +  P  u  Yr  -  ^[r2  Yr]r  -  8a'Z  RR,  (2.5) 

p  xt  +  p  u  xr  -  ^(f2  xr)r  -  «a‘2  RR;  (2.6) 

continuity  requires  that 

r2  pt  +  (r2  u  P)r  -  0.  (2.7) 

It  is  advantageous,  both  in  the  development  of  the  analysis  and  in  the  veri¬ 
fication  of  the  quasisteady  hypothesis,  to  treat  first  those  mixtures  for  which 
of  ■  <7q  »  1.  For  such  mixtures, 

L(f  +  Y)  -  L(*  +  X)  -  0.  (2.8) 

Under  the  initial  conditions  *(r,0)  »  *u  -  1,  X(r,0)  »  Xu  -  b,  and  Y(r,0)  »  Yu  * 
a(<  b,  for  the  fuel-lean  case),  equation  (2.8)  implies 

tf(r,t)  +  Y(r,t)  -  1  +  a,  (2.9) 

#(r,t)  +  X(r,t)  -  1  +  b,  (2.10) 


and  (2.4)  becomes  [J*  is  a  known  factor  that  varies  with  the  cold-mixture 

properties  --  see  (2.46)  in  Carrier  et  al.  (1990)] 

* 

L(«)  -  ea'z  (1  +  a  -  0(1  +  b  -  e)  exp[-E/(*  -  1)]  -  ea'z  RR.  (2.11) 
uref 

Using  (2.1),  and  also  using  the  hypothesis  that  the  length  and  time  scales  of  e 
jn  ih£  reaction  zone  are  such  that  only  diffusion  and  reaction  are  important 
mechanisms,  we  have,  as  in  the  planar- flame  counterpart, 

exx  -  -  e°‘2  RR  .  (2.12) 

In  short,  at  the  outer  edge  of  the  reaction  zone, 

e(0,R)  -  1  +  a,  (2.13) 

ex(Q,R)  »  -  0,  (2.14) 

where  (see  (2.15)  and  (2.5)  in  Carrier  et  al.  (1990),  and  recall  that  k(f)  •  1] 

* 

p2  -  2(1  +  a)°'2  -x—  G(a,b  -  a,E,l,l),  (2.15a) 

uref 
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S(a,b-a,E,c7Q,<7p) 


fl+a  r 

(b'a) (?_ a~T)  +  (^a'T)  exp[-  yq-l  dT  .  (2.15b) 

1  [  ff0  °F  J  L  J 


In  the  preheat  zone,  x  >  0,  the  reaction  rate  is  negligible,  and  (2.4)  and 
(2.7)  can  be  combined  to  give 


(r2  p  *)t  +  (r2  p  u  f)f  -  Jr2  *r]r  -  0. 


(2.16) 


We  have  confined  our  attention  to  isobaric  configurations  whose  Implication  for 
the  equation  of  state  In  the  present  approximations,  p  t  -  1,  assures  that  (2.16) 
reduces  to 


4  -  *r]  -  f(t), 


(2.17) 


where  the  function  of  integration  is  identified  as  follows.  We  anticipate  [see 
Carrier  et  al.  (1990),  Section  4]  that  the  burnt  gas  in  r  <  R(t)  [where  R(t)  is 
recalled  to  be  the  flame  position]  will  be  virtually  motionless,  so  by  (2.17)  we 
have 


u(r,t)  -  «r(r,t)  -  -j  0r[R(t),t] 

r 


(2-18) 


Equations  (2.4),  (2.13),  (2.14),  and  (2.18)  imply  that,  for  the  preheat  zone, 


p  R  eR  -  P  R  -  ex(x,R)  -  p  ^  ex 


•  Mr2  v), 


(2.19) 


e(0, R)  ■  1  +  a,  6x(0,R)  ■  -  0,  e(»,r)  ■  1, 


(2.20) 


At  this  point,  we  explore  the  possibility  that  (2.19)  can  be  modelled  (i.e., 
replaced)  by  a  much  simpler  equation  whose  analysis  would  be  more  elementary  and 
whose  interpretation  and  implementation  would  be  extremely  advantageous.  While 
(2.2)  and  (2.3)  furnish  some  of  the  tentative  anticipations,  it  is  easiest  just  to 
present  the  replacement  equation  and  list  its  attributes  than  it  Is  to  construct 
the  manner  in  which  we  “found"  It.  The  replacement  equation  is 

P  «*]x  ♦  rri r*  9x  ■  °-  (22i) 

and,  again,  it  Is  to  be  solved  under  the  conditions  given  in  (2.20).  A  series 
expansion  of  e(x,R)  in  powers  of  x,  using  (2.19),  is  precisely  the  same  as  that 
obtained  using  (2.21),  through  terms  of  order  x3.  Furthermore,  whichever  of 
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these  equations  is  used,  and  its  slope  6px  are  each  zero  at  x  -  0.  It  is 
clear  that  the  macroscopic  features  of  the  structure  of  e(x,R)  as  defined  by 
(2.19)  will  be  preserved  by  the  e(x,R)  associated  with  (2.21).  Accordingly,  we 
shall  adopt  (2.21),  as  a  suitable  approximation  to  (2.19),  but,  when  (2.21)  is 
solved,  we  shall  analyze  critically  the  accuracy  of  its  solution. 

The  solution  of  (2.21)  is 


f  ,  R2 

*  «  * 

Rx, 

^  exp 
Jx  r 

1  +  a_ 

dxj  -  1  +  I (x;£,a,R,R) . 


(2.22) 


constrained  by 

I(0;/J,a,R,R)  -  a.  (2.23) 

The  constraint  may  be  rewritten  as 


fii?  («p[rr/(1  +  a)]}  j“ 


RR/(1  +  a) 


SS^dz 


a. 


(2.24) 


Under  the  definition  of  the  exponential  integral  (e.g.,  Gautschi  and  Cahill  1967) 


E  (z )  -  f  dt  [n  -  0,  1,  2,...  ;  Real(z)  >  0],  (2.25) 

n  Jl  tn 

(2.24)  may  be  rewritten  as 

*[exp(q)]  E2(q)  •  1,  (2.26) 

where  a,  the  ratio  of  the  flame  radius  to  the  e-folding  thickness  of  the 
diffusive-convective  zone  of  the  one-dimensional  flame,  is  given  by 

A  -  0R/a,  (2.27) 

and  q/A,  the  ratio  of  the  spherical  flame  speed  (relative  to  the  burnt  gases)  to 
the  one-dimensional  flame  speed  (relative  to  the  burnt  gases),  is  given  by 

*  *(1  ♦  a)7a  *  (2>28) 

[These  characterizations  with  respect  to  the  one-dimensional  flame  are  evident 
from  remarks  in  Carrier  et  al.  (1990):  see  above  (2.!)  and  (2.15),  (2.19),  and 
(2.20)  in  that  text.)  The  quantity  q/A  is  plotted  as  a  function  of  A  in  Figure 
1.  This  figure  Indicates  that  (q/A)  -*  0  as  a  -*  1;  i.e.,  quasisteady  spherical 
flame  propagation  appears  to  lead  to  a  self-contradiction  for  a  <  1,  so,  for 
R  <  (*/£)»  *n  inherently  transient  ignition  process  seems  involved.  As  a  -  ■>, 
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« 

q/A  -  1  -  (2/a) ,  or  R/[(£/a)(l  +  a))  -  1  -  2/(£R/a);  for  Lewis-Semenov  number 
unity,  Frankel  and  Sivashinsky  (1983),  who  present  results  for  the  R  -  «  limit 
only,  give,  in  present  notation,  R/[(£/a)(l  +  a)]  -  1  -  2{[in(l  +  a) ]/a)/ (^R/a ) . 
The  two  expressions  agree  only  in  the  weakly  exothermic  limit  a  -  0.  For  future 
reference,  we  also  present  q/A  »  F(q)  -  q[exp (q) ] Eg (q)  vs.  q  in  Figure  2. 

For  i  -  1,  0  -  a.  Substitution  of  the  approximate  solution,  (2.22)  and 
(2.23),  to  the  boundary-value  problem  (2.19)- (2.20)  yields  residuals  [from 
(2.19)]  as  a  function  of  x,R.  If  the  absolute  value  of  the  residual  is  normal¬ 
ized  by  the  magnitude  of  the  diffusional  term  evaluated  at  x  •  0  and  R,  then,  for 
the  typical  value  a  •  6,  the  normalized  residual  is  equal  to  2.8xl0'5  at  x  -  0.01 
and  rises  to  a  maximum  of  about  0.38  at  x  -  5,  for  R  »  1.5;  is  equal  to  2 . 9x10" 5 
at  x  •  0.01  and  rises  to  a  maximum  of  0.092  at  x  -  1.5,  for  R  •  3;  and  i„  equal  to 
8x10”®  at  x  -  0.01  and  rises  to  a  maximum  of  about  0.075  at  x  -  2,  for  R  «  10. 

As  a  check  on  (2.3),  this  approximate  solution  gives  the  magnitude  of  ex  to  be  at 
least  ten  times  the  magnitude  of  ©p  for  all  values  of  x  fro*  zero  to  about  0.75, 
for  R  »  1.5;  to  about  1.75,  for  R  ■  3;  to  3.0,  for  R  »  5;  and  to  infinity,  for 
R  -  10. 

When  one  or  both  of  op  and  oq  differ(s)  from  unity,  another  difficulty 
arises.  If  one  proceeds  as  though  e,  Y,  and  X,  each  evaluated  on  the  hot  side  of 
the  reaction  zone,  took  on  the  values  (1  +  a),  0,  and  (b  -  a),  respectively,  the 
foregoing  procedures  applied  to  each  of  e(x,R),  Y(x,R),  X(x,R)  would  lead  to  a 
mathematical  system  for  which  no  solution  exists.  Physically,  the  imbalance  that 
underlies  that  obstacle  is  as  follows.  The  temperature  distribution  to  which  our 
procedures  would  lead  is  given  schematically  in  Figure  3,  but  the  Y  distribution 
that  satisfies  Its  boundary  condition  at  x  -  0  cannot  satisfy  Y(®,R)  -  a  because 
the  determination  of  R(R)  has  already  been  co-opted  in  the  analysis  leading  to 
e(x,R).  However,  there  is  a  temperature,  1  +  a  +  g(R),  that  can  be  assigned  to 
e(0,R),  such  that,  with  Y(0,R)  »  0  (the  reaction  zone  is  defined  to  be  the  region 
In  which  the  reaction  is  brought  to  completion),  Y(®,R)  can  have  the  value  a. 
Analogously, 

X(0.R)  -  b  -  a  +  h(R).  (2.29) 

Since  g  and  h  change  with  time,  it  is  clear  that  the  temperature  in  the 
burnt-gas  region  is  not  uniform.  It  is  also  apparent  that  the  slow  and  weak 
variation  of  ©  in  that  region  would  diffuse  innocuously  with  no  impact  on  the 
continuing  phenomenology  at  larger  values  of  r. 

Accordingly,  in  the  reaction  zone, 
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cpY+e«l+a+  g(R), 


(2.30) 


CQX  +  e«l  +  a  +  g(R)  +  ^[b  -  a  +  h(R)]. 

When  ffp  ■  1,  g  ■  0;  when  oq  »  1,  g  ■  -h;  when  •  oq  -  1,  g  ■  h  «  0;  when 

Cq  •  cp,  h  »  0. 

The  Integration  corresponding  to  (2.14)  and  (2.15)  Implies 

a  +  g(R),b  -  a  +  hfE,*0,*F] 


e*(0,R) 


2(1  +  a  +  g(R)]1 


»-2  J 


•  •> 

-  u>  r. 


Also,  by  (2.30)  and  (2.31), 

eZ(0,R)  -  c\  YZ(0,R)  -  o\  XZ(0,R). 

The  conservation  equations  in  the  preheat  zone  are  ^ex(0,R)  -  -  p  j 

,2 


p  R  -  P  R  +  p 


e<j)(x,R)  +  \ 

rc 


e 


<j)  _  fr2  6<j)). 


.  f 

r2  '• 


0, 


(2.31) 


(2.32) 

(2.33) 


(2.34) 


(2.35) 


where 


is  e,  Y,  or  X,  and  correspondingly  is  1,  ov,  or  a 


(j) 


0* 


(2.36) 


The  arguments  used  to  justify  (2.21)  are  unchanged,  so  that  (2.29)  and  (2.30)  are 
replaced  by  [a  ■  a  +  g(R)] 


P  ex]x +  777  P  sx]  ■  °- 


(2.37) 


-fP  vx)x  ♦  7T7  P  yJ  ’ °' 


(2.38) 


•oP  Xx)x  +  777  P  Xx)  ■  »•  <2'39> 

Each  of  e,  Y,  and  X  can  now  be  subjected  to  the  procedures  used  to  obtain  e 
for  the  of  •  oq  •  1  case.  One  gets 


1 1  -  7 


e  ■  1  + 


(2.40) 


The  boundary  conditions  at  x  -  0  then  imply 


it  • 


where  [see  (2.25)  and  Figure  2] 


(2.41) 

(2.42) 

(2.43) 

(2.44) 

(2.45) 


F(z)  -  z2  [exp(z)]  f  SiEtil  ds  (2.46a) 

Jz  s 

.  z  [exp(z))  E2(z).  (2.46b) 

For  any  premixture  characterized  by  given  values  of  the  parameters  E,  0O,  *F» 
a,  and  b,  the  calculation  of  R,  recalled  to  be  the  spherical  flame  speed  (relative 
to  the  burnt  gases),  as  a  function  of  R,  recalled  to  be  the  flame  radius,  can  pro¬ 
ceed  as  follows.  [It  may  be  useful  also  to  recall  that  this  section  has  been  de¬ 
veloped  under  the  fuel-lean  (or  stoichiometric)  scenario  a  s  b;  for  the  fuel -rich 
scenario  a  >  b,  the  symbols  a  and  b  should  be  interchanged,  and  the  symbols  oq 
and  of  should  be  Interchanged.)  For  each  of  a  sequence  of  real  positive  values 
assigned  to  the  quantity  q,  where  now  [In  correspondence  with  the  definition 
given  in  (2.28)] 


Q 


R  R 
1  +  a 


(2.47) 
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one  finds  the  value  of  the  function  F  [defined  in  (2.46)]  for  each  of  the  quan¬ 
tities  q,  (q/«70) *  and  (q/pp).  By  dividing  (2.43)  by  (2.44),  one  obtains  the 
quantity  g;  by  dividing  (2.45)  by  (2.44),  one  obtains  the  quantity  h.  From 
(2.15),  (2.21),  (2.32),  and  (2.33),  one  now  can  obtain  p' .  From  the  definition 

given  just  above  (2.37),  one  obtains  a'  (if  4  s  1,  b'  if  4  >  1).  From  (2.43),  one 
« 

obtains  R;  from  (2.47),  one  obtains  R.  Hence  the  sequence  of  values  assigned  to  q 
yields  R  as  a  function  of  R,  and  thence,  by  one  more  quadrature,  R  is  obtained  as 
a  function  of  time  (elapsed  since  tref,  when  R  •  Rfe^,  the  smallest  R  obtained 
from  the  sequence  of  values  of  q).  The  range  of  R  reasonably  examined  is  limited 
only  by  the  occurrence  of  a  singularity  incurred  in  computation,  or  by  the  viola¬ 
tion  of  a  condition  under  which  the  above  solution  Is  derived.  From  examination 
of  the  asymptotic  behavior  of  (2.43)- (2.47) •  we  find  that  R  -»  0  at  R  »  (a//9  )pr 

/  r 

for  4  <  1,  and  at  R  «  [b/p  )pQ  for  4  >  1,  as  generalizations  of  results  previously 
found  for  Pp  •  oQ  -  1.  It  is  reiterated  that  these  results  establish  the  scale 
smaller  than  which  an  explicitly  unsteady  initiation  process  is  inherently  in¬ 
volved  in  spherical  propagation  of  a  laminar  flame.  As  R  -*  «,  for  4  <  1,  g  and  h 
vanish  asymptotically,  and  R  goes  to  the  burned-gas  speed  of  the  planar  flame: 


R  -  |(1  +  a). 


_iI*F 


WaT* 


a.i. 


2|  a 


°F  •  po] 
(0/a)R 


(2.48) 


for  4  >  1,  the  same  formulas  hold  provided  a  -»  b,  Pp  -»  pq,  and  pq  -  Pp. 

2.2  Discussion  of  Experimental  Data 

Experimental  data  are  presented  by  Strehlow  (1984)  for  the  spherical -flame- 
* 

front  speed  R  as  a  function  of  spherical -flame-front  position  R,  for  nine  stoichi¬ 
ometries  distributed  over  the  range  0.78  <  4  <  1.68,  for  spark-ignited,  atmo¬ 
spheric-pressure  mixtures  of  propane  and  air;  three  of  the  values  of  the  equiv¬ 
alence  ratio  are  mislabeled,  the  proper  labeling  being  given  by  Law  (1989,  p. 
1394).  The  data  are  reproduced  as  Figure  4.  There  is  the  possibility  of  resid¬ 
ual  artifacts  of  the  ignition  event  in  the  data.  Frankel  and  Sivashinsky  (1983), 
in  work  cited  above  [just  after  (2.28)],  carried  out  an  analysis  valid  ffit  large 
B  only,  and  found  for  4  <  1,  in  the  present  terminology. 
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r-/n(l  +  a)  .  al(a)  2 

a/(i  +  a)  2(i  +  a)J  (1  +  a)U/a)R  * 


(2.49) 
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where 


1(a)  «  [*  if>tl  *  ^  d2‘  " 
J0  2 


(2.50) 


For  4  >  1,  a  -  b  and  *F  -  ^q,  In  (2.49)  and  (2.50).  According  to  (2.49)  and  (2.50), 

for  >  cFc,  where  o^c  permits  the  expression  in  square  brackets  in  (2.49)  to 

vanish,  R  approaches  its  asymptotic  value  from  above,  as  R  -*  •;  for  a F  <  *Fc’  R 

approaches  its  asymptotic  value  from  below.  For  +  >  1,  similar  statements  concern- 

ing  the  asymptotic  approach  of  R  can  be  made  for  oQ  >  aQc  and  <  oQc.  For  the 

propane-air  system,  oxygen  is  the  more  rapidly  diffusing  reactant,  and,  indeed,  In 

Strehlow's  data  R  does  approach  Its  asymptotic  value  from  above  for  sufficiently 

« 

rich  mixtures,  whereas  R  approaches  Its  asymptotic  value  from  below  for  other 
stoichiometries.  On  this  basis,  without  numerical  comparison  with  Strehlow's  data 
(which  is  reproduced  in  their  paper  with  the  typographical  errors),  Frankel  and 
Sivashinsky  (1983,  p.  134)  state  that  their  "theoretical  estimates  yield  both 
qualitative  and  quantitative  agreement  with  experiment".  Among  the  reasons  for 
reservations  are  the  following:  (1)  the  transition  In  whether  R  decreases  or 
increases  to  its  limiting  value  occurs  in  the  range  1.3  <  4>  <  1-43  in  the  data, 
whereas  It  occurs  at  4  -  1  In  the  theory;  (2)  the  theory  is  derived  on  the  basis 
that  m,  defined  in  (2.50),  is  constant,  but,  since  especially  a  and  even  a  vary 
significantly  with  4>,  necessarily  the  dimensionless  activation  temperature  E  must 
vary  with  4,  a  rather  academic  formulation;  (3)  the  treatment  in  the  theory  of  a 
as  constant  is  unsuitable  (as  an  example)  for  hydrogen/air  mixtures  but  no  re¬ 
strictions  on  the  theory  are  stated;  and  (4)  stoichiometric  conditions  are  not 
treated  by  the  theory. 

Frankel  and  Sivashinsky  (1983)  adopt  stoichiometrically  Invariant  values 
(for  the  pertinent  diffusivities)  that  they  ascribe  to  Pelce  and  Clavin  (1982), 
who  in  turn  state  that  the  values  are  from  Frlstrom  and  Westenberg  (1965);  we 
ourselves  were  unable  to  locate  the  cited  values  in  that  reference.  Our  calcula¬ 
tions  for  the  Lewis -Semenov  numbers  yield  different  values  for  these  numbers  (Fig¬ 
ure  5);  accordingly,  our  computations,  based  on  the  formulae  derived  in  Section 
2,  give  the  firefront  speed  to  be  a  monotonlcally  Increasing  function  of  radius 
for  all  flammable  propane/air  mixtures  (Figure  6).  The  range  of  flamefront  radii 
examined  In  Figure  6  is  far  less  than  that  in  Figure  4,  and  we  shall  return  to 
this  matter.  In  fact.  If  we  substitute  the  values  that  we  believe  to  be  appro¬ 
priate  for  the  transport  properties  of  atmospheric-pressure  propane/air  mixtures 
of  various  equivalence  ratios,  we  find  that  the  behavior  of  the  flamefront  speed 


11-10 


as  a  function  of  the  flamefront  radius  obtained  from  the  expressions  of  Frankel 
and  Sivashlnsky  (1983)  is  about  the  same  as  that  obtained  from  the  present  re¬ 
sults,  and  the  behavior  Is  quite  at  odds  with  that  observed  by  Palm-Leis  and 
Strehlow  (1969)  (Figures  7  and  8). 

If  indeed  It  were  pertinent  to  compare  the  results  of  a  theory  addressing  ef¬ 
fects  associated  with  the  Interplay  of  equivalence  ratio,  differing  diffusi vi ties , 
and  curvature,  then  It  may  be  useful  to  consider  other  data  taken  In  the  same 
apparatus  by  the  same  investigators  (Palm-Leis  1966;  Palm-Leis  and  Strehlow  1969, 
p.  118,  Figure  5)  for  atmospheric-pressure  methane-air  flames  (reproduced  here  as 
Figure  9).  If  the  role  of  a  stolchiometrically  deficient,  relatively  quickly  dif¬ 
fusing  reactant  is  the  explanation  for  an  approach  to  asymptotic  outward-spher- 
ical-flamefrcnt-propagation  speeds  from  higher  values  (that  otherwise  would  be 
from  smaller  values  --  since  the  preheating  demands  seem  geometrically  more  chal¬ 
lenging  at  smaller  radii),  then  one  might  expect  such  decreasing-firefront-speed 
behavior  in  methane/air  mixtures  to  arise  at  some  fuel -lean  stoichiometries. 
Although  the  precise  equivalence  ratios  tested  were  not  ascertained,  the  experi¬ 
mentalists  do  present  results  for  the  fi refront  speed  as  a  function  of  fi refront 
radius  for  five  mixtures  that  are  stated  to  span  a  wide  range  of  equivalence 
ratios  (from  a  near-lean-flammability-limit  mixture  to  a  stoichiometric  mixture  to 
a  near-rich-flammability-limit  mixture).  In  all  cases,  the  flamefront  speed 
clearly  Increases  monotonlcally  with  flamefront  radius,  over  the  entire  range  of 
radii  examined  (0.5  cm  to  7.5  cm).  For  completeness,  our  theoretical  results  for 
spherical  methane/air  flames  are  given  in  Figure  10. 

However,  Important  issues  arise  in  connection  with  the  very  attempt  to  use 
the  Strehlow  data  (Figure  4)  tc  validate  any  model  discussing  the  effect  of  dif¬ 
fering  diffusivities  on  outwardly  propagating  spherical  flames.  The  effect  of 
differing  diffusivities  on  planar-flame  propagation  has  been  examined  separately 
(see,  for  example,  Carrier  et  al.  1990).  Any  additional  effects  owing  to  signif¬ 
icant  flame  curvature  are  probably  absent  from  the  Strehlow  data  for  the  following 
reason.  The  spatial  scale  characterizing  the  thickness  of  a  laminar  flame  Is 
the  diffusive  scale  i  •  nu/uref#,  for  the  cold  hydrocarbon-air  mixture,  *u  • 

0(0.1  cma/s)  *nd  «ref  *  °(40  cm/$)»  so  2  is  a  few  hundredths  of  a  millimeter. 

Even  If  we  use  burned-gas  emotions,  the  estimated  thickness  but  doubles,  or  so, 

since  */*u  varies  as  (T/Tu)  ,  and  u/uref  varies  as  T/Tu,  so  i/iu  varies  as 
i/a 

(T/Tu)  ,  a  factor  equal  to  no  more  than  about  two  or  three  for  physically 
interesting  conditions.  The  Strehlow  data  are  for  flamefront  radii  of  one-to- 
eight  centimeters  or  so  (there  are  a  very  few  Isolated  readings  for  fuel -lean 
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stoichiometries  at  one-half-centimeter  radius).  Thus,  the  data  describe  phenom¬ 
ena  occurring  at  flamefront  radii  of  about  one-hundred  flame  thicknesses  or  more, 
and  thus  describe  effectively  planar  conditions.  What  Is  disconcerting  about  the 
data  is  that,  even  at  flamefront  radii  of  six-to-eight  centimeters,  many  of  the 
tests  have  not  obtained  asymptotic  conditions  --  it  may  be  argued  that  some  tests 
show  no  clear  sign  of  approaching  an  asymptote.  However,  we  may  assign  asymptotic 
values  for  each  test  by  fitting  each  curve  with  the  form  R  »  +  (ca/R  )  --  a 

procedure  that  appears  to  be  reasonable  in  practice,  although  we  offer  no  reason 
why  the  adopted  form  is  pertinent.  {Actually,  we  obtain  effectively  the  same 
estimate  for  the  asymptotic  value  for  R  if  we  simply  adopt  the  flrefront  speed 
at  the  largest  value  of  R  for  which  data  are  given.)  We  then  take  uu*,  the 
experimentally  measured  planar-flame  speeds  for  the  propane-air  mixtures  over  a 
wide  range  of  equivalence  ratios  (Yu  et  al.  1986),  and  calculate  appropriate 
burned-gas  values  (for  comparison  with  the  Strehlow  data): 
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where  it  is  recalled  that  subscript  u  refers  to  the  cold  mixture  and  subscript  f, 
to  the  burned  gas.  The  quantity  m^  is  the  molecular  weight  of  species  1  of  the 
species  present;  the  mass  fractions  for  the  burned  gas  are  incidental  Infor¬ 
mation  obtained  from  the  equilibrium-mixture-composition  computer  codes  (e.g., 
Gordon  and  McBride  1976)  in  the  course  of  calculating  the  burned-gas  temperature 
Tf.  The  molecular-weight  factor  was  omitted  from  the  previous  approximate  state¬ 
ment  of  the  equation  of  state  for  the  isobarlc  phenomena  of  Interest,  but  here  we 
seek  to  be  a  little  more  meticulous.  We  find  (Figure  11)  that  the  asymptotic  R 
values  ascribed  to  the  Strehlow  data  are  only  roughly  comparable  to  the  R  values 
inferred  from  planar-flame  data  --  and  the  Strehlow-data  values  particularly  un¬ 
dershoot  the  planar-flame-derived  values  for  the  richest  mixtures  studied  (l.e., 
for  the  largest  equivalence  ratios  examined).  We  might  ascribe  this  undershoot 
to  radlatlonal  heat  loss  from  the  burned  gas  In  the  spherical  flame  propagation; 
we  note  that  such  radlatlonal  loss  may  arise  in  connection  with  large-molecular- 
weight  hydrocarbons  and  other  gaseous  species,  and  need  not  involve  necessarily 
the  occurrence  of  carbonaceous  particulate.  One  might  ascribe  this  undershoot 
the  variation  In  sign  of  dR  /dR  with  +  for  propane-air  mixtures,  but  not  for 
methane-air  mixtures,  Is  attributable  to  the  greater  radlatlonal  loss  from  suffi¬ 
ciently  fuel -rich  gas  In  the  propane-air  case.  The  amount  of  radlatlonal  loss 
would  be  expected  to  Increase  with  the  volume  of  burned  gas.  However,  the 


11-12 


flamefront  speed  may  be  approaching  an  asymptotic  value  as  the  fi refront  radius 
increases,  but,  if  so,  it  seems  unlikely  that  the  burned-gas  volume  in  the  exper¬ 
iments  is  of  sufficient  size  that  the  burned  gas  is  becoming  opaque.  Thus  the 

•*  * 

reason  for  the  change  in  sign  of  dR  /dR  with  <t>  seems  unresolved.  The  asymptotic 
results  for  spherical  methane-air  flames  (that  are  analogous  to  the  large- flame- 
front-radius  results  for  propane-air  flames)  are  given  in  Figure  12,  together 
with  more  recent  data. 

The  just-mentioned  more  recent  data  (Ronney  and  Wachman  1985,  p.  115,  Figure 
4)  for  initially-atmospheric-pressure,  methane-air  mixtures  (over  the  range  of 
equivalence  ratios  from  stoichiometry  to  the  lean  limit)  present  results  for 
flamefront  radius  of  3  cm  or  greater.  Hence,  the  same  comments  presented  In  con¬ 
nection  with  the  Strehlow  data  again  are  pertinent.  Over  the  range  of  flamefront 

radius  of  3  cm  to  6  cm,  for  which  an  isobarlc  approximation  is  appropriate  (P. 

•  * 

Ronney,  private  communication) ,  R  Is  constant  for  all  equivalence  ratios  exa¬ 
mined,  except  for  those  mixtures  in  the  immediate  vicinity  of  the  lean-flammabil- 
ity  limit.  (The  inferred  values  for  the  planar-flame  speed  are  about  25%  below 
those  measured  by  other  techniques,  but  the  discrepancy  concerns  issues  not  perti¬ 
nent  to  the  present  discussion.)  In  contrast  to  the  data  of  Ronney  and  Wachman 

for  methane/air  mixtures,  the  data  of  Palm-Leis  and  Strehlow  for  methane/air  mix- 

•  * 

tures  does  not  report  R  to  be  constant  over  approximately  the  same  range  of 

•  * 

flamefront -radius  values.  The  explanation  for  nonconstant  values  of  R  in  this 

flamefront-radius  range  for  methane/air  mixtures  remains  unidentified  (as  does 

•* 

the  reason  for  the  nonconstant  values  of  R  in  this  range  for  propane/air  mix¬ 
tures). 


3.0  CONCLUDING  REMARKS 

We  suggest  that  a  need  remains  for  careful  experiments  on  isobarlc  spherical 
flame  propagation,  for  a  variety  of  simple -hydrocarbon/air  mixtures  over  a  wide 
range  of  equivalence  ratios,  with  emphasis  on  (1)  minimal,  transient,  localized 
energy  deposition  for  Ignition;  (2)  data  collection  at  small  flamefront  radii, 
for  purposes  of  elucidating  curvature  effects  in  the  presence  of  differing 
diffusl vities;  and  (3)  observation  whether  flamefront  speeds  consistent  with 
planar-flame  speeds  are  achieved  at  appropriate  flrefront  positions.  In  the 
absence  of  such  data,  especially  for  the  flamefront  speed  as  a  function  of 
flamefront  radius,  the  clarification  of  key  Issues  In  the  theoretical  Investiga¬ 
tion  of  spherical  flame  propagation  seems  to  be  stymied. 
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1.0 


Figure  1 


A  plot  of  q/A  vs.  x,  where  A[exp(q)]  Ea(q)  -  1,  and  Ea(q)  Is  an 
exponential  Integral. 
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Figure  3.  A  schematic  for  discussion  of  the  outward  propagation  of  a  fuel-lean 
spherical  flame,  where  e(x,R)  is  the  (nondimensional ized)  tempera¬ 
ture,  Y(x,R)  is  the  energetically  adjusted  mass  fraction  (with  value 
a  In  the  cold  mixture),  x  •  r  -  R(t),  and  r  *  R(t)  is  the  (nondimen¬ 
sional  ized)  radial  position  of  the  flame  front  at  time  t.  Hence, 
burned  gas  lies  in  R(t)  >  r  >  0  and  unburned  gas  in  ■>  >  r  >  R(t). 
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Figure  4.  Experimental  results  reported  for  atmospheric-pressure  propane/air 
mixtures,  for  the  spherical -flamefront  speed  R*  as  a  function  of 
spherical -flamefront  radius  R  ,  for  several  values  of  the  equiva¬ 
lence  ratio  4  (Palm-Lels  1966;  Palm-Lels  and  Strehlow  1969; 
Strehlow  1984;  Law  1989). 
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Figure  5.  The  Lewis-Semenov  number  for  the  fuel  species  op  and  for  oxygen  o0, 
as  a  function  of  the  equivalence  ratio  *.  The  solid  symbols  denote 
the  stolchlcmetrlcally  Invariant  values  adopted  by  Frankel  and 
Slvashlnsky  (1983)  for  propane/air  mixtures;  only  the  oj  of  the 
stolchlometrlcally  deficient  species  enters  their  model.  The  open 
symbols  denote  the  values  adopted  here  for  propane/air  and  methane/ - 
alr  mixtures;  variability  with  *  arises  from  the  variability  of  the 
thermal  dlffuslvity.  Mixtures  are  iu  atmospheric  pressure. 
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R*  (cm) 


Figure  6.  Variation  of  the  normalized  spherical -flamefront  speed  R*  with 

spherical -fl amefront  radius  R  ,  for  atmospheric-pressure  propane/air 
mixtures  of  various  equivalence  ratios  according  to  the  present 
th-.ory.  We  adopt  E  -  25,  DFu  -  0.099  cma/$,  D0u  -  0.22  cm2/s. 
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Figure  7.  The  spherical -flamefront  speed  R  as  a  function  of  spherical  - 
flamefront  radius  R  ,  for  fuel -lean  equivalence  ratios  4,  for 
atmospheric-pressure  propane/air  mixtures.  The  observations  of 
Palm-Lels  and  Strehlow  (1969)  are  compared  with  the  results  of  the 
present  theory,  and  with  the  results  obtained  from  the  formulae  of 
Frankel  and  Slvashlnsky  (1983)  fcy  substituting  the  same  parameter 
values  1M  M&  adopt. 
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Figure  8.  The  same  presentation  as  In  Figure  7,  but  here  for  fuel -rich 
equivalence  ratios. 
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Figure  9.  Experimental  results  reported  for  atmospheric-pressure  methane/air 
mixtures,  for  the  spherical -flamefront  speed  R  as  a  function  of 
spherical -flamefront  radius  R  ,  for  several  estimated  values  of  the 
equivalence  ratio  4  (Palm-leis  1966;  Palm-leis  and  Strehlow  1969). 


11-23 


Figure  10.  Variation  of  the  normalized  spherical -flamefront  speed  R*  with 

spherical -flame front  radius  R  ,  for  atmospheric-pressure  methane/air 
mixtures  of  various  equivalence  ratios  4,  according  to  the  present 
theory.  We  adopt  E  -  25,  Dpu  -  0.23  cma/s,  D0u  -  0.22  cma/s. 
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Figure  11.  For  atmospheric-pressure  propane/air  mixtures,  the  dimensionless 
spherical -flamefront  speed  R  holding  at  large  spherical -flamefront 
radius,  as  a  function  of  the  equivalence  ratio  4.  The  spherical  - 
flame  data  of  Palm-Lels  and  Strehlow  (1969),  and  the  planar-flame 
data  of  Yu,  Law,  and  Wu  (1986)  --  adjusted  with  the  burned-gas- 
temperature/cold-mixture-temperature  ratio  (Gordon  and  McBride 
1976),  appear  with  results  of  the  present  theory. 
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Figure  12.  Same  as  Figure  11,  but  for  methane/air  mixtures.  The  spherical  - 
flame  data  of  Ronney  and  Wachman  (1985)  [obtained  for  fuel-lean 
mixtures  only]  also  are  presented.  See  the  caption  of  Figure  9 
concerning  the  data  of  Palm-Leis  and  Strehlow. 
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Abstract 

Preliminary  theoretical  studies  are  undertaken  of  the  feasibility  of 
an  air-breathing  supersonic  combustor  based  on  a  stabilized,  conically 
configured  oblique  detonation  wave.  The  conical  wave  is  the  resultant  of 
the  interaction  of  a  train  of  spherical  detonation  waves,  each  directly 
Initiated  by  a  very  rapidly  repeatedly  pulsed  laser,  which  is  tightly 
focused  on  a  fixed  site  (taken  to  be  the  origin  of  coordinates)  in  a  steady 
uniform  supersonic  stream  of  combustible  gaseous  mixture.  Downwind  of  the 
Chapman-Jouguet  detonation  is  a  supersonic  isentropic  inviscid  flow.  This 
expanslonal  flow  Is  selfslmllar,  and  describable  entirely  in  terms  of  the 
spherical  polar  angle  0  (where  6  =  0  is  the  axis  of  symmetry  downwind  of 
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the  nonintruslve-energy-deposition  site,  and  0  -  p,  (jt/2)  >  p  >  0, 
identifies  the  locus  of  the  conical  detonation  wave).  The  combustor  is 
idealized  as  a  circular-cross-section  pipe  upwind  of  the  axial  position  at 
which  the  conical  detonation  wave  Interacts  with  the  wall.  Downwind  of 
this  axial  position,  the  reacted-gas  flow  is  no  longer  selfsimilar,  and  we 
seek  to  Identify  a  not-impractical ly-long,  small-drag-incurring, 
axisymmetric  nozzle  configuration  such  that  the  method-of-characteristics- 
computed  pressure  field  at  the  nozzle-exit  plane  Is  nearly  uniform  at  the 
ambient-atmosphere  value. 

1.  Introduction 

Most  designs  for  an  air-breathing  supersonic  combustor  (i.e.,  a 
combustor  in  which  the  flow  is  everywhere,  except  in  wall  boundary  layers, 
supersonic)  are  based  on  the  mixing-controlled  burning  of  diffusion 
flames1.  Here,  however,  we  examine  the  feasibility  for  propulsion  of  an 
alternative  design  based  on  an  oblique  detonation  wave  (ODW) .  Not  only  is 
the  entropy  Increase  for  an  oblique  wave  less  than  that  for  a  normal  wave 
in  the  same  flow,  so  more  of  the  combustion  heat  is  available  to  do  useful 
work,  but  also,  since  the  postdetonation  flow  is  supersonic  for  the  oblique 
wave,  downwind  disturbances  cannot  propagate  upwind  to  disturb  the  wave. 

While  ODW  combustors  have  been  discussed  theoretically  and 
experimentally  for  decades,  virtually  all  are  predicated  on  the 
stabilization  of  the  ODW  by  means  of  a  pointed  impervious  body  that  lies 
Intrusively  and  symmetrically  In  the  rapidly  flowing  combustible 
mixture1"6.  Use  of  an  Intrusive  body  without  a  positive  Ignition  device 
Incurs  frictional  drag,  and  may  Incur  the  large  entropy  losses  associated 
with  a  normal  shock  If  the  body  becomes  blunted  or  misaligned  with  respect 
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to  the  flow;  more  Importantly,  in  the  rarefied  flow  at  high  altitudes  for 
which  such  a  combustor  is  likely  to  be  utilized,  the  pressure  of  the 
(weakly)  shocked  gas  may  be  modest  enough  that  chemical  reaction  rates  are 
relatively  slow,  ignition-delay  times  are  significant,  and  the  conversion 
of  reactants  to  products  may  occur  undesirably  far  downwind  of  the  oblique 
shock4*7-10. 

2.  Laser-Initiated  Oblique  Detonation  Waves 
Because  of  constraints  on  combustor  length  in  practical  aircraft,  we 
examine  the  replacement  of  the  intruding  body  by  a  nonintrusive  positive 
Ignition  device.  Specifically,  we  consider  the  use  of  a  rapidly  repeatedly 
pulsed  laser,  tightly  focused  on  a  fixed  site  in  the  midst  of  the  flowing 
combustible  mixture,  to  effect  the  direct  photochemical  Initiation  of  a 
train  of  spherical  detonation  waves  (Figure  1).  The  combustible  mixture 
may  (or  may  not)  need  to  be  sensitized  to  the  incident  laser  radiation  by 
trace-level  doping  with  a  suitable  additive  species.  As  the  time  Interval 
between  successive  laser  pulses  approaches  zero,  the  interaction  of  the 
spherical  detonations  (Figure  2)  produces  virtually  the  same  stabilized 
conically-configured  detonation  wave  as  would  be  produced  by  a  tightly 
focused,  continuous  laser  source;  that  such  is  the  case  follows  from  the 
fact  that  the  same  element  of  mixture  can  be  converted  from  reactants  to 
products  but  once,  and  the  fact  that  the  scalloped  Interface  between 
reacted  and  nonreacted  mixture  (Figure  3)  Is  smoothed  as  the  Interval 
between  pulses  Is  decreased.  The  entropy  rise  across  reflected  shocks  that 
are  part  of  the  Interactlng-detonatlon  phenomena  Is  a  major  source  of 
losses  In  the  pulsed-detonatlon  configuration;  these  losses  decrease  as  the 
pulse  spacing  is  shortened. 
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Another  motivation  for  preferring  a  rapidly  repeatedly  pulsed  laser 
arises  from  the  desirability  of  depositing  Into  the  reactive  mixture  only 
the  minimal  energy  and  power  (over  the  minimal  volume)  necessary  to  induce 
the  direct  Initiation  of  detonation.  An  appreciably  larger  energy 
deposition,  possibly  Into  already-burned  mixture,  would  be  required  to 
stabilize  an  ODW  by  a  continuous  optical  source.  Furthermore,  continuous 
deposition  might  give  rise  to  a  locally  strong  (normal)  detonation  wave. 
Even  If  existing  pulsed  lasers  do  not  possess  sufficiently  rapid  repetition 
rate,  candidate  lasers  may  soon  be  developed,  particularly  If  the 
motivation  provided  by  a  proof-of-prlncipal  study  can  be  established  at 
this  time.  We  note  that  successful  demonstration  of  the  direct  photo¬ 
chemical  initiation  of  a  detonation  in  appreciably-subatmospheric-pressure 
mixtures  has  been  reported  by  Lee  and  coworkersll,  by  use  of  xenon  flash 
lamps. 

Observational ly,  spherical  detonation  waves  develop  a  cellular 
structure,  so  that  the  ODW  may  have  finite  thickness,  possibly  on  the  order 
of  a  centimeter. 12-14  However,  while  our  treatment  (below)  of  the  ODW  as  a 
surface  without  structure  may  not  be  meticulously  precise,  this  detail 
would  not  seem  to  be  of  critical  Importance  to  our  examination  of  ODW-based 
propulsion.  Nevertheless,  we  emphasize  the  desirability  of  laboratory 
experiments.  Highest  priority  Is  experimentation  to  establish  the  pulse- 
Initlated-detonatlon  limits  on  deposition  energy,  time  Interval  for 
deposition,  seeding  (if  necessary),  and  deposition  volume.  Furthermore, 
the  Interaction  of  two  spherical  detonation  waves  set  off  In  close  spatial 
proximity,  one  just  after  the  other,  seems  worthy  of  study  to  be  sure  that 
unforeseen  transitions  to  stronger  detonations  with  their  attending  larger 
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entropy  rise  do  not  occur.  Finally,  the  Interaction  of  a  spherical 
detonation  with  a  container  wall  ought  to  be  examined  to  be  certain  that  no 
unexpected  noise,  vibration,  or  structural,  consequences  arise  that  can  lead 
to  stronger  detonation.  It  seems  unlikely  that  detonation  will  propagate 
upwind  into  the  undisturbed  flowing  reactive  mixture  in  the  subsonic 
sublayer  Immediately  contiguous  to  the  combustor  walls. 

We  remark  that  it  is  possible  to  conceive  of  the  nonintrusive  local 
deposition  of  reduced  energy  and  power  by  a  rapidly  pulsed  repeatedly  laser 
such  that  a  train  of  spherical  deflagration  waves,  as  opposed  to  spherical 
detonation  waves,  are  induced  (Appendix  A).  The  lateral  rate  of  spread  of 
a  spherical  detonation  wave  tends  to  be  comparable  to  the  convective  speed 
of  a  supersonic  stream,  so  that  a  train  of  spherical  detonation  waves  can 
span  a  combustible  stream  with  flame  within  a  relatively  short  streamwise 
distance  downwind  of  the  spatially  fixed,  periodic-energy-deposition  site. 
In  contrast,  the  highly  subsonic  rate  of  lateral  spread  of  a  spherical 
deflagration  wave  is  so  miniscule  relative  to  a  supersonic  convective  speed 
that  the  combustible  stream  would  be  spanned  with  flame  only  impractical ly 
far  downwind.  The  only  recourse  is  to  deposit  deflagration-inducing  energy 
at  many  spanwise  sites,  but  then  the  possibility  of  transition  to  a  large- 
entropy-loss-lncurring  normal -detonation-wave  structure  arises.  Finally, 
while  deflagratlon-to-detonatlon  transition  Is  observed  in  many 
contexts^,  15f  the  residence-time  and  entropy-increase  constraints  In  the 
present  propulsive  context  makes  such  a  design  highly  problematic. 

3.  Outline  and  Objectives  of  the  Investigation 

In  a  uniform  channel  flow  at  a  speed  that  is  several  times  the 
Chapman-Jouguet  (CJ)  detonation  speed,  the  weak  detonation  (for  which  the 
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downstream  flow  is  uniform  and  supersonic)  implies  an  entropy  rise  across 
the  detonation  that  is  greatly  in  excess  of  the  entropy  rise  across  a  CJ 
detonation  in  the  same  flow;  the  strong  normal  detonation  (for  which  the 
downstream  flow  Is  uniform  and  subsonic)  implies  an  entropy  rise  that  is 
huge  compared  with  that  of  the  weak  detonation.12*15  [The  term  weak 
detonation  is  used  here  and  elsewhere  in  this  text  in  the  sense  defined 
(Ref.  15,  p.  29)  in  connection  with  situations  in  which  a  uniform  one¬ 
dimensional  burning  of  the  combustible  mixture  occurs.  A  weak  normal 
detonation  is  much  stronger  than  an  oblique  CJ  detonation  in  the  same 
flow.]  We  believe  that  we  can  select  a  channel  of  cross-section  A(x)  [with 
A ' (x)  >  0]  so  that  the  flow  downstream  of  the  quasi-conical-detonation 
region  has  very  low  losses  and  gives  rise  to  a  quasi-uniform  exit  flow 
whose  averaged  entropy  weakly  exceeds  that  emerging  from  the  detonation 
Itself. 

For  the  conical  detonation  produced  in  the  limit  of  very  short  pulse 
spacing,  the  flow  configuration,  upstream  of  the  cross-section  at  which  the 
detonation  front  Intersects  the  channel  wall,  will  be  selfsimilar  (Figure 
4).  [Background  (included  here  as  Appendix  B  for  the  reader's  convenience) 
Is  furnished  by  the  long-known  selfsimilar  property  of  the  expansional  flow 
behind  a  spherical  detonation  propagating  radially  outward  in  an  unbounded 
expanse  of  Initially  stagnant  mixture.15'18]  We  find  this  similarity 
solution  for  the  present  configuration,  and  then  continue  the  supersonic- 
flow  analysis  downstream  of  that  intersection.  Upon  analyzing  a  modest 
number  of  channel  shapes,  we  expect  to  have  a  rather  good  estimate  of  the 
performance  (thrust-to-drag  ratio)  to  be  expected  from  the  type  of 
configuration  under  study. 
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£. _ Tne  Postdetonation  Similarity  Flow 

Since  the  Rankine-Hugoniot  statements  (relating  conditions  upwind  and 

immediately  downwind  of  the  oblique,  conical,  structureless  detonation 

shown  in  Figure  4)  are  known,  we  proceed  to  quantify  that  portion  of  the 

shocked,  chemically  reacted,  supersonic  flow  that  may  be  treated  es 

selfsimilar.  This  portion  contains  fluid  elements  to  which  information 

concerning  the  finite  (cylindrical-)  radial  dimension  of  the  container  has 

not  reached.  For  the  steady  compressible  axisymmetric  inviscid  flow  of  an 

ideal  gas,  we  adopt  a  spherical-polar-coordinate  system  (r *,o,f)  anchored 

at  the  laser-energy-deposition  site;  since  dimensionless  quantities  are  to 

be  Introduced,  a  super  asterisk  denotes  a  dimensional  quantity  (i_n  Sections 

4  and  5  and  iji  Appendix  C  only) ;  r*  is  the  spherical  radial  coordinate,  9 

is  the  polar  angle  (where  9  =  0  is  directly  downwind),  and  <p  is  the 

azimuthal  angle  (of  which  the  description  is  independent).  We  shall  have 

occasion  (in  treating  the  nozzle  flow)  iv  Introduce  a  coanchored 

(alternative)  cylindrical -polar-coordinate  system  (a*,x*,^),  where  a*  is 

the  cylindrical  radial  coordinate  [0  <  a*  <  aw(x*),  with  the  container 
*  * 

radius  aw  Independent  of  x  upwind  of  the  circle  demarking  the  detonation- 

pipe  "Intersection"]  and  x*  is  the  axial  coordinate  (positive  downwind). 

In  other  sections  and  appendices,  no  dimensionless  quantities  are 

introduced,  and  the  super  asterisk  is  dispensed  with. 

We  let  subscript  o  refer  to  the  known  uniform  conditions  in  the  cold 

mixture,  and  subscript  1  refers  to  conditions  immediately  downwind  of  the 

detonation.  If  v*  Is  the  velocity,  p*  is  the  pressure,  p*  is  the  density, 

and  T*  is  the  temperature,  then  (Figure  4)  for  the  Ideal  gas  (if  R*  Is  the 
★ 

gas  constant  and  a0  is  the  speed  of  sound  in  the  cold  flow) 
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lie)  -  tyl  -  ^  r  +  ^1  e  -  m  v  +  w(*)  2  ; 


p(0)  «  .  pw  =  ^  .  T(0)  =  ■  «  T  ^  ; 


P0  /(p0  R  )  M*) 
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1/2 
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*  (  *  itr'l-1/4-  */■**■]  Hi  />] 

ao  =  iTPo^oJ  '  al  =  bpl^lJ  '  pl  *  “  '  Pi  =  “  •* 
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2  al  P1  u0  uo  s1n  P 

a,  »  ~  ,  M  -  -S-  f  M  -  -2 - 7 -  so  M  =  M  sin  p, 

*  Pl  an 

ao  0  0 


‘I 


(4) 


yt 

where  u0  is  the  (purely  axial )  cold-mixture  velocity  and  M  is  the  Chapman- 
Jouguet  Mach  number  (the  Mach  number  of  the  component  of  the  cold-mixture 
speed  that  is  normal  to  the  oblique  wave).  For  completeness,  the  Rankine- 

Hugonio*  relations  for  the  detonation  identify  [7  =  cp/cv,  a  =  7  q  /a0  , 

★  * 
where  q  Is  the  heat  of  reaction  per  mass  of  mixture,  cp  is  the  (constant) 

heat  capacity  at  constant  pressure,  and  cp  -  cv  =  R*] 


pj  ■  1  +  a  (7  -  1)  |l  + 


I_  =  1  +  a  (7  -  1) 

P\  7 
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(1* 

1/2' 
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1  + 

1  +  — y1 — 

1/2' 
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ail  ~  !)- 

d 

a  ii2  -  1) 

1/2 

A. 

L  (12  -  d] 

2/ 

T 

[“  27 
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(5) 

(6) 

(7) 


Since  (7)  gives  M  in  terms  of  input  quantities,  and  M0  is  given,  (4) 
identifies  the  angle  p. 

The  component  (of  the  cold  flow)  that  is  parallel  to  the  detonation 
and  is  unaltered  across  the  wave  is  identified: 
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(8) 


Uj  *  4  =  -5  cos  p  =  M0  cos  P  =  [m02  -  M2]172  , 

3  3 

0  0 

by  (4).  The  component  of  the  just-detonated  flow  that  is  perpendicular  to 
the  wave  is  at  the  local  sound  speed  (the  Chapman-Jouguet  condition): 

wj  *  *  =  -  (p •  (9) 

Thus,  the  streamlines  are  deflected  away  from  the  axis  of  symmetry  upon 
crossing  the  oblique  wave. 

In  the  selfsimilar  postdetonation  flow,  conservation  of  radial 
momentum  gives  (if  subscript  6  denotes  ordinary  derivative) 

ue  *  w.  (10) 


a  statement  that  the  azimuthal  component  of  vorticity  vanishes,  as  do  the 
radial  and  pclar-angle  components  of  vorticity.  Continuity  gives 


(p  w  sin  e)Q 


+  2  pu  «  0  so 


sin  e  ~  r~  - p 

The  Bernoulli  equation  gives 


p q  2u  +  w  cot  6  +  Wg 

_  ”  • 


w 


(11) 


Uu.w.-a^.Uj.Wj) 


(12a) 


with 

P/Px  ■  ( p/pi )7  •  (12b) 

Differentiating  (12a)  and  substituting  for  {pg/p)  in  (11)  gives 


2u  +  w  cot  9  + 

w2u/I 

f 

u,w;a1(7 

1  - 

w2/I 

u,w;a1( 

-T.Wj.wJ 

] 

(13) 


The  coupled  equations  (10)  and  (13)  may  be  integrated  numerically  from  6  = 
fi  to  6-6],  p  >  >0,  upon  assignment  of  input  values  to  the  parameters 
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Mq,  a,  and  7;  explicit  dimensional  results  are  obtained  by  specifying  p0 
* 

to  complement  T0,  implicitly  selected  in  assigning  a  value  to  M0  and  a. 

Near  8  -  p  the  equations  are  singular,  so  the  following  expansions  may  be 
used  to  "start"  the  solution: 

u(0)  =  Uj  +  a x  (p  -  6)  -  (2/3) e  ( p  -  8)3/2  +  ....  (14) 

1/2  2al  111/2 

w(0)  =  Wj  +  e  (p  -  8)u  -  [yj-T  [ui  '  ai  cot  P}\  *  (15) 

where  physical  considerations  select  the  negative  root  for  the  quadratic 
for  e2,  and  uj  and  wj  are  given  by  (8)  and  (9).  The  integration  may  be 
terminated  at  8  -  0j,  where 

tan  8  -  -  w(0)/u(0),  >  8  >  0,  (16) 

a  (physically  anticipated  and  mathematically  confirmed)  condition  that 
states  that  the  flow  in  0-j  >  8  >  0  is  entirely  axial  and  constant,  and  has 
no  finite  velocity  component  in  the  cylindrical  radial  direction.  In 
cylindrical  polars, 

v=Ux+va,U=u  cos  6  -  w  sin  6,  v  =  w  cos  8  +  u  sin  8  .  (17) 

Under  (16),  v  =  0,  U  =  u/cos  8.  But  if  pg  =  0,  (11)  and  (16)  give  u  =  - 
w g;  this,  with  (10),  implies  that  w  =  -B  sin  8,  B  const.;  hence,  u  =  B  cos  8, 
or  U  *  B.  Thus,  a  constant,  purely  axial  flow  in  the  core  is  consistent  with 
a  thermodynamically  uniform  core:  in  8]  >  6  >  0,  p{8 )  =  p{8\),  p(0)  =  p(0i). 

Further,  since  the  core  gas,  having  transversed  the  detonation  further 
upwind  of  where  the  gas  near  the  pipe  wall  crosses  the  detonation,  has  had 
an  opportunity  to  expand,  we  expect  the  density,  pressure,  and  temperature 
to  increase  monotonically  as  a  increases  from  o\  to  1  at  x  =  L.  Here  [if  o* 
Is  recalled  to  be  the  cylindrical  radial  coordinate,  with  0  <  a*  <  a^(x*), 
where  a#  is  recalled  to  be  the  wall  position] 
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(18) 


*  * 

'  L  S  =  tan ~p  *  =  1  tan  5i  ; 

rpipe  "pipe 


x  =  — — 
rpipe 


*>*> 

.  *w(x)  =  W - 


* 

rpipe 


(19) 


rp1pe  Is  radius  of  the  pipe;  as  previously  noted,  aw  =  1  for  0  <  x  <  L. 

Thus  L,  the  dimensionless  axial  length  from  the  energy-deposition  site  to 

the  detonation-wall  interaction,  Is  known  since  p  is  known.  It  is 

★ 

convenient  to  Introduce  a  translated  axial  coordinate  (x1  =  x’*/r  .  ) 

pipe' 


x‘  =  x  -  L  , 


(20) 


so  x’  =  0  is  the  axial  plane  containing  the  detonation-wall  interaction. 

We  present  some  results  for  the  similarity  flow  in  terms  of  6  (Figure 
5),  but  other  results  are  presented  (for  engineering  insight)  in  terms  of 
the  dimensional  cylindrical  radial  coordinate  a *  (Figures  6  and 
7)  where  at  x'*  =  0 

'  W  °  =  rp1pe  L  tan  9  *  rp1pe  <ta"  «/tan  ■  (2D 

This  form  Is  convenient  for  initiating  the  nozzle  calculation  at  x*  =  L*, 
l.e.,  at  x'*  =  0;  of  course,  It  is  necessary  to  specify  a  value  for  the 
pipe  radius.  These  figures  Indicate  that  p  and  6\  decrease  as  7  decreases 
(with  o,M0  fixed),  as  a  decreases  (with  7,M0  fixed),  and  as  M0  Increases 
(with  e, 7  fixed).  For  fixed  7,  M0,  and  a,  both  u  and  w  algebraically 
Increase  (but  p,  p,  and  T  decrease)  as  6  decreases  from  the  conical  locus 
of  the  detonation  (6  =  p)  to  the  conical  locus  of  constant-property-core 
edge  ( B  *  0j).  We  below  compute  nozzle  results  for  the  conditions  based  on 
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pQ  *  7  kPa,  T0  -  201.7  K,  but  such  near-tropopause  conditions  are  probably 
not  of  primary  practical  interest. 

5.  Nozzle  Design 

Downwind  of  the  similarity  flow,  we  undertake  an  optimization  problem 
of  inverse  type:  rather  than  seeking  the  implications  for  the  (axially- 
decreasing)  pressure  field  of  given  container  shape,  we  seek  to  identify  a 
boundary  configuration  that  gives  a  nozzle-exit-plane  pressure  field  that 
is  virtually  constant  (with  the  cylindrical  radial  coordinate  a)  at  the 

it 

ambient  value  p0. 

An  informative  illustrative  calculation  for  the  nozzle  wall  (Figure  8) 
is  the  “streamline  wall"--the  downwind  continuation  of  the  similarity-flow 
streamsurface  that  passes  through  the  detonation-pipe  intersection, 

r*  *  rpipe/sin  P'  6  =  P  (°r  a  =  lf  x'  =  °)*  If  r  =  r*/rpiP*  and  the 

ilf  it 

streamfunction  y  =  f*/(p0  a0  rp-jpe) ,  then,  if  subscripts  r  and  $  denote 
partial  differentiation, 

dy(r,0)  *  yr  dr  +  yj  r  d0  =  0,  p($)v(6)  =  curl  }}  ;  (21) 


dr  _  r  _  -  p  u  sin  0  _  uffll 
r  d0  ~  yf  r  -  /i  w  sin  ‘  w(0) 


(22) 


If  r  (0)  denotes  the  nozzle  configuration  just  defined  as  the  streamline 

5  W 

wall. 


rsw(5) 


» 

\t» 

sTn/i  exp 

. 

■P  ’ 

si 

n  e 

irsw^1^  sin  6 

«,<«</! 

0  <  e  <  ei 


(23) 


In  cylindrical -polar  coordinates,  the  equivalent  parametric  prescription 
for  <7w(x ' )  is 
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ffw(0)  ■  sin  6  [r$w(0)]  ,  x'(0)  =  cos  6  [r$w(0)]  -  L  .  (24) 

Since  U  >  0,  w  <  0,  0  <  6]  <  p,  the  streamline  wall  expands  out  radially, 

then  becomes  a  pipe  (of  circular  cross-section)  of  (cylindrical) 

radius  [sin  tf-,-]  [rsw(0-{)3  at  x'  =  x^  =  [cos  0i3  [rsw(0j)]  -  L.  While  the 

pressure  at  x'  >  x^  remains  uniform  in  a,  the  value  in  general  is  far 

★ 

above  the  ambient  value  p^. 

Thus,  a  reasonable  family  of  nozzle  configurations  to  investigate  (by 
numerical  integration  of  the  conservation  equations  for  steady  isentropic 
inviscid  axisymmetric  supersonic  flow  within  an  impervious  container,  with 
initial  conditions  at  x*  =  0  given  as  a  function  of  a  by  the  results  of 
Section  4)  are  those  that  follow  the  “streamline-wall"  shape  to  some 
distance  x'(>  0)  at  which  daw/dx'  equals  a  preselected  value;  thenceforth, 
i.e.,  for  all  larger  x',  <7w(x ' )  Is  taken  to  be  a  cone  with  angle  equal  to 
the  daw/dx'  value  selected.  The  suitable  nozzle  wall  terminates  at 
x'  =  Xg,  where  p(x^)  =  1  [I.e.,  p*(x^*)  =:  p*,  the  ambient  pressure]. 
"Tracking"  the  streamline  wall  Initiates  the  uniformization  of  the 
postdetonation  pressure  field;  "tracking"  the  conical  extension  aids  the 
pressure-reduction  process.  Rough  estimation  of  pertinent  values  for  x^  is 
given  in  Appendix  C. 

The  thrust  T*  is  given  by,  if  aw*(0)  denotes  aw*  evaluated  at  x1*  =  0, 


* 
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2  ir 


(x^*) 


*  *2 
+  P  u  ) 
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(25a) 
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l 


xe  *  *  do*(x'*) 

[»w(x'*)]{p  [*w(x'*),x'*3)  - 

0  w  w  dx'* 


dx” 


(25b) 


In  the  Newtonian  approximation  (that  molecules  encountering  a  surface 
Impact  and  slide  along  it  without  rebounding,  so  that  their  momentum  normal 
to  the  surface  is  reduced  to  zero)  appropriate  for  hypersonic  flow,  the 
wave  drag  D*  owing  to  the  flow  past  the  container  is  given  by 


*  *  r 

D  *  2*  pQ  u0  J 


*2  rxe*  0*(x'*)]  [da*(x'*)/dx'*]3 


w' 


1  +  [dff*(x'*)/dx'*]2 


dx” 


(26) 


For  a  conical  nozzle  of  (constant)  angle  6, 


★ 


D  *  *  Pq  u0  sin2^  ^ffw{xe*^2 


C^*(o)32)  . 


(27) 


We  seek  a  nozzle  configuration  such  that  the  ratio  T*/D*  sufficiently 
exceeds  unity  so,  even  with  accounting  for  the  (relatively  small) 
frictional  drag,  the  combustor  would  be  practical  for  propulsion. 

Actually,  (27)  suggests  that  a  roughly  conical  sheath  emanating  from  x*  =  0 
might  usefully  envelop  the  pipe-nozzle  configuration,  in  that  the  sheath 
might  flare  near  x  -  0  (so  that  sln^fl  is  large  where  a  is  small),  and 
evolve  to  a  pipe  at  larger  x*  (so  that  s1n20  Is  small  where  a*2  is  large). 

The  nozzle  flows  here  are  calculated  by  a  modified  method-of- 
characteri sties  (MoC)  procedure.  Whereas  conventionally  compatibility 
relations  for  the  two  families  of  characteristics  are  Introduced,  here  use 
Is  made  of  a  computational  mesh  of  streamlines  and  the  right-running  family 
of  characteristics.  This  formulation  permits  convenient  treatment  of 
nonuniform  values  of  the  molecular  weight  and  nonuniform  values  of  the 
ratios  of  the  specific  heats;  however,  In  the  present  approximation  this 
capability  was  not  utilized.  The  calculations  reported  here  typically 
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involved  a  mesh  of  70  streamlines  and  250-to-300  characteristics;  70  of  the 

characteristics  emanated  from  points  along  the  "starting  line"  (i.e., 

emanated  from  values  of  a*  at  x'*  =  0),  and  the  remainder  emanated  from 

* 

points  along  the  nozzle  wall  [i.e.,  emanated  from  values  of  <r,(x'*)].  The 

w 

execution  time  was  typically  15  minutes  on  a  personal  computer  with  the 
performance  of  an  IBM  AT.  Several  cases  were  checked  by  computing 
essentially  identical  results  with  the  VNAP2  computer  program,  which  adopts 
the  unsplit  MacCormack  scheme19;  however,  In  other  cases,  extensive 
“tuning"  was  required  before  VNAP2  yielded  mathematically  and/or  physically 
plausible  results,  so  eventually  th«i  MoC  procedure  was  used  exclusively. 

The  MoC  calculations  indicate  that,  for  the  conditions  at  x'*  =  0 
given  by  the  solid  lines  In  Figures  7  and  8,  a  nozzle  contour  that  adopts 
the  streamline  wall  to  a  distance  x'*  =  12.67  cm,  and  then  smoothly  is 
continued  as  a  5°  right-circular  cone,  results  In  a  pressure  of  6.94  kPa 
very  nearly  uniformly  at  all  radial  positions  0  <  a*  <  aw(x'*),  where 

W  0 

★ 

x'*  =  3.0  m  and  a  =  47.5  cm  at  this  axial  distance.  Meticulously,  the 
pressure  at  the  axis  of  symmetry  does  slightly  exceed  the  cited  pressure, 
computed  to  hold  at  the  wall.  As  an  Indication  of  the  sensitivity  of  the 
calculations,  we  present  results  (Figures  9,  10,  and  11)  for  a  nozzle 
contour  that  follows  the  streamline  wall  to  x'*  =  10  cm,  then  smoothly  is 
continued  as  a  5.93°  cone  to  a  distance  x^*  2  2.5  m;  actually,  at  axial 
distance  2.5  m,  for  which  the  nozzle-wall  radius  is  46.2  cm,  the  pressure 
at  the  wall  has  expanded  past  7  kPa  to  6.48  kPa. 

6.  Concluding  Comments:  Effect  of  a  Nonuniform  Mixture 
Since  achieving  intimate  mixing  of  fuel  injected  into  high-speed  air 
is  a  major  topic  of  current  high-speed-propulsion  research  (and  of  high- 
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energy-chemical -laser  research),  In  practice  the  cold  stream  Is  expected  to 
consist  of  imperfectly  mixed  fuel  and  air.  We  discuss  here  briefly  a 
throughput  flow  consisting  of  a  fraction  composed  of  well -mixed  fuel  and 
air,  through  which  is  Interspersed  imperfectly  mixed  blobs  with  fuel  and 
air  in  stoichiometric  proportion.  We  do  not  delve  here  into  the  creation 
of  such  a  mixture,  and  the  scale  of  the  blobs  is  to  be  determined  by 
further  research;  nevertheless,  a  millimeter  scale  may  be  plausible.  The 
well-mixed  fraction  of  the  throughput  is  converted  to  product  gases  across 
the  weak  detonation  (Figure  12). 

The  "lumpy"  blobs  are  compressed  during  passage  through  the 
detonation,  but  these  blobs  remain  unreacted  over  the  centimeter-or-so 
scale  of  the  (cellular)  oblique  wave.  They  do,  however,  react  in  the  hot 
flow  downstream  of  the  front,  as  an  isolated  diffusion-controlled  reaction 
occurs  in  each  blob.  If  the  blobs  are  of  dimension  a,  their  burn-up  time 
Is  about  a2/D,  where  D  is  the  Ficklan  diffusion  coefficient  of  the  lighter 
reactant.  In  fact,  If  one  takes  the  liberty  of  Introducing  the  concept  of 
an  equivalent  deflagration,  the  "lumpy"  blobs  have  a  relatively  slow  flame 
speed.  The  fact  that  unreacted  pockets  of  reactive  mixture  occur  behind  a 
real  detonation  front  Is  well  appreciated  in  a  solid  detonable  mixture.20 
However,  this  phenomenon  In  solid  explosives  is  a  consequence  of 
uncontrollable  departure  from  Ideal  behavior.  The  distinction  here  Is  that 
we  raise  the  possibility  of  intentionally  tailoring  a  gaseous  mixture  to 
control  the  nature  of  the  departure  from  perfect  mixing.  For  example,  If 
the  total  heat  release  per  unit  mass  of  mixture  is  denoted  r,  and  If  the 
heat  release  per  unit  mass  of  mixture  across  the  weak  detonation  Is  a,  then 
we  seek  to  extract  the  heat  (f  -  a)  without  Incurring  any  shocks,  or  at 
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most  minimal  shocks.  While  the  entropy  at  the  downwind  plane  at  which  the 
mixture  is  fully  converted  to  product  must  be  higher  than  the  entropy  just 
after  the  weak  detonation,  we  seek  to  achieve  the  lowest  possible  entropy 
state  in  the  fully-burned  gas  (given  the  upwind  thermodynamic  state  and 
flow  speed,  and  given  a  quasi-one-dimensional  flow  at  the  downwind  exit 
plane  with  cross-sectional  area  about  the  same  as  that  for  the  upwind  cold 
flow).  Thus,  research  into  whether  attempting  to  tailor  imperfect  mixing 
offers  a  possibility  for  enhancing  performance  may  be  worthwhile. 

Appendix  A.  Energy  Deposition  for  the  Initiation  of  Detonation 
and  Deflagration 

Consider  a  gaseous  mixture  with  composition  such  that  the  mixture  lies 
within  the  limits  of  detonation  and  deflagration.  In  order  to  Initiate  a 
deflagration  wave  with  a  nonlntruslve  local  deposition  of  energy  to  the 
combustible  mixture,  one  must  supply  an  appropriate  amount  of  energy  to  the 
blob  of  mixture,  irradiated  during  a  time  Interval  that  is  small  compared 
with  the  time  required  for  the  loss  of  much  of  that  heat  by  diffusion  to 
the  contiguous  medium.  These  requirements,  and  their  counterparts  for  the 
initiation  of  a  detonation,  may  be  only  loosely  quantifiable,  but  the 
effort  Is  a  good  place  to  start. 

Let  the  radius  of  the  blob  to  receive  the  deposited  energy  be  denoted 
by  b,  so  that  the  mass  of  the  (say,  roughly  spherical)  blob  Is  (4/3)*y>0b3, 
where  p0  Is  the  density  of  the  ambient  (unreacted)  premixture,  and  b  refers 
to  the  unheated  state.  While  the  minimum  energy  deposition  required  to 
Initiate  a  continuing  deflagration  cannot  be  found  trivially,  it  Is 
sufficient  to  deposit  an  amount  of  energy,  E,  equal  to  the  amount  of  heat 
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that  complete  combustion  of  that  mass  of  premixture  would  provide;  i.e.,  E 
=  (4/3)*7?0b3q,  where  q  Is  the  heat  of  reaction  per  mass  of  mixture. 
Diffusion  to  contiguous  material  would  occur  on  a  time  scale  Tq  2  b2/x, 
where  c  characterizes  the  thermal  diffusivity  of  the  ambient  medium;  thus, 
the  deposition  time  T  should  obey  T  «  b2//c. 

A  detonation  will  be  directly  initiated  if,  in  addition  to  the 
foregoing,  the  deposition  time  T  Is  so  brief  that  the  acoustic  transit  time 
across  the  blob  (or  half  of  it)  Is  large  compared  to  T.  That  Is,  If  a0  is 
the  ambient  speed  of  sound,  the  energy  deposition  may  initiate  a 
deflagration  (but  not  a  detonation)  If  (a0T/b)  »  1;  but  a  detonation  may 
be  initiated  if  the  deposition  time  T  (and  also  the  time  required  for  local 
reaction)  obeys  T  «  b/a0. 

Thus,  if  (as  an  example)  it  -  0.1  cm2/s,  b  =  0.1  cm,  and  a0  =  3 - 104 
cm/s,  for  a  deflagration  (with  T  In  seconds)  3-105  T  »  1  »  10-1  T.  Thus 
T  *  10"2  s  and  E  *  (4/3)*/>0b3q  would  Imply  Initiation  of  a  deflagration 
with  very  little  chance  of  either  failure  to  Ignite  or  initiation  of 
detonation.  For  the  same  Illustrative  parameters,  one  could  expect  with 
confidence  to  initiate  a  detonation  if  the  deposition  of  E  units  of  energy 
occurred  in  10”?  s.  Experiments  are  needed  to  obtain  a  more  accurate 
criteria  for  the  initiation  of  either  a  deflagration  or  detonation.  It  is 
easily  possible  that  the  minimum  time  for  energy  deposition  for  initiation 
of  a  detonation  could  be  only  10-6  s  and  that  the  minimum  energy  could  be 
as  small  as  (4/3)i/>ob3q/10.  The  point  Is  that  somewhere  between  deposition 
times  of  10-2  s  and  10"7  s  the  situation  gets  "delicate"  and  the  need  of 
study  Is  evident. 
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The  minimum  (normalized)  amount  of  deposited  energy  required  to 
achieve  direct  initiation  of  a  detonation  plausibly  varies  with  the 
(normalized)  time  of  deposition  as  sketched  in  Figure  13.  That  is, 
E/[(4/3)ir/>0b3q]  Increases  with  relatively  large  values  of  (a0T/b),  simply 
because  some  of  the  earlier-deposited  energy  is  carried  away  from  the  blob 
prior  to  the  achievement  of  detonation  onset  by  later-deposited  energy;  for 
smaller  values  of  (a0T/b),  no  such  dispersing  of  deposited  energy  occurs, 
and  the  minimum  energy  requirement  for  direct  initiation  suffices.  There 
is  experimental  evidence  (albeit  presented  in  the  dimensional  form  of 
critical  initiation  energy  vs.  time  of  deposition)  that  is  fully  consistent 
with  this  discussion  [Ref.  12,  p.  407,  Figure  13-7(a)]. 

For  Initiation  of  detonation,  the  heated  gas  will  be  confined  to 
approximately  Its  original  volume,  so  that  Its  early  temperature  rise  (AT) i 
*  q/cv,  where  cv  Is  the  specific  heat  at  constant  volume  of  the  gas. 
Specifically,  for  a  spherical  detonation  wave,  the  heated  gases  are 
strongly  confined  by  the  dynamics  of  the  propagation  mechanism,  so  that, 
just  behind  the  wave  front,  the  temperature  has  a  maximum  whose  order  of 
magnitude  Is  (AT)j,  and  the  pressure  near  the  wave  front  is  of  order  p-j , 
where 

Pi  =  Po  R[(AT)i  +  T0] , 

and  R  is  the  gas  constant  and  T0  is  the  ambient  temperature.  These 
estimates  are  based  on  the  above-noted  possibility  (confirmed  in  some  cases 
Involving  deflagrations)  that  the  total  energy  absorbed  from  the  laser  beam 
into  the  sample  need  be  but  a  small  fraction  of  the  heat  of  combustion.  In 
contrast,  In  a  deflagration  the  time  of  deposition  Is  long  enough  that  the 
burnt  gases  will  have  expanded  during  the  burning  to  a  new  volume  at  which 
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the  pressure  is  virtually  the  same  as  the  ambient  pressure  (in  the  unburnt 
gas).  For  a  deflagration,  the  rather  uniform  temperature  rise  in  the  burnt 
domain  is  approximately  q/cp. 

Appendix  B.  Selfsimilar  Treatment  of  the  Shocked  and  Reacted  Gas  in  a 
Spherical  Detonation 

The  state  of  the  shocked  and  reacted  gas  just  downstream  of  a 
spherical  Chapman-Jouguet  detonation  (Ref.  17,  p.  195,  Figure  43)  is 
characterized  by  pressure  pf,  density  pf,  sound  speed  a-f,  and  particle 
speed  uf,  such  that  Uf  relative  to  the  speed  of  the  detonation  wave  is 
sonic,  l.e.,  Is  equal  to  af[=  i'jPf/pf)^^,  for  the  perfect-gas  model 
adopted  here].  If  R (t)  is  the  radial  position  of  the  outwardly  propagating 
spherical  wave,  then  in  0  <  r  <  R(t),  for  time  t  >  0,  the  pertinent 
equations  are  (7  =  cp/cv) 

r2Pt  +  (Pr2 u)r  s  0  *  (B.l) 

/jut  +  />uur  +  pr  =  0  ,  (B.2) 

P/Pf  =  (/>//>f)7-  (B.3) 

•  • 

The  position  R(t)  =  Rt,  where  R  is  a  constant,  specifically,  the  speed  of  a 
Chapman-Jouguet  detonation.  That  is,  we  focus  on  constraints  that  would  be 
the  result  of  an  extrapolation  back  in  time  of  a  selfsimilar  solution16-18 
of  the  above  conservation  laws  subject  to  the  thermodynamic  properties  of 
the  shocked  and  reacted  gas.  Hence,  at  r  =  R(t), 

• 

P  *  Pit  P  *  Pft  Uf  =  R  -  af  .  (B.4) 

The  selfsimilar  variable  Is  necessarily  of  the  form  r/  -  r/(kt),  where 
k  is  a  parameter  of  the  problem  that  has  the  dimensions  of  velocity. 
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Without  loss  of  generality,  we  adopt  k  =  af.  We  define 

p(r,t)  =  pf  G(j 7),  p(r,t)  =  pf  H (??),  u(r,t)  =  af  F(»?)  .  (B.5) 

Under  (B.5),  the  statements  (B.l)  -  (B.4)  become,  since  H  =  G^,  two  coupled 


quasilinear  ordinary  differential  equations: 

*7 [ (F  -  V)2  ~  g7_1]G,  s  *2  F(F  -  »?)G*  (8-6) 

f}[(t  -  rj) 2  -  G?"1^  =  2  FG^-1  ,  (B.7) 

« 

over  0  <  17  <  X,  X  =  R/af,  where 

G(X)  =  1,  F (X)  =  X  -  1  .  (B.8) 

This  system  may  be  reduced  to  a  single  autonomous  equation  by  judicious 
substitution16  but  that  reformulation  does  not  appear  to  aid  present 
purposes  of  extracting  quantitative  profiles.  The  system  (B.6)  -  (B.8)  is 
singular  at  tj  =  X,  and  it  helps  to  deduce  that,  near  rj  =  X, 

F  a  X  -  1  ■  r(X  •  rj)112  ,  G  a  1  •  r(X  ■  rj)l/2  ,  (B.9) 

r  *  2 { (X  -  1 ) / [ (7  +  1)X]}1/2  .  (B. 10) 

With  the  aid  of  (B.9)  and  (B.10),  numerical  integration  of  (B.6)  -  (B.8) 
from  r\  -  X  toward  smaller  rj  may  proceed;  the  integration  proceeds  to  i;  = 
fj\  (<  X),  where  F(^^)  =  0  —  in  fact,  F(> 7)  =  0,  0  <  rj  <  rj\ ,  so  G {rj)  =  G(rj]) 


in  0  i  rj  S  rj].  Some  numerical  results  are  graphed  by  Sedov17,  but  these 
seem  not  to  be  easily  accessible,  so  some  results  are  presented  here 
(Figure  14).  Barenblatt^1  has  dwelled  upon  the  not-enti rely-smooth 
behavior  earlier  noted  by  Zeldovich  to  hold  at  rj  =  rjy,  for  our  purposes, 
this  behavior  is  very  familiar  In  hyperbolic  systems  and  violates  no 
criteria  with  significant  physical  Implications.18 
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Appendix  C.  A  Crude  Estimate  of  Nozzle  Lenqth 


Adopting  the  definitions  and  relations  of  (l)-(9),  we  introduce  the 

★ 

stagnation  speed  of  sound  for  the  ambient  flow,  a$0,  and  for  the  state  just 

★ 

downwind  of  the  oblique  detonation,  a$^: 


*  2  *2  *2*2*2  * 

Ho  ■  ■»  +  [<7  -  ««]  %  •  Hi  'Ho  +  (7  -  1)  fl  • 


(C.1) 


The  speed  of  the  flow  just  behind  the  detonation,  Iv^l,  is  known  from  (4) 


and  (7),  since 


♦*  o  *c  2 

IVjl  *  aj  +  uQ  cos  p  . 


(C.2) 


Since  subscript  e  denotes  conditions  holding  at  the  exit  plane,  by 
definition  or  deduction  it  follows  that 


**  2  *2  *2  *  *  *2  *2  *  * 

Ivj  =  +  \ia  ,  pa  =  p„  ,  ac/1  =  ae1  ,  pco  =  p. 

e  e  e  re  ro  se  si  rse  rsl 


(C.3) 


since  the  exit  plane  flow  is  an  Isentropic  expansion  of  the  flow  just 

* 

downwind  of  the  detonation;  p$1-  is  defined  to  be  the  stagnation  pressure 
at  state  1. 


From  Bernoulli's  equation  and  (C.3),  psg  is  obtained: 


pse  _  psl 


»  1  - 


I  v*l2l  -7/(7-D 

-  i  'V 


(C.4) 


Also  from  Bernoulli's  equation  and  (C.3), 

Po*  f.  ,.1W(,-1). 


(C.5a) 


by  rearranging,  we  obtain  lvgl: 
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<4* 

'V  . . 

2 

* 

asl 

7  -  1 

k 

Since 

*  1  pse 

p$e*  *  2  * 

1  - 


1(7-1) /7 


'seJ 


1/2 


(C.5b) 


(C.6) 


se 


from  (C.5a)  we  obtain  pQ: 


e 

* 

°$e 


i  - 

1  2 


lsl 


l/(7-D 


(C.7) 


If  the  flow  is  nearly  one-dimensional  with  cross-sectional  area  A  ,  we 
obtain  a  crude  estimate  of  the  cylindrical  radius  of  the  nozzle  at  the  exit 
plane  a*[s  a*(x^*)]: 


*2 


*  * 
PP  IvJ  A 


3.  Un 
0  0 


A  so  —  =  — 

An  / 

o  r, 


*  * 
pn 

ro  o 


pipe 


*  .** 
IvJ 


(C.8) 


Whereas  (C.8)  results  in  general  in  an  implicit  equation  for  x^*,  for 
a  conically  shaped  nozzle  (of  selected  angle  0)  joining  the  pipe  radius 


rpipe  t0  the  ex^  radl'us  °e»  the  nozzle  length  x^*  is  roughly  estimated 
to  be 


)/tan  6  . 
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Figure  1.  Each  of  the  periodic  pulses  nonlntruslvely  deposits  sufficient 
energy  [In  a  uniform  supersonic  stream  (Mach  number  M0)  of  combustible 
mixture,  flowing  In  an  axisymmetric  container]  to  initiate  directly  a  CJ 
detonation,  which  expands  radially  as  It  Is  convected  downwind.  The 
envelope  of  detonation-front  positions  Is  Indicated. 


Rl-477-90 


III -27 


— 1 "0*2  “*1 
*  *  I  I 

»o«r*i  ! 


Preflected 


knunu\\n\\\uvi\\uuV/*\\u%\u\\uuu\UHU%u%Mw* 


J _ REFLECTED 

SHOCK 


!tmwu\mw\wmww«www«™i 


FUEL-AIR  ESSMEip 

pSlURE  «G'°N  0F  INTERACTI°N 


FOCAL  POINT  OF  IRRADIATION 


Figure  2.  A  simplified  schematic  showing  three  detonation  fronts  Initiated 
in  a  premixed  supersonic  stream  which  is  travelling  faster  than  the 

it 

detonation  fronts  (which  travel  at  speed  u0  sin  p).  The  downwlndmost  front 
has  Interacted  with  the  nonflared  wall,  and  a  resulting  shock  travels 

it  it  it  it  it 

through  already-reacted  gas.  Time  t3  >  t2  >  tj,  such  that  (t3  -  t£)  = 

(*2  “  tl)* 
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Figure  3.  The  "scalloped"  front  constituted  by  portions  (actually,  zonal 
strips)  of  the  Interacting  spherical  CJ  detonation  waves  initiated  by 
pulsing  at  a  fixed  finite  frequency. 
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Figure  4.  Schematic  (not  to  scale)  of  a  supersonic  combustor  (x  =  x'  +  L 
n,  t  are  normal,  tangential  coordinates).  The  selfslmllar  solution 
actually  holds  for  x*  >  0  near  the  axis,  as  determined  by  the  "Inward¬ 
running"  characteristic  from  the  nozzle-pipe  joining. 
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6  (rad) 


Figure  5.  Dimensionless  (spherlcal-)radial  velocity  component  u  and  polar- 
angle  velocity  component  w,  as  functions  of  the  polar  angle  6,  for  various 
values  of  the  ambient  Mach  number  M0l  heat-capacity  ratio  7,  and  chemical- 
heat-to-static-enthalpy  ratio  a.  The  range  depicted  Is  0-|  <  0  <  p,  i.e., 
from  the  uniform-core-edge  cone  to  the  detonation-wave  cone. 
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Figure  6.  Axial  velocity  component  U*,  (cyllndrlcal-)radlal  velocity 

component  v*,  and  density  p* ,  computed  from  the  similarity  solution  and 

presented  as  functions  of  the  (cyllndrlcal-)radlal  coordinate  a*.  Here, 

★ 

M0  *  10,  a  *  30,  7  *  1.4,  r  =  0.2  m ,  P  -  49°.  For  the  solid  curves, 
p*  *  7  kPa,  T*  =  201  K;  for  the  dashed  curves,  p*  =  0.1  kPa,  T*  =  271  K. 
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Figure  7.  Temperature  T* 
solution  and  expressed  as 
a*,  for  the  conditions  of 


and  pressure  p*t  computed  from  the  similarity 
functions  of  the  (cyl indri cal-) radial  coordinate 
Figure  6. 


Rl-477-90 


1 11-33 


Figure  9.  The  nozzle-wall  geometry  for  the  so-called  streamline  wall 
(drawn  to  the  axial  distance  at  which  the  pressure  Is  unlformized 
radially),  and  also  for  the  wall  shape  (examined  In  Figures  10  and  11)  of 
10-cm  span  along  the  streamline  wa' 1  fc! lowed  smoothly  by  a  5.94°  cone. 
The  conditions  are  those  in  Figure  6  (solid  curves). 


R1 -477-90 


1 1 1 -35 


Figure  10.  Pressure  p*  at  the  axis  of  symmetry  a*  *  0  and  at  the  nozzle 

•k  k 

wall  a  *  aw,  as  a  function  of  the  axial  distance  x'*.  For  the  geometry, 
see  the  "cone"  In  Figure  9;  for  the  starting  conditions  at  x'*  =  0,  see  the 
solid  curves  In  Figures  6  and  7.  Results  are  from  a  method-of- 
characterl sties  (MoC)  code  and  from  the  VNAP2  code19. 
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Figure  11.  Pressure  p*  as  a  function  of  radial  coordinate  a*  at  the  exit 
plane,  for  the  case  discussed  In  Figure  10.  The  exit  plane  x^*  =  2.47  m 
for  VNAP2  results,  *  2.5  m  for  MoC  results. 
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Figure  12.  The  half-shaded  spheroids  of  scale  a  denote  the  (I-N)%  of  the 
mixture  constituted  by  blobs  consisting  of  fuel  (shaded)  and  air  (unshaded) 
In  stoichiometric  proportion.  The  well-mixed  N%  Is  converted  to  product 
gas  In  the  oblique  detonation. 
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Figure  13.  A  rough  sketch  of  how  the  minimum  deposited  energy  required  for 
the  direct  onset  of  detonation,  E,  normalized  against  the  chemical  heat  of 
reaction  available  in  the  Irradiated  blob  (taken  for  specificity,  to  be  a 
sphere  Initially  of  radius  b,  of  density  p0,  and  of  exothermicity  per  mass 
q),  Is  anticipated  to  vary  with  the  deposition  Interval  T,  normalized 
against  (b/a0),  where  a0  is  the  speed  of  sound  In  the  ambient  detonable 
gaseous  mixture. 
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Figure  14.  For  a  spherical  detonation,  the  dimensionless  flow  speed  F  of 
the  shocked  and  reacted  gas,  as  a  function  of  the  similarity  variable  rj 
defined  below  (A. 4),  for  the  ratio  of  specific  heats  7  =  1.4,  for  several 
values  of  the  ratio  X  of  the  CJ-wave  speed  to  the  just-detonated-gas  sound 
speed. 
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Abstract 


We  have  achieved  volumetric  ignition  of  H2/O2,  H2/air  and  CH4/O2 
mixtures  in  an  open  flow  system  initially  at  1  atm  and  room  temperature  via 
an  ArF  excimer  laser  (193  nm)  photolysis  of  a  small  amount  of  NH3  doped  in 
the  flow  mixtures.  The  ignition  appears  to  be  homogenous  since  the 
ignition  delay  times  measured  at  several  locations  are  close  in  value.  The 
minimum  Ignition  energy  density  is  measured  to  be  137  ±  8,  190  ±  20  and  300 
±  30  mJ/cm3  for  stoichiometric  mixtures  of  H2/02/NH3,  H2/air/NH3  and 
CH4/O2/NH3  flows,  respectively,  and  appears  to  be  insensitive  to  the  fuel 
equivalence  ratio  between  0.35  and  3.0.  The  ignition  delay  time  depends 
strongly  on  the  initial  NH3  concentration  and  the  laser  energy  deposition 
density,  however  is  also  nearly  Invariant  to  the  equivalence  ratio.  The 
required  minimum  ignition  energy  density  appears  to  be  substantially  lower 
than  that  expected  for  thermal  Ignition.  Preliminary  kinetic  modeling 
calculations  suggest  that  the  hot  H  atoms  and/or  electronically  excited 
NH2 (A)  from  the  photolysis  of  NH3  may  play  an  important  role  In  enhancing 
the  Ignition. 

I.  Introduction 

Photochemical  Ignition  of  combustion  provides  an  opportunity  to  study 
details  of  the  associated  chemical  processes,  since  energy  can  be  directed 
homogeneously  and  within  a  short  duration  into  a  specific  mode  of  a 
reactant,  to  prepare  a  well-defined  Initial  condition.  Earlier  experiments 
by  Farkas  et  al.  [1]  have  shown  that  H2/O2  mixtures  with  a  small  amount  of 
NH3  become  explosive  at  a  moderate  temperature  (~420°C)  under  ultraviolet 
radiation.  The  reaction  appears  to  proceed  much  more  rapidly  via  NH3 
photolysis,  rather  than  by  other  means  of  photogeneration  of  H  atoms,  such 


Rl-003-89 


IV-1 


as  via  photosensitized  mercury  atoms  that  dissociate  H2  into  2H  [2].  It 
has  been  speculated  that  NH2  may  play  an  important  role  in  the  ignition 
chemistry;  however,  the  detailed  reaction  mechanism  remains  poorly 
understood  [3]. 

Recent  advances  in  high-energy  ultraviolet  lasers  have  renewed 
interest  in  photochemical  ignition,  since  these  lasers  can  generate 
sufficient  radicals  by  photolysis  to  initiate  combustion  in  an 
exceptionally  controlled  way.  Laser-spark  has  been  used  to  ignite  fuel- 
oxidizer  mixtures,  either  by  gaseous  breakdown  [4,5]  or  multiphoton 
processes  [6,7],  within  a  tightly  focused  spot.  Absorption  of  laser 
photons  which  produce  reactive  radical  species  in  the  absence  of  spark 
formation  has  also  been  used  to  initiate  combustion.  For  example,  Lavid 
et  al.  [8]  have  used  a  F2  laser  at  157  nm  to  dissociate  O2  to  ignite  H2/O2 
mixtures,  and  Lucas  et  al .  [9]  have  used  a  KrF  laser  at  248  nm  to 
dissociate  O3  to  ignite  H2/O2-,  CH4/O2-,  and  C3H8/O2  mixtures.  We  are 
primarily  Interested  in  the  latter  kind  of  photochemical  Ignition.  The 
related  subjects  of  Ignition  by  the  addition  of  radicals  have  been  studied 
theoretically  by  Gulrguis  and  Oppenheim  [10],  -loane  [11-13],  and  Sloane 
and  Schoene  [14]. 

In  this  paper,  we  describe  a  study  of  Ignition  initiated  by  photolysis 
of  NH3  by  the  use  of  an  ArF  laser  at  193  nm  In  mixtures  of  H2/O2/NH3  and 
H2/a1r/NH3  and  CH4/O2/NH3.  We  have  also  performed  kinetic  modeling 
simulation  on  the  ignition  of  H2/O2  by  radical  addition.  Our  kinetic 
modeling  calculations  suggest  that  the  highly  energetic  photo-fragments, 
such  as  hot  H  atoms  and  electronically  excited  NH2 (A) ,  may  play  an 
important  role  in  the  initiation  of  combustion. 
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II.  Review  of  Relevant  Photochemistry 


Light  of  193  nm  wavelength  excites  NH3  into  v^-  6  of  the  ft*A 2"  «•  XUj ' 

transition  [15,16],  and  the  excited  NH3  dissociates  into  NH2  and  H  with 
nearly  unit  quantum  efficiency  [17,18].  The  molar  extinction  coefficient 
for  NH3  at  193  nm  is  measured  to  be  3250  liter/mole-cm  [19],  which  Is  much 
higher  than  that  of  O2  [20].  The  absorption  by  O2  is  estimated  to  be  less 
than  about  one  percent  of  that  by  NH3  under  our  experimental  conditions. 

The  bond  dissociation  energy  for  NH3,  D0(H-NH2),  is  -446.4  kJ/mole  [21,22]; 
hence  there  is  about  171.5  kJ/mole  residual  energy  left  from  the  193  nm 
photons.  Based  on  the  distribution  of  the  translational  energy  of  the  H 
atom  photofragments  from  the  photolysis  of  NH3,  Biesner  et  al .  [23] 
conclude  that  about  21  percent  of  the  residual  energy  is  in  the  form 
translational  excitation  and  the  rest  is  In  the  excitation  of  NH2.  These 
investigators  also  derive  a  branching  ratio  of  0.33  for  the  yield  of  NH2(A) 
to  NH2(X).  However,  a  lower  branching  ratio  (0.025)  was  reported  earlier 
by  Donnelley  et  al ,  [24]  based  on  the  NH2  emission  following  the  193  nm 
photolysis  of  NH3. 

At  higher  laser  intensity,  NH(A)  can  also  be  produced  by  a  sequential 
two-photon  absorption  process  via  an  Internally  excited  Intermediate 
state  of  NH2 (^)  [25,26].  The  molar  extinction  coefficient  for  the 
absorption  of  the  intermediate  NH2(X)  at  193  nm  Is  measured  to  be 
-10.4  Z  mole”1  cm-1  by  Hofzumahaus  et  al.  [27]. 

Quenching  of  the  Internally  excited  photolysis  products  of  NH2(X)  by 
02  and  H2  appear  to  be  relatively  slow  [25].  This  implies  that  chemical 
reactions  with  the  two  collision  partners  are  not  an  important  channel  for 


Rl-003-89 


IV-3 


the  decay  of  NH2 (X) -  The  photolysis  products  of  NH2 (A)  appear  to  have  two 
components:  ~85%  In  a  short-lived  component  and  ~15%  in  a  long-lived 
component  [28].  The  former  is  postulated  to  be  pure  NH2(A),  while  the 
latter  to  be  a  mixed  state  of  NH2 (A)  and  NH2 (X) .  The  long-lived  component 
Is  not  very  reactive  with  O2  [25].  However,  the  reactivity  of  the  short¬ 
lived  component  with  O2  has  not  been  reported  in  the  literature. 

The  reaction  of  NH2  and  O2  has  been  recently  reviewed  by  Hanson  et  al . 
[29].  Fujii  et  al .  [30]  proposed  that  the  reaction  proceed  as  follows, 

NH2  +  O2  *  HNO  +  OH.  (1) 

The  rate  constant  is  determined  to  be  1012*25  exp(-15,000/RT)  with  an 
activation  energy  of  ~63  kJ/mole.  The  NH2(A)  state  has  an  excitation 
energy  of  122  kJ/mole  above  the  NH2 (X)  ground  state.  It  is  possible 
that  the  short  lived  component  of  the  NH2$)  from  the  photolysis  of  NH3  may 
enhance  the  rate  of  Reaction  1.  However,  quenching  cross  sections  for  the 
short-lived  component  are  generally  quite  high  with  collisional  partners, 
Including  NH3,  H2,  Ar  and  He  [25,31],  the  loss  of  the  short-lived  component 
is  likely  to  be  dominated  by  collisional  quenching  rather  than  by  the 
reaction  with  O2. 

The  hot  H  atoms  from  the  photolysis  of  NH3  may  enhance  the  rate  of 
reaction  with  O2  [32-36]: 

H  +  O2  ■*  OH  +  0.  (2) 

According  to  a  theoretical  study  by  Miller  [32],  the  reaction  cross 

section  has  a  threshold  energy  of  about  52  kJ/mole  and  rises  gradually 

0 

to  a  broad  maximum  of  ~0.4  A2  for  translation  energies  between  105  and  209 
kJ/mole.  The  distribution  of  H  atom  kinetic  energy  from  the  photolysis  of 
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NH3  has  been  measured  by  Biesner  [22,23].  From  this  one  estimates  that  ~20 
percent  of  the  H  atoms  should  have  a  kinetic  energy  above  the  threshold. 
Using  the  measured  H  atom  distribution  profile  and  the  calculated  reaction 
cross  section  [32] ,  one  derives  an  effective  rate  constant  (k2)  of  ~9.4  x 
10*2  cm^/mole-s  for  those  H  atoms  with  translational  energy  above  the 
threshold.  A  major  competing  reaction  for  the  H  atoms  is  Reaction  3,  which 
may  be  considered  as  a  chain-breaking  reaction,  since  HO2  is  not  very 
reactive: 

H  +  02  +  M  -»  H02  +  M.  (3) 

The  product  of  the  rate  constant  of  k3  [37]  and  the  gas  density  (M)  is 
calculated  to  be  about  3.6  x  10*2  cm^/mole-s  at  one  atmosphere.  This 
Indicates  that  ~20  percent  of  the  H  atoms  produced  from  NH3  photolysis  may 
react  with  02  favorably  via  Reaction  2. 

III.  Experimental  Procedure 

The  experimental  arrangement  Is  shown  In  Figure  1.  A  Lambda  Physik 
Exclmer  Laser  (150  EST),  which  provides  ArF  laser  radiation  at  193  nm,  Is 
used  as  the  photolysis  source.  This  exclmer  laser  system  consists  of  two 
laser  units.  The  output  from  an  unstable  resonator  of  the  first  unit  is 
amplified  through  the  second  unit  to  achieve  a  high  energy  output  with  a 
low  beam  divergence.  A  UV  grade  MgF2  lens  with  a  50  cm  focal  length  is 
used  to  deliver  the  beam  with  a  rectangular  profile  of  about  2.1  mm  wide  x 
7.1  mm  high  at  the  leading  edge  and  1.2  mm  wide  x  3.3  mm  high  at  the 
trailing  edge  of  a  6-cm  diameter  flat  flame  burner  (McKenna  Products).  The 
focus  of  the  laser  beam  Is  far  past  the  burner  to  avoid  any  gaseous 
breakdown  In  the  mixtures.  The  laser  beam  size  Is  determined  by  burn  spots 
on  Polarlod  films. 
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The  Incident  as  well  as  the  transmitted  laser  energies  are  monitored 
by  the  use  of  quartz  flat  beam  splitters  and  laser  energy  ratiometer  probes 
(Laser  Precision  Model  7200).  A  typical  incident  laser  energy  is  in  the 
range  of  100-290  mJ  per  pulse.  The  OH  emission  is  monitored  at  three 
locations  along  the  beam  path:  2.2  cm,  3.0  cm  and  3.8  cm  from  the  leading 
edge  of  the  burner.  The  center  location  i^3  cm)  is  observed  with  a 
monochromator  and  photomultiplier  tube,  whereas  the  other  two  locations  are 
monitored  with  a  band-pass  filter  (at  310  nm  with  10  nm  bandwidth)  and 
photomultiplier. 

The  premixed  gas  mixtures  which  are  prepared  in  a  stainless  steel 
manifold  are  flowed  through  the  flat  flame  burner.  Coaxial  N2  shroud  gas 
is  used  to  minimize  any  mixing  of  the  combustible  flow  mixtures  with  the 
surrounding  air.  Matheson  anhydrous  grade  ammonia  (99.99%),  hydrogen 
(99.995%)  and  O2  (99.5%)  are  used.  The  total  gas  flow  is  3.75  £/min  for 
H2/O2/NH3  and  CH4/O2/NH3,  and  8.75  £./min  for  H2/air/NH3  mixtures.  These 
correspond  to  cold  gas  linear  velocities  of  1.3  and  3.0  mm/mi n, 
respectively.  Ignition  experiments  are  conducted  on  the  premixed  gases  in 
an  open  system  at  ~1  mm  immediately  above  the  flat  flame  burner.  All  the 
data  are  collected  at  1  atm  and  at  room  temperature. 

Majority  of  the  experiments  are  carried  out  under  conditions  that  the 
incident  excimer  laser  energy  is  sufficiently  high  to  bleach  out  the  NH3 
absorption  along  the  irradiated  pathlength,  but  is  not  too  high  to  cause 
any  significant  multiphoton  effects.  This  should  yield  a  relatively 
uniform  laser  energy  deposition  density  and  a  well  defined  production  of 
Initial  radical  density  from  the  photolysis  in  spite  of  any  non-uniformity 
or  hot  spots  In  the  laser  beam  profile. 
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We  have  made  numerous  observations  with  an  optical  multichannel 
analyzer  (OMA).  The  emission  during  ArF  laser  Irradiation  is  found  to  be 
mainly  composed  of  the  NH  band  at  336  nm,  NH2  at  500-800  nm  and  OH  bands  at 
308  nm  and  283  nm.  No  atomic  or  ionic  emission  lines  indicative  of  plasma 
formation  are  observed.  This  precluse  ignition  caused  by  any  laser  induced 
gaseous  breakdown. 

IV.  Results  and  Discussion 
H?/Q?/NH^  Mixtures 

The  results  for  the  ignition  of  H2/O2/NH3  mixtures  by  photolysis  of 
NH3  are  summarized  in  Figures  2-7.  Figures  2a-2c  show  typical  OH  emission 
temporal  profiles  in  three  different  time  scales  to  illustrate  various 
events.  Shown  in  Figure  2a  is  a  prompt  OH  emission  with  a  rapid  rise  and 
fall  which  Is  nearly  coincident  with  the  excimer  laser  pulse  (which  also 
has  two  maxima).  The  prompt  OH  emission  Is  present  Independently  of 
whether  the  mixture  is  Ignited  or  not.  Figure  2b  shows  that  the  OH 
emission  decreases  nearly  to  the  background  level  after  the  prompt  OH 
emission,  and  then  increases  rapidly  after  a  certain  delay  time.  Figure  2c 
shows  that  the  delayed  OH  emission  reaches  a  maximum  and  then  decreases  to 
a  relatively  low  steady-state  value.  The  peak  intensity  of  the  prompt  OH 
Is  In  this  case  about  5  times  of  that  of  the  delayed  OH  emission  which  In 
turn  Is  about  an  order  of  magnitude  higher  than  the  steady  state  value  at 
the  late  time. 

The  presence  of  the  strong  prompt  OH  emission  Implies  that  the  initial 
photolysis  products  may  be  converted  rapidly  into  0  and  OH,  since  the  OH 
emission  Is  formed  mainly  via  a  three  body  recombination  of  H  +  0  +  M  or  an 
energy  transfer  from  some  excited  species  to  the  ground  OH  molecules.  It 
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Is  known  that  thermal  H  atoms  and  NH2 (X)  are  not  reactive  with  O2  at  a 
relatively  low  temperature,  this  suggests  that  the  hot  H  atoms  and  possibly 
NH£(A)  from  the  photolysis  of  NH3  may  indeed  react  with  O2  rapidly  to 
convert  into  the  more  reactive  species  of  0  and  OH  as  discussed  in 
Section  II. 

The  delayed  OH  emission  as  shown  in  Figure  2b  is  indicative  of 
combustion  since  it  occurs  only  if  the  mixture  is  ignited.  An  ignition 
delay  time  may  therefore  be  assigned  based  on  the  delayed  emission  by 
extrapolating  the  linear  portion  of  the  OH  emission  rise  to  the  zero 
intensity  base  line.  Figures  3-5  summarize  the  measured  ignition  delay 
times  as  function  of  incident  laser  energy  monitored  at  three  different 
locations  for  various  initial  NH3  mole  fractions.  Several  noticeable 
features  are  summarized  below:  (1)  The  ignition  delay  times  at  the  three 
locations  appear  in  general  to  be  relatively  close  in  value  for  the 
Incident  excimer  laser  energy  above  that  required  to  bleach  out  the  NH3 
absorption  as  shown  for  all  the  data  in  Figures  3  and  4  and  part  of  the 
data  with  incident  laser  energy  above  the  bleaching  threshold  (~155  mJ)  in 
Figure  5.  The  variation  in  the  ignition  delay  times  among  the  three 
locations  under  bleaching  conditions  are  not  enough  to  be  accounted  for  by 
flame  propagation  from  a  particular  ignition  site.  This  suggests  that  the 
ignition  sites  are  probably  well  distributed  In  the  entire  irradiated 
volume.  The  variation  may  only  reflect  some  statistical  instability  of  the 
local  Ignition  events.  (2)  The  ignition  delay  time  is  in  general 
relatively  Insensitive  to  the  Incident  laser  energies  for  those  above  the 
bleaching  threshold  as  shown  In  Figures  4  and  5.  This  is  probably  expected 
since  the  laser  energy  deposition  should  remain  constant  under  bleaching 
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conditions  Irrespective  of  increases  in  the  Incident  laser  energy.  However 
this  Independency  on  the  Incident  laser  energy  appears  to  break  down  more 
noticeably  for  a  lower  NH3  mole  fraction  as  shown  in  Figure  3.  Here  the 
ignition  delay  time  decreases  more  substantially  with  increasing  incident 
laser  energy.  The  cause  for  this  may  be  attributed  to  the  fact  that  as  the 
NH3  concentration  Is  reduced  to  that  near  a  critical  value  needed  for 
Ignition,  the  Ignition  becomes  more  sensitive  to  any  multiphoton  processes. 
One  of  the  Important  multiphoton  processes  Is  a  sequential  two-photon 
absorption  which  dissociates  NH3  into  NH(A)  and  2  H  as  discussed  in  Section 
II.  The  quantum  yield  for  this  process  is  estimated  to  be  ~  1.8  x  10~3  at 
140  mJ  and  ~  3.5  x  10~3  at  260  mJ  of  the  incident  laser  energy.  Although 
these  values  are  only  moderate  as  compared  to  the  nearly  unity  for  the 
single  photon  process.  It  is  conceivable  that  any  mildly  Increase  In  the 
radical  concentration  may  drastically  enhance  the  ignition  process  for 
those  with  low  NH3  mole  fraction  near  that  required  for  ignition. 

Figure  6  shows  the  ignition  delay  time  as  function  of  Initial  NH3  mole 
fraction  for  several  fuel  equivalence  ratios.  Also  shown  in  the  abscisa  is 
the  laser  energy  deposition  density  which  is  derived  experimentally  from 
the  difference  In  the  Incident  and  transmitted  laser  energy  divided  by  the 
total  irradiated  volume.  The  derived  laser  energy  deposition  density 
agrees  well  with  that  expected  to  bleach  out  the  NH3  Initially  In  the 
mixtures.  The  Ignition  delay  time  plotted  here  Is  an  average  value  for  all 
the  data  with  incident  laser  energies  above  a  bleaching  threshold.  Some 
ambiguity  may  remain  for  those  events  with  low  Initial  NH3  mole  fraction  by 
using  the  average  value  due  to  the  observed  variation  with  Incident  laser 
energy.  This  ambiguity  however,  appears  to  be  relatively  unimportant  in 
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view  of  the  steeper  variation  of  the  ignition  delay  time  with  the  initial 
NH3  mole  fraction.  The  error  bar  given  for  the  cases  with  the 
stoichiometric  ratio  =  1.0)  Indicates  one  standard  deviation  in  the 
spread  of  the  measured  data  points.  The  ignition  delay  times  which  depends 
strongly  with  the  Initial  NH3  mole  fraction,  is  however,  nearly  Independent 
of  the  fuel  equivalence  ratio  between  0.35  and  3.0  within  the  experimental 
uncertainty.  Also  shown  in  Figure  6  Is  an  approximate  threshold  value  of 
Initial  NH3  mole  fraction  or  laser  energy  deposition  density  needed  for 
ignition.  Below  this  threshold  no  ignition  occurs.  Figure  7  summarizes 
the  derived  minimum  NH3  mole  fraction  and  the  corresponding  minimum  laser 
energy  deposition  density  needed  for  ignition  for  several  equivalence 
ratios.  It  Is  noted  that  the  required  minimum  laser  energy  deposition 
density  Is  also  Insensitive  to  the  fuel  equivalence  ratio  between  0.35  and 
3.0.  We  have  also  performed  additional  experiments  using  a  larger  laser 
beam  with  a  nominal  size  of  ~2.1  mm  x  8.5  mm.  The  ignition  delay  time  and 
the  minimum  laser  energy  deposition  density  for  this  condition  are  not 
significantly  different  from  that  shown  above. 

The  observation  that  the  minimum  Ignition  energy  density  and  Ignition 
delay  time  are  Insensitive  to  the  equivalence  ratio  is  of  interest.  Lucas 
et  al.  [9]  have  also  shown  that  the  minimum  laser  energy  deposition  density 
Is  insensitive  to  the  equivalence  ratio  between  0.7  and  3.3  by  the 
photolysis  of  O3,  and  Lavid  et  al.  [8]  have  shown  a  minimum  near  0.6  and 
increasing  only  moderately  at  0.4,  1.0  and  1.6  by  the  photolysis  of  O2.  On 
the  contrary,  the  minimum  Ignition  energy  Increases  very  drastically  on 
both  sides  of  the  equivalence  ratios  from  that  near  the  stoichiometric 
ratio  for  ignition  Induced  by  electric  sparks  [3],  laser-sparks  [4-5]  and 
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multlphoton  excitation  processes  [6-7].  One  of  the  main  differences  here 
Is  that  the  radicals  are  produced  at  a  relatively  low  temperature.  In  a 
theoretical  study,  Sloane  [11]  has  predicted  that  the  minimum  Ignition 
energy  is  nearly  insensitive  to  the  equivalence  ratio  between  0.4  and  1.2 
for  the  ignition  of  Cfy/air  by  the  addition  of  atomic  oxygen.  This 
suggests  that  a  photochemical  method  or  other  suitable  means  of  radical 
additions  may  be  useful  for  ignition  of  extreme  lean  on  rich  fuel  mixtures 
which  may  not  be  readily  achieved  by  a  conventional  thermal  or  spark 
ignition. 

Maas  and  Warnatz  [38]  have  recently  modeled  the  thermal  Ignition  of 
H2/O2  premixtures  and  predicted  that  a  minimum  energy  density  of  ~430 
mJ/cm3  needed  for  Ignition  of  a  stoichiometric  mixture  Initially  at  1  atm 
and  room  temperature.  The  observed  threshold  for  ignition  (137  ±  8  mJ/cm3) 
by  the  NH3  photolysis  is  significantly  lower  than  that  predicted  for  the 
thermal  ignition.  This  suggests  that  the  method  of  radical  addition  by 
photolysis  may  be  more  efficient  in  the  initiation  of  combustion.  The 
radical  pair  density  of  NH2  +  H  needed  for  ignition  by  the  photolysis  of 
NH3  is  estimated  to  be  ~  1.3  x  1017  cm"3.  This  value  is  very  close  to  the 
atomic  oxygen  concentration  needed  for  ignition  based  on  the  photolysis  of 
O3  reported  by  Lucas  et  al .  [9].  This  implies  that  the  photolysis  products 
of  NH2  and  H  are  as  effective  as  the  reactive  0  atoms  in  initiating 
combustion.  Our  kinetic  modeling  calculations,  to  be  described  later, 
Indicate  that  if  the  NH2  and  H  were  thermal  1  zed  before  reaction, 
ignition  should  not  have  occurred.  This  suggests  that  the  hot  H  atoms 
and  NH2$)  produced  from  the  NH3  photolysis  may  enhance  the  Initiation 
of  combustion. 
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H?/a1r/NH-3  Mixtures 

The  air  flow  Is  simulated  by  the  use  of  an  appropriate  ratio  of  O2  and 
N2  flows.  The  results  for  the  Ignition  of  H2/a1r/NH3  stoichiometric  flow 
mixtures  by  NH3  photolysis  are  summarized  in  Figure  8.  The  ignition  delay 
time  plotted  here  Is  an  average  value  for  all  the  data  with  Incident  laser 
energies  above  the  bleaching  threshold.  The  error  bar  indicates  one 
standard  deviation.  The  Ignition  delay  time  depends  strongly  on  the 
Initial  NH3  mole  fraction  or  the  laser  energy  deposition  density.  In 
comparing  to  the  H2/O2/NH3  mixtures,  the  minimum  laser  energy  deposition 
density  needed  for  ignition  Is  higher,  and  the  Ignition  delay  time  is 
longer  for  the  same  energy  deposition  density.  This  Is  probably  expected 
due  to  the  Increase  in  the  total  heat  capacity  with  excess  N2.  However, 
the  Ignition  threshold  appears  to  Increase  only  moderately  from  137  ±  8 
mJ/cm3  for  the  H2/O2/NH3  to  ~190  ±  20  mJ/cm3  for  the  H2/air/NH3  flow 
mi xtures . 

CHa/0?/NHi  Mixtures 

The  results  for  the  Ignition  of  CH4/O2/NH3  stoichiometric  flow 
mixtures  by  NH3  photolysis  are  summarized  in  Figure  9.  The  Ignition  delay 
time  plotted  here  Is  an  average  value  for  all  the  data  with  incident  laser 
energies  above  the  bleaching  threshold.  The  error  bar  Indicates  one 
standard  deviation.  The  Ignition  delay  depends  strongly  with  the  initial 
NH3  mole  fraction  or  the  laser  energy  deposition  density.  The  minimum  NH3 
mole  fraction  needed  for  Ignition  Is  measured  to  be  1.2  ±  0.1%  and  the 
corresponding  minimum  laser  energy  deposition  density  Is  300  ±  30  mJ/cm3. 
These  values  are  higher  compared  to  that  for  the  H2/O2/NH3  mixtures. 
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Kinetic  Modeling  of  H?/0?/NHt  Ignition 

Ignition  of  H2/O2/NH3  mixtures  by  NH3  photolysis  is  simulated  by  using 
the  zero-dimensional  version  of  the  HCT  code  of  Lund  [39],  and  also  the 
CHEMKIN  code  of  Kee  et  al.  [40].  The  zero  dimensional  modeling  which 
assumes  no  loss  of  energy  and  species  out  of  the  irradiated  volume  should 
yield  an  upper  bound  on  the  Ignition  delay  time  and  the  mixture 
ignltability.  The  H-O-N  reactions  and  rate  constants  are  adopted  from 
those  complied  with  the  codes.  Both  codes  yield  similar  results  for  all 
the  conditions  reported  here. 

The  experimental  results  with  an  Initial  NH3  mole  fraction  of  0.012 
under  bleaching  conditions  are  considered  for  modeling  comparison.  Table  1 
summarizes  the  results  for  four  different  Input  conditions.  The  gas 
temperature  Is  estimated  by  assuming  that  all  the  available  residual  energy 
of  the  photon  energy  after  breaking  the  N-H  bond  from  the  photolysis  of  NH3 
Is  utilized  In  heating  the  gas  mixture  initially  at  room  temperature.  This 
estimated  gas  temperature  should  be  an  upper  bound  value  since  some  of  the 
residual  energies  are  known  to  retain  in  the  photolysis  products  of  hot  H 
atoms  and  Internally  excited  NH2. 

In  case  1,  if  all  the  Initial  NH3  molecules  are  assumed  to  be 
photolyzed  Into  H  and  NH2  and  that  the  energetics  In  the  photolysis 
products  are  quenched  before  reaction  proceeds,  one  predicts  no  Ignition 
within  1  second  time  span.  This  may  be  explained  by  the  fact  that  thermal 
H  atoms  are  likely  to  recombine  with  O2  to  form  unreactive  HO2  and  that  NH2 
are  not  reactive  with  O2  at  a  relatively  low  temperature.  This  calculation 
suggests  that  the  energetics  contained  In  the  photolysis  products  should  be 
considered  to  explain  the  experimental  observation. 
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As  discussed  In  Section  II,  20%  of  the  H  atoms  produced  from  the  NH3 
photolysis  [22-23]  should  have  translational  energy  above  the  threshold  for 
Reaction  2  [28].  In  case  2,  20%  of  the  H  atoms  are  assumed  to  convert  into 
OH  and  0.  However  no  Ignition  Is  predicted  within  1  second  time  span. 

Also  discussed  In  Section  II,  the  short-lived  component  of  the  NH2(A)  from 
NH3  photolysis  is  likely  to  be  collislonally  quenched  rather  than  to  react 
with  03.  Nevertheless,  In  case  3,  21%  yield  of  NH2  [23]  are  assumed  to 
convert  Into  HNO  and  OH  via  Reaction  1.  This  yields  a  predicted  Ignition 
delay  time  of  1.1  ms  which  Is  too  long  compared  to  the  observed  value  of  27 
±  8  /is.  In  case  4,  if  all  the  H  atoms  from  the  NH3  photolysis  are  assumed 
to  convert  Into  OH  and  0  via  Reaction  2,  one  calculates  an  ignition  delay 
time  of  130  /is  which  is  still  too  long  compared  to  the  observed  value. 
Although  we  have  not  been  able  to  predict  an  Ignition  delay  time  In 
agreement  with  that  observed  even  under  a  very  optimistic  condition,  these 
calculations  suggest  that  the  energetics  in  the  photolysis  products, 
Including  hot  H  atoms  and  possibly  NH2 (A) ,  may  play  an  Important  role  to 
enhance  the  initiation  of  the  combustion. 

In  addition,  we  have  Investigated  the  dependence  of  the  Ignition-delay 
time  on  both  the  fuel -equivalence  ratio  and  the  Initial  radical 
concentration,  by  computer  calculations.  We  again  have  used  the  Chemkin 
code  for  a  batch-process  (zero-spatlal-dimenslon)  model.  The  Initial 
concentrations  of  0,  OH,  and  NH2  are  taken  to  be  Identical.  This  choice 
Implies  that  all  of  the  H  atoms  from  the  photolysis  of  NH3  are  converted 
rapidly  Into  0  and  OH.  Although  this  may  be  an  oversimplifying  assumption, 
It  serves  the  objective  of  predicting  the  variation  of  the  Ignition-delay 
time  on  the  equivalence  ratio  and  the  initial  radical  concentration. 
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The  calculated  Ignition-delay  time  (Fig.  10)  for  H2/O2  mixtures, 
initially  at  1  atm  and  369  K,  depends  quite  sensitively  on  the  initial 
radical  concentration,  as  we  have  observed  (Fig.  6).  Furthermore,  the 
calculated  Ignition-delay  time  appears  to  be  relatively  insensitive  to  the 
fuel -equivalence  ratio,  again  as  we  have  observed.  Within  the  experimental 
uncertainty,  we  cannot  resolve  the  slight  variation  predicted. 

Since  the  ignition-delay  time  Is  relatively  insensitive  to  the 
equivalence  ratio,  we  expect  that  the  minimum  ignition-energy  density  also 
should  be  insensitive  to  the  equivalence  ratio.  This  is  Indeed  what  Is 
observed  (Fig.  7).  For  an  electric  spark,  the  minimum  energy  Increases 
drastically  for  both  fuel-lean  and  fuel-rich  mixtures.  A  major  reason  for 
the  very  different  ignition  behavior  for  the  two  methods  of  ignition  may 
concern  the  fact  that  the  electric  spark  creates  a  small  ignition  kernel, 
which  is  strongly  Influenced  by  gas-transport  phenomena,  as  now  discussed. 

According  to  Lewis  and  von  Elbe  [3],  spark  ignition  occurs  only  if  the 
energy  added  to  the  gas  Is  sufficient  to  heat  a  volume  of  the  gas  of  about 
the  thickness  of  steadily  propagating  flame  to  the  adiabatic  flame 
temperature.  By  relating  the  flame  thickness  to  the  laminar  burning 
velocity  Sl,  the  minimum  spark-ignition  energy  may  be  given  approximately 
as  [41,  42] 

Ei»1n  =  <4/3>  [k/V  SL>31  f  Cp  <T.  -  V' 
where  k  is  the  thermal  conductivity,  Cp  Is  the  specific  heat,  p  is  the  gas 
density,  T«  Is  the  adiabatic  flame  temperature,  and  T0  is  the  Initial  gas 
temperature.  This  equation  Indicates  a  strong  dependence  of  the  minimum 
ignition  energy  on  the  laminar  burning  velocity.  For  both  H2/O2  and  CH4/O2 
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mixtures,  the  laminar  burning  velocity  Is  maximum  near  the  stoichiometric 
mixture,  and  falls  off  rapidly  away  from  that  mixture.  This  observation 
qualitatively  explains  the  sensitive  dependence  of  spark-ignition  energy  on 
equivalence  ratio.  Recently  Sloane  [13]  and  Sloane  and  Schoene  [14]  have 
modeled  the  Ignition  of  a  small  volume  of  thermally  heated  Cfy/air  mixture, 
and  have  qualitatively  explained  spark-ignition  behavior.  Maas  and  Warnatz 
[38]  also  have  predicted  the  minimum  Ignition-energy  density  for  H2/O2 
mixtures  when  thermal  energy  is  added  to  a  volume  of  about  1-mm  radius 
within  a  relatively  large  volume  of  gas  mixture.  Their  work  shows  that  the 
spark-ignition  energy  Is  strongly  Influenced  by  gas-transport  phenomena, 
and  that  the  Ignition  energy  becomes  strongly  dependent  on  equivalence 
ratio  when  the  Ignition-spot  size  Is  small.  When  the  spot  size  Is  larger, 
the  Ignition  Is  Influenced  strongly  by  chemistry,  and  the  minimum  Ignition- 
energy  density  becomes  Insensitive  to  the  equivalence  ratio,  independently 
of  whether  the  initial  energy  is  added  in  the  form  of  heat  or  of  radical 
dissociation. 

V.  Conclusion 

We  have  achieved  volumetric  homogeneous  ignition  of  H2/O2/NH3, 
H2/air/NH3  and  CH4/O2/NH3  mixtures  In  an  open  flow  system  initially  at  1 
atm  and  room  temperature  via  photolysis  of  NH3.  The  emission  spectrum 
monitored  by  an  OMA  shows  no  discernable  atomic  or  ionic  emission, 
precluding  Ignition  caused  by  gaseous  breakdown.  The  ignition  delay  times 
based  on  the  308  nm  OH  emission  band  at  three  locations  appear  to  be  close 
in  value,  suggesting  the  presence  of  well  distributed  Ignition  sites  within 
the  entire  irradiated  volume.  The  Ignition  delay  time  depends  strongly  on 
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the  Initial  NH3  mole  fraction  and  nearly  Independent  of  the  fuel 
equivalence  ratio  between  0.35  and  3.0  In  H2/O2/NH3  mixtures.  The  minimum 
Ignition  energy  density  is  also  Insensitive  to  the  equivalence  ratio. 
Furthermore,  the  required  minimum  ignition  energy  density  (137  ±  8  mJ/cm3) 
appears  to  be  substantially  less  than  that  by  a  thermal  process  (~  430 
mJ/cm3)  for  a  stoichiometric  H2/O2  mixture. 

Our  kinetic  modeling  calculations  show  that  If  the  hot  H  atoms  and 
NH2 (A)  from  the  photolysis  of  NH3  were  collislonally  quenched  rapidly, 
Ignition  should  not  have  occurred.  This  implies  that  the  excitation  energy 
contained  in  the  initial  photolysis  products  may  play  an  Important  role  in 
enhancing  the  ignition.  However,  the  detailed  mechanism  leading  to  the 
enhanced  Ignition  with  low  required  Ignition  energy  density  and  short 
ignition  delay  time  remains  not  well  understood. 

We  acknowledge  E.Y.  Wang  (TRW)  and  J.Y.  Chen  (Sandia  National 
Laboratories)  for  assistance  In  computer  modeling  calculations  and  T. 

Sloane  (General  Motors)  for  helpful  discussions.  This  work  was  sponsored 
by  the  Air  Force  Office  of  Scientific  Research,  Air  Force  System  Command, 
USAF,  under  Contract  F49620-87-C-0081. 
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Table  1.  Computer  Modeling  of  H2/O2  Ignition  by  NH3  Photolysis(a) 


Case 

No. 

Gas  Composition 
(mole  fraction) 

Gas 

Temperature 

Ignition- 
Delay  Time 

Remarks  On 

Gas  Composition 

1 

H2  0.651 

02  0.325 

H  0.012 

NH2  0.012 

369 

>  1  s 

Energetics  in  the 
photolyzsis 
products  of  H  and 
NH?  were 
collisionally 
quenched. 

2 

H2  0.651 

02  0.324 

H  0.0096 

NH2  0.0012 

0  0.0024 

OH  0.0024 

369 

>  1  S 

20%  of  H  reacted 
with  02  to  form 

OH  +  0. 

3 

H2  0.648 

02  0.323 

H  0.0096 

NH2  0.0094 

0  0.0024 

OH  0.0049 

HN0  0.0025 

369 

1.1  ms 

20%  of  H  and  21% 
of  NH2  reacted 
with  02  to  form 

OH  +  0  and  HN0  + 

0,  respectively. 

4 

H2  0.643 

02  0.321 

NH2  0.012 

O  0.012 

OH  0.012 

369 

130  ns 

All  H  atoms 
reacted  with  02 
to  form  OH  +  0. 

(a)  Assuming  that  the  Initial  NH3  mole  fraction  Is  0.012  and  all  NH3 
molecules  are  photolyzed  into  H  +  NH2  under  bleaching  condition. 
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Figure  Captions 


1.  Schematic  of  experimental  arrangement  for  the  ignition  of  premixed 
combustible  flows  by  ArF  excimer  laser  photolysis  of  NH3,  M:  mirror, 
QF:  quartz  flat  and  PMT:  photomultiplier  tube. 

2.  Typical  OH  Emission  profiles  at  three  different  time  scales  to  show 
various  events  following  ArF  excimer  laser  irradiation:  (a)  prompt  OH 
emission,  (b)  rise  of  delayed  OH  emission,  and  (c)  rise  and  fall  of 
the  delayed  OH  emission. 

3.  Ignition  delay  time  in  stoichiometric  H2/O2/NH3  mixtures  with  0.0065 
mole  fraction  of  NH3  as  function  of  incident  laser  energy  based  on  the 
OH  emission  measurements  at  three  locations:  (V)  2.2  cm,  (0)  3.0  cm, 
and  (A)  3.8  cm  downstream  of  the  burner  leading  edge. 

4.  Ignition  delay  time  In  H2/O2/NH3  stoichiometric  mixtures  with  0.0086 
mole  fraction  of  NH3  as  function  of  Incident  laser  energy  based  on  the 
OH  emission  at  three  locations:  (7)  22  cm,  (0)  3.0  cm  and  (A)  3.8  cm 
downstream  of  the  burner  leading  edge. 

5.  Ignition  delay  time  In  H2/O2/NH3  stoichiometric  mixtures  with  0.0135 
mole  fraction  of  NH3  as  function  of  incident  laser  energy  based  on  OH 
emission  at  three  locations:  (7)  2.2  cm,  (0)  3.0  cm  and  (A)  3.8  cm 
downstream  of  the  burner  leading  edge.  The  arrow  indicates  the 
minimum  incident  laser  energy  needed  to  bleach  out  the  NH3  absorption 
up  to  the  3.8  cm  location. 

6.  Ignition  delay  time  In  H2/O2/NH3  mixtures  as  function  of  Initial  NH3 
mole  fraction  and  laser  energy  deposition  density  for  several  fuel 
equivalence  ratios:  (A)  0.35,  (7)  0.5,  (0)  1.0,  (•)  2.0,  (  )  3.0. 

The  error  bar  given  for  the  cases  with  fuel  equivalence  ratio  of  1.0 
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Indicates  one  standard  deviation  in  the  spread  of  the  measured  data 
points.  The  shaded  area  indicates  the  region  of  no  Ignition. 

7.  Minimum  laser  energy  deposition  density  and  minimum  initial  mole 
fraction  of  NH3  needed  for  ignition  of  H2/O2/NH3  mixtures  at  several 
fuel  equivalence  ratios. 

8.  Ignition  delay  time  In  H2/a1r/NH3  stoichiometric  mixtures  as  function 

of  Initial  NH3  mole  fraction  and  laser  energy  deposition  density.  The 

error  bar  Indicates  one  standard  deviation  in  the  spread  of  the 

measured  data  points.  The  shaded  area  Indicates  the  region  of  no 
ignition. 

9.  Ignition  delay  time  in  CH4/O2/NH3  stoichiometric  mixtures  as  function 

of  initial  NH3  mole  fraction  and  laser  energy  deposition  density.  The 

error  bar  indicates  one  standard  deviation  in  the  spread  of  measured 

data  points.  The  shaded  area  indicates  the  region  of  no  ignition. 

10.  Calculated  dependence  of  the  Ignition-delay  time  for  H2/O2  mixtures  as 
a  function  of  fuel -equivalence  ratio,  for  an  Initial  temperature  of 
369  K  and  an  Initial  pressure  of  1  atm.  The  initial  concentrations  of 
the  radicals  NH2,  0,  and  OH  are  taken  to  be  equal,  and  to  be  equal  to 
mole  fraction  0.01,  0.012,  and  0.014  in  three  calculated  cases. 
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FIGURE  1. 
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FIGURE  8. 
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